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Abstract 

We prove global existence of appropriate weak solutions for the compressible 
Navier-Stokes equations for more general stress tensor than those covered by P.- L. 
Lions and E. Feireisl’s theory. More precisely we focus on more general pres¬ 
sure laws which are not thermodynamically stable ; we are also able to handle some 
anisotropy in the viscous stress tensor. To give answers to these two longstanding 
problems, we revisit the classical compactness theory on the density by obtaining 
precise quantitative regularity estimates: This requires a more precise analysis of 
the structure of the equations combined to a novel approach to the compactness 
of the continuity equation. These two cases open the theory to important physical 
applications, for instance to describe solar events (virial pressure law), geophysical 
flows (eddy viscosity) or biological situations (anisotropy). 
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1 Introduction 

The question of global in time existence of solutions to fluid dynamics’ models goes 
back to the pioneering work by J. Leray (1933) where he introduced the concept 
of weak (turbulent) solutions to the Navier-Stokes systems describing the motion of 
an incompressible fluid; this work has become the basis of the underlying math¬ 
ematical theory up to present days. The theory for viscous COMPRESSIBLE fluids 
in a barotropic regime has, in comparison, been developed more recently in the 
monograph by P. L. Lions ^T] (1993-1998), later extended by E. Feireisl and 
collaborators [37! (2001) and has been since then a very active field of study. 

When changes in temperature are not taken into account, the barotropic Navier- 
Stokes system reads 


f ®tP + div(pu) = 0, x 

[ dt(pu) -f- div(pu (*) u) — div o — pf 1 

where p and u denote respectively the density and the velocity field. The stress 
tensor a of a general fluid obeys Stokes’ law cr = S — PId where P is a scalar 
function termed pressure (depending on the density in the compressible barotropic 
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setting or being an unknown in the incompressible setting) and S denotes the viscous 
stress tensor which characterizes the measure of resistance of the fluid to flow. 

Our approach also applies to the Navier-Stokes-Fourier system, as explained in 
section [TOl which is considered more physically relevant. But our main purpose here 
is to explain how the new regularity method that we introduce can be applied to a 
wide range of Navier-Stokes like models and not to focus on a particular system. For 
this reason, we discuss the main features of our new theory on the simpler (11.11) . We 
indicate later in the article after the main ideas how to extend the result to some 
Navier-Stokes-Fourier systems. 

In comparison with Leray’s work on incompressible flows, which is nowadays 
relatively ’’simple” at least from the point of view of modern functional analysis (and 
in the linear viscous stress tensor case), the mathematical theory of weak solutions to 
compressible fluids is quite involved, bearing many common aspects with the theory 
of nonlinear conservation laws. 

Our focus is on the global existence of weak solutions. For this reason we will not 
refer to the important question of existence of strong solutions or the corresponding 
uniqueness issues. 

Several important problems about global existence of weak solutions for com¬ 
pressible flows remain open. We consider in this article the following questions 

• General pressure laws, in particular without any monotonicity assumption; 

• Anisotropy in the viscous stress tensor which is especially important in geo¬ 
physics. 

In the current Lions-Feireisl theory, the pressure law P is often assumed to be of 
the form P(p) = ap 7 but this can be generalized, a typical example being 

P € C 1 ([0, +oo)), P(0) = 0 with 

ap 7-1 — b < P'ip) < -p 7-1 + b with 7 > d/2 ( 1 . 2 ) 

a 

for some constants a > 0, b > 0: See B. Ducomet, E. Feireisl, H. Petzeltova, 
I. Straskraba m or E. Feireisl [34] for slightly more general assumptions. 
However it is always required that P(p) be increasing after a certain critical value 
of p. 

This monotonicity of P is connected to several well known difficulties 

• The monotonicity of the pressure law is required for the stability of the thermo¬ 
dynamical equilibrium. Changes in monotonicity in the pressure are typically 
connected to intricate phase transition problems. 

• At the level of compressible Euler, i.e. when 5 = 0, non-monotone pressure 
laws may lead to a loss of hyperbolicity in the system, possibly leading to 
corrected systems (as by Korteweg in particular). 

In spite of these issues, we are able to show that compressible Navier-Stokes sys¬ 
tems like (HD are globally well posed without monotonicity assumptions on the pres¬ 
sure law ; instead only rough growth estimates are required. This allows to consider 
for the first time several famous physical laws such as modified virial expansions. 
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As for the pressure law, the theory initiated in the multi-dimensional setting by 
P.-L. Lions and E. Feireisl requires that the stress tensor has the very specific 
form 

a = 2 pD(u) + Adivwld — P(p) Id 

with D{u) = (V?i +Vu t )/2, y, and A such that X + 2p/d > 0. The coefficients A and 
p do not need to be constant but require some explicit regularity, see for instance 
[55] for temperature dependent coefficients. 

Unfortunately several physical situations involve some anisotropy in the stress 
tensor; geophysical flows for instance use and require such constitutive laws, see for 
instance [53] and m with eddy viscosity in turbulent flows. 

We present in this article the first results able to handle more general viscous 
stress tensor of the form 


g = A(t)X7u + Adivuld — P(p) Id 

with a d x d symmetric matrix A with regular enough coefficients. The matrix A can 
incorporate anisotropic phenomena in the fluid. Note that our result also applies to 
the case 

g = A(t)D(u) + Adivrt Id — P(p) Id 
where D(u) = (Vu + V?t T )/2 still. 

Our new results therefore significantly expand the reach of the current theory 
for compressible Navier-Stokes and make it more robust with respect to the large 
variety of laws of state and stress tensors that are used. This is achieved through 
a complete revisiting of the classical compactness theory by obtaining quantitative 
regularity estimates. The idea is inspired by estimates obtained for nonlinear conti¬ 
nuity equations in [7] , though with a different method than the one introduced here. 
Those estimates correspond to critical spaces, also developed and used for instance 
in works by J. BOURGAIN, H. Brezis and P. MiRONESCU and by A.C. Ponce, see 
HU and [S3]. 

Because of the weak regularity of the velocity field, the corresponding norm of the 
critical space cannot be propagated. Instead the norm has to be modified by weights 
based on a auxiliary function which solves a kind of dual equation adapted to the 
compressible Navier-Stokes system under consideration. After proving appropriate 
properties of the weights, we can prove compactness on the density. 

The plan of the article is as follows 

• Section [3] presents the classical theory by P. L. Lions and E. Feireisl, with 
the basic energy estimates. It explains why the classical proof of compactness 
does not seem able to handle the more general equations of state that concern 
us here. We also summarize the basic physical discussions on pressure laws and 
stress tensors choices which motivates our study. This section can be skipped 
by readers which are already familiar with the state of the art. 

• In Section [3] we present the equations and the corresponding main results 
concerning global existence of weak solutions for non-monotone pressure law 
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and then for anisotropic viscous stress tensor. Those are given in the barotropic 
setting. 

• Section [4] is devoted to an introduction to our new method. We give our 
quantitative compactness criterion and we show the basic ideas in the simplex 
context of linear uncoupled transport equations and a very rough sketch of 
proof in the compressible Navier-Stokes setting. 

• Section [5] states the stability results which constitute the main contribution of 
the paper. 

• Section ED states technical lemmas which are needed in the main proof and 
are based on classical harmonic analysis tools: Maximal and square functions 
properties, translation of operators. 

• Section [T] and |8] constitutes the heart of the proof. Section [7] is devoted to 
renormalized equation with definitions and properties of the weights. Section 
[8] is devoted to the proof of the stability results of section [5] both concerning 
more general pressure laws and concerning the anisotropic stress tensor. 

• Section ED concerns the construction of the approximate solutions. It uses the 
stability results of section [5] to conclude the proof of the existence theorems of 
section ED 

• The extension for non-monotone pressure laws (with respect to density) to 
the Navier-Stokes-Fourier system is discussed in Section [10] It contains a 
discussion of the state of the art complementing section [2] in that case. It 
follows some steps already included in the book [36] but also ask for a careful 
check of the estimates at each level of approximation. 

• We present in Section [TT| some models occurring in other contexts where the 
new mathematical techniques presented here could be useful in the future. 

• Section m is a list of some of the notations that we use. 

• Section [13] is an appendix recalling basic facts on Besov spaces which are used 
in the article. 

2 Classical theory by E. Feireisl and P.-L. Lions, 
open problems and physical considerations 

We consider for the moment compressible fluid dynamics in a general domain fl 
which can be the whole space K. d , a periodic box or a bounded smooth domain 
with adequate boundary conditions. We do not precise the boundary conditions and 
instead leave those various choices open as they may depend on the problem and we 
want to insist in this section on the common difficulties and approaches. We will 
later present our precise estimates in the periodic setting for simplicity. 
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2.1 A priori estimates 


We collect the main physical a priori estimates for very general barotropic systems 
on R + x 

f d t p + div(pit) = 0, , , 

\ d t (pu) + di v(pu ®u) — T>u + VP(p) = pf, 


where T> is only assumed to be a negative differential operator in divergence form on 
u s.t. 

/ u-Vudx~ — \\7u\ 2 dx, (2-2) 

JQ JQ 

and for any (/> and u 


<f> -Vudx < C ||V0||i2 ||Vu|| L 2. 


(2-3) 


The following estimates form the basis of the classical theory of existence of weak 
solutions and we will use them in our own proof. We only give the formal derivation 
of the estimates at the time being. 

First of all, the total energy of the fluid is dissipated. This energy is the sum of 
the kinetic energy and the potential energy (due to the compressibility) namely 


E(p,u) 




dx, 


where 

e(p) = j P(s)/s 2 ds 

J Pref 

with p re f a constant reference density. Observe that formally from (ED 
Id 2 / ld 2 \ 

d t(p—) + div \ pu 2 ) ~ u ' Vu + u ' VP ( p ) = Pf 'd 

and thus 

4 [ d44 - / u-Vu- f P(p) div u = f pf ■ u. 

az Jn z Jn Jo, Jq 

On the other hand, by the definition of e, the continuity equation on p implies that 


d t {pe(p)) + div (pe(p) u ) + P(p) div it = 0. 


Integrating and combining with the previous equality leads to the energy equality 

4;E{p,u)- [ u-Vu= [ pf-u. (2.4) 

at Jq jq 

Let us quantify further the estimates which follow from (12.411 . Assume that P(p) 
behaves roughly like p 1 in the following weak sense 

C“V - C < P(p) < C p' 1 + C, (2.5) 
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then pe(p ) also behaves like p 1 . Note that (12.51) does not imply any monotonicity 
on P which could keep oscillating. One could also work with even more general 
assumption than m Different exponents 7 on the left hand side and the right- 
hand side for instance... But for simplicity we use m - 

Assuming that / is bounded (or with enough integrability), one now deduces 
from (12. dl) the following uniform bounds 


sup f p\u\ 2 dx < C + E(p°,u°), 
t J n 


sup 

t 




dx < C, 


|Vit | 2 dx < C. 


( 2 . 6 ) 


We can now improve on the integrability of p , as it was first observed by P.-L. Lions. 
Choose any smooth, positive x(^) with compact support, and test the momentum 
equation by X9 = X& P a where B is a linear operator (in x) s.t. 


div g = {p a -p a ) 1 ||V.g||iP < C p \\p a -p a \\ L p, \\B4 >\\lp < C p ||0|| L p, V 1 <p < 00 , 


where we denote by p a the average of p a over f 1. Finding g is straightforward in the 
whole space but more delicate in bounded domain as the right boundary conditions 
must also be imposed. This is where E. Feireisl et al . introduce the BOGOVSKI 
operator. We obtain that 

/ x(t) / p a P(p)dxdt < / \{t) / g (d t (pu) + div (pu ® u) — V u — p f) dx dt 
J Jq J J q 

+ J x(t) J pvPip). 

By (12.51) , the left- hand side dominates 

X ( t ) f p a+1 dxdt . 

J n 

It is possible to bound the terms in the right-hand side. For instance by (12731) 


- J x(t) gPudxdt < C7||Vu|| L 2 ([0i T ],l 2 (n)) ||x Vff||i, 2 ([0> t], L 2 (n)) 

< C || Vu|| L 2 ([0i T]) L 2 (n)) Hx ( P a - p a )||L2 ([0j T]i L 2 (n)) , 

by the choice of g. Given the bound (12.61) on Vu, this term does not pose any 
problem if 2 a < a + 7 . Next 


J Xgdt{pu) dxdt = - J {gx' (i) + X{t) ‘ 
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The first term in the right-hand side is easy to bound; as for the second one, the 
continuity equation implies 


x(t)B(d t (p a - p a )) pudxdt = - J X [£> (div {up a ))\ pu 


(a — 1) B (p a div u — di v(up a ) + (a — l)p a divit)) 


( 2 . 8 ) 


pu. 


Using the properties of B and the energy estimates (12.611 . it is possible to control 
those terms as well as the last one in (12.711 . provided a < 2j/d — 1 and 7 > d/2 which 
leads to 

r T r 


/ / p 7+a dx dt < C(T, E(p°, u 0 )). 

J 0 J n 


(2.9) 


2.2 Heuristic presentation of the method by E. FEIREISL and 
P. L. Lions 

Let us explain, briefly and only heuristically the main steps to prove global exis¬ 
tence of weak solutions in the barotropic case with constant viscosities and power 
7 pressure law. Our purpose is to highlight why a specific form of the pressure or 
of the stress tensor is needed in the classical approaches. We also refer for such a 
general presentation of the theory to the book by A. Novotny and I. Straskraba 
[55] , the monograph Etats de la Recherche edited by D. Bresch |T3; or the book 
by P. Plotnikov and J. Sokolowski [59:. 

Let us first consider the simplest model with constant viscosity coefficients p and 
A. In that case, the compressible Navier-Stokes equation reads on K.+ x fl 

f &tP + div(pu) = 0 , . . 

dt(pu) + div(pu ® u) — pAu — (A + p)S7divu + VP(p) = pf, 

with P(p) = ap 1 . For simplicity, we work in a smooth, Lipschitz regular, bounded 
domain with homogeneous Dirichlet boundary conditions on the velocity 

w| 9 fi = 0. (2.11) 

A key concept for the existence of weak solutions is the notion of renormalized 
solution to the continuity equation, as per the theory for linear transport equations 
by R. J. DiPerna and P.-L. Lions. Assuming p and u are smooth and satisfy the 
continuity equation, for all b £ C([0,+oo)), one may multiply the equation by b'(p) 
to find that (p, u ) also 

dtb(p) + div(b(p)u) + (b'(p)p — b(p))divu = 0. (2.12) 

This leads to the following definition 

Definition 2.1 For any T £ (0,+oo), /, po, mo satisfying some technical assump¬ 
tions, we say that a couple (p, u) is a weak renormalized solution with bounded energy 
if it has the following properties 
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. p 6 L°°(0, T; LT'(fi)) nC°([0, T], L^n)), p > 0 o.e. in (0, T) x to, p | t=0 = p 0 
a.e. in fi; 

• u £ L 2 (0, T; 7Jq(S 1)) ; p|z *| 2 £ L°°(0, T; I/(f2)) ; pu is continuous in time with 
value in the weak topology of L^jJJ^ l+1 \u), (pit)|t=o = mo a.e. tt; 

• (p, u) extended by zero out of solves the mass and momentum equations in 
R d , in V'((0,T) x R d ); 

• For any smooth b with appropriate monotony properties, b(p) solves the renor¬ 
malized Eq. (12.121) . 

• For almost all r £ (0 ,T), ( p,u ) satisfies the energy inequality 

E{ P ,u){t)+ f f (p\Vu\ 

Jo Jn 

In this inequality, E(p,u)(r) = J n (p\u\ 2 /2 + pe(p))(r), with e(p) = f P{s)/s 2 ds 
(p re f being any constant reference density), denotes the total energy at time r and 
E o = f Q |mo| 2 / 2 po + Poo(po) denotes the initial total energy. 

Assuming P(p) = ap 7 (in that case e{p) may equal to ap 7 _ 1 /( 7 — 1)), the theory de¬ 
veloped by P. L. Lions to prove the global existence of renormalized weak solution 
with bounded energy asks for some limitation on the adiabatic constant 7 namely 
7 > 3 d/(d+ 2). E. Feireisl et al. have generalized this approach in order to cover 
the range 7 > 3/2 in dimension 3 and more generally 7 > d/2 where d is the space 
dimension. 

We present the initial proof due to P. -L. Lions and indicate quickly at the end 
how it was improved by E. Feireisl et al.. The method relies on the construction 
of a sequence of approximate solution, derivation of a priori estimates and passage 
to the limit which requires delicate compactness estimate. We skip for the time 
being the construction of such an approximate sequence, see for instance the book 
by A. Novotny and I. straskraba for details. 

The approximate sequence, denoted by (pk, Uk), should satisfy the energy in¬ 
equality leading to a first uniform a priori bound , using that p > 0 and A+ 2 p/d > 0 

sup f (p k |it/c | 2 /2 + ap]./{'y - 1)) dx + p f f |Vufc | 2 dxdt < C, 
t Jfl Jo Jft 

for some constant independent of n. 

For 7 > d/2, we also have the final a priori estimate (12.91) explained in the 
previous subsection namely 

J™ J^pl +a <C(i?,T) for a < |y - 1. 

When needed for clarification, we denote by U the weak limit of a general sequence 
Uk (up to a subsequence). Using the energy estimate and the extra integrability 


2 + (X + p) |divu| 2 ) < E 0 + [ [ pf-u. 

Jo Jn 
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property proved on the density, and by extracting subsequences, one obtains the 
following convergence 

p fe -pinC°([0,T];ii; eak (fi)), 
pi -^inL^ +o ^(( 0 ,T) x Q), 
p k u k - pu in C°([0,T]-L^/^ +1 \n)), 

Pku\v? k —*• pu % v? in V'((0,T) x fi) for i,j = 1,2,3. 


The convergence of the nonlinear terms p k u k and p k u k ® u k uses the compactness 
in time of p k deduced from the uniform estimate on dtp k given by the continuity 
equation and the compactness in time of v fphU k deduced form the uniform estimate 
on dt(p k u k ) given by the momentum equation. This is combined with the L 2 estimate 
on Vrifc. 

Consequently, the extensions by zero to (0,T) x R 3 /£l of the functions p,u,~jf f, 
denoted again p,u,~jfi, satisfy the system in R + x U. d 

( dtp + div (pu) = 0, _ (2 13 n 

( d t (pu) + div(pu ® u) — pAu — (A + /x)Vdivu + aVp 7 = pf. 


The difficulty consists in proving that ( p , u) is a renormalized weak solution with 
bounded energy and the main point is showing that pi = p 7 a.e. in (0, T) x Q. 

This requires compactness on the density sequence which cannot follow from the 
previous a priori estimates only. Instead P. L. Lions uses a weak compactness of 
the sequence {ap k — A + 2pdivu k } k eN* which is usually called the viscous effective 
flux. This property was previously identified in one space dimension by D. Hoff 
and D. Serre. More precisely, we have the following property for all function 
b £ C 1 ([0, Too)) satisfying some increasing properties at infinity 

lim / / (ap)) — (2/x T A)di vu k )b(p k )ipdxdt 

fc -H- oo J Q J n 

= f f (api — (2/i T X)d\\u)b(p)pdxdt (2-14) 

Jo Jn 

where the over-line quantities design the weak limit of the corresponding quantities 
and ip £ T>((0,T) x fl). Note that such a property is reminiscent of compensated 
compactness as the weak limit of a product is shown to be the product of the weak 
limits. In particular the previous property implies that 


pdivu — pdivrt 


P(P)P ~ P(P) P 

A 


Taking the divergence of the momentum equation, we get the relation 


(2.15) 


A[(2/x T A)di vu k - P(p k )] = di v[d t (p k u k ) + di v(p k u k ® u k )]. 


Note that here the form of the stress tensor has been strongly used. From this 
identity, P.-L. Lions proved the property (12.1411 based on harmonic analysis due 
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to R. CoiFMAN and Y. Meyer (regularity properties of commutators) and takes 
the observations by D. Serre made in the one-dimensional case into account. The 
proof by E. Feireisl is based on div-curl Lemma introduced by F. Murat and 
L. Tartar. 

To simplify the remaining calculations, we assume 7 > 3 d/{d + 2) and in that 
case due to the extra integrability on the density, we get that p k £ L 2 ((0,T) x fl). 
This lets us choose b(s) = s log s in the renormalized formulation for p k and p and 
take the difference of the two equations. Then pass to the limit n — > +00 and use the 
identity of weak compactness on the effective flux to replace terms with divergence 
of velocity by terms with density using (12.151) . leading to 



Observe that the monotonicity of the pressure P(p) = ap 7 implies that 


p(p)p- p {p)p < 0. 


This is the one point where the monotonicity assumption is used. It allows to show 
that the defect measure for the sequence of density satisfies 


dft [p k -p](t)= / plogp{t) - plogp(t) dx < dit[p k - p\(t = 0). 


a 


On the other hand, the strict-convexity of the function s 1—► slogs, s > 0 implies 
that dft[/9fc — p] > 0. If initially this quantity vanishes, it then vanishes at every later 
time. 

Finally the commutation of the weak convergence with a strictly convex function 
yields the strong convergence of the density pk in L\ oc . Combined with the uniform 
bound of pk in L 7+a ((0,T) x SI), we get the strong convergence of the pressure term 

pI- 

This concludes the proof in the case 7 > 3d/(d + 2). The proof of E. Feireisl 
works even if the density is not a priori square integrable. For that E. Feireisl 
observes that it is possible to control the amplitude of the possible oscillations on 
the density in a norm L p with p > 2 allowing to use an effective flux property with 
some truncature. Namely he introduced the following oscillation measure 


osc P [p k - p]= sup [limsup \\T n (p k ) - T n (p) ||LP (( 0 ,T )x n], 


n> 1 k —)-+oo 

where T n are cut-off functions defined as 


T n (z) = nT(-), n> 1 

with T £ C 2 {R) 


T(z ) = 2 for z < 1, T(z) = 2 for z > 3, T concave on R. 
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The existence result can then obtained up to 7 > d/2 : See again the review by 

A. Novotny and I. Straskraba [55] . 

To the author’s knowledge there exists few extension of the previous study to 
more general pressure laws or more general stress tensor. Concerning a generaliza¬ 
tion of the pressure law, as explained in the introduction there exist the works by 

B. Ducomet, E. Feireisl, H. Petzeltova, I. Straskraba m and E. Feireisl 
[34] where the hypothesis imposed on the pressure P imply that 

P(z) = r 3 (z ) - r 4 (z) 

where r 3 is non-decreasing in [0,+oo) with r 4 £ C 2 ([0,+00)) satisfying r 4 > 0 and 
r 4 (z) = 0 when z > Z for a certain Z > 0. The form is used to show that it is possible 
respectively to continue to control the amplitude of the oscillations osc p [pk — p\ and 
then to show that the defect measure vanishes if initially it vanishes. The two papers 
[34] and [3D] we refer to allow to consider for instance the two important cases : Van 
der Waals equation of state and some cold nuclear equations of state with finite 
number of monomial (see the subsection on the physical discussion). 


2.3 The limitations of the Lions-Feireisl theory 

The previous heuristical part makes explicit the difficulty in extending the global 
existence result for more general non-monotone pressure law or for non-isotropic 
stress tensor. First of all the key point in the previous approach was 

p(p)p- p(p)p< 0. 


This property is intimately connected to the monotonicity of P(p) or of P(p) for 
p > p c with truncation operators as in [34] or [301 . Non-monotone pressure terms 
cannot satisfy such an inequality and are therefore completely outside the current 
theory. 

The difficulty with anisotropic stress tensor is that we are losing the other key 
relation in the previous proof namely (12.1511 . For non-isotropic stress tensor with an 
additional vertical component and power pressure law for instance, we get instead 
the following relation 


pdivw — pdivu < a 


pKp 1 - pA^pZ 
Px + A 


with some non-local anisotropic operator A M = (A — (p z — p ix )d‘/)~ l d‘i where A is 
the total Laplacian in terms of (x, z) with variables x = (xi, • • • , Xd-i), 2: = Xd- 

Unfortunately, we are again losing the structure and in particular the sign of the 
right-hand side, as observed in particular in m • Furthermore even small anisotropic 
perturbations of an isotropic stress tensor cannot be controlled in terms of the defect 
measure introduced by E. Feireisl and collaborators: Remark the non-local behav¬ 
ior in the right-hand side due to the term A For this reason, the anisotropic case 
seems to fall completely out the theory developed by P.-L. Lions and E. Feireisl. 

Those two open questions are the main objective of this monograph. 
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2.4 Physical discussions on pressure laws and stress tensors 

The derivation of the compressible Navier-Stokes system from first principles is 
delicate and goes well beyond the scope of this manuscript. In several respects the 
system is only an approximation and this should be kept in mind in any discussion 
of the precise form of the equations which should allow for some uncertainty. 

2.4.1 Equations of state 

It is in general a non straightforward question to decide what kind of pressure law 
should be used depending on the many possible applications: mixtures of fluids, 
solids, and even the interior of stars. Among possible equation of state, one can find 
several well known laws such as Dalton’s law of partial pressures (1801), ideal gas 
law (Clapeyron 1834), Van der Waals equation of state (Van De Waals 1873), virial 
equation of state (H. Hamerlingh Onnes 1901). 

In general the pressure law P(p,'d) can depend on both the density p and the 
temperature d. We present the temperature dependent system in full later but 
include the temperature already in the present discussion to emphasize its relevance 
and importance. 

Let us give some important examples of equations of state 

• State equations are barotropic if P(p) depends only on the density. As ex¬ 
plained in the book by E. Feireisl [35] (see pages 8-10 and 13-15), the simplest 
example of a barotropic flow is an isothermal flow where the temperature is 
assumed to be constant. If both conduction of heat and its generation by dissi¬ 
pation of mechanical energy can be neglected then the temperature is uniquely 
determined as a function of the density (if initially the entropy is constant) 
yielding a barotropic state equation for the pressure P(p) = ap 7 with a > 0 
and 7 = (R + c v )/c v > 1. Another barotropic flow was discussed in [50] . 

• The classical Van der Waals equation reads 

(P + ap 2 )(b- p) =cp’d , 

where a , b, c are constants. The pressure law is non-monotone if the temper¬ 
ature is below a critical value, $ < i? c , but it satisfies (11.21) . In compressible 
fluid dynamics, the Van des Waals equation of state is sometimes simplified by 
neglecting specific volume changes and becomes 

{P + a) (b- p) = c pfl, 

with similar properties. 

• Using finite-temperature Hartree-Fock theory, it is possible to obtain a tem¬ 
perature dependent equation of state of the following form 

P(p , $) = a 3 (l + a)p 2+a - a 0 p 2 + kti ^ B n p n , (2-16) 

n> 1 

where k is the Boltzmann’s constant, and where the last expansion (a simplified 
virial series) converges rapidly because of the rapid decrease of the B n . 
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• Equations of state can include other physical mechanism. A good example 
is found in the article where radiation comes into play: a photon assem¬ 
bly is superimposed to the nuclear matter background. If this radiation is in 
quasi-local thermodynamical equilibrium with the (nuclear) fluid, the resulting 
mixture nucleons+photons can be described by a one-fluid heat-conducting 
Navier-Stokes system, provided one adds to the equation of state a Stefan- 
Boltzmann contribution of black-body type 

PrW = a with a > 0, 

and provided one adds a corresponding contribution to the energy equation. 
The corresponding models are more complex and do not satisfy CE3 in general. 

• In the context of the previous example, a further simplification can be intro¬ 
duced leading to the so-called Eddington’s standard model. This approxima¬ 
tion assumes that the ratio between the total pressure P(p,'d) = Pg(Pj^) + 
Pr($) and the radiative pressure Pr($) is a pure constant 

_ W'1 _ = 

Pg(p,#) + Pr(0) P ’ 

where 0 < /? < 1 and Pq is given for instance by (12.161) . Although crude, this 
model is in good agreement with more sophisticated models, in particular for 
the sun. 

One case where this model leads to a pressure law satisfying m is when one 
keeps only the low order term into the virial expansion. Suppose that er = 1 
and let us plug the expression of the two pressure laws in this relation, 

2 a 3 p 3 - a 0 p 2 + kBrfp = ^ ^ 

By solving this algebraic equation to leading order (high temperature), one 
gets 

leading to the pressure law 

T}( Q\ ^ a 3 3 2,7 TD /6a3(l — /3) N. 7/4 

PyPifi) = —P - a 0 p +kB i(-—- )p‘ l , 

which satisfies m because of the constant coefficients. 

However in this approximation, only the higher order terms were kept. Consid¬ 
ering non-constant coefficient or keeping all the virial sum in the pressure law 
was out of the scope of m and leads to precisely the type of non- monotone 
pressure laws that we consider in the present work. 
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• The virial equation of state for heat conducting Navier-Stokes equations can 
be derived from statistical physics and reads 

P(p^) = pV(j2 B nWp n ) 

n> 0 

with B 0 = cte and the coefficients B n ({f) have to be specified for n > 1. 

While the full virial pressure law is beyond the scope of this article, we can 
handle truncated virial with appropriate assumptions or pressure laws of the 
type P{p, d) = P e (p) + P th (p, d) 

• Pressure laws can also incorporate many other type of phenomena. Compress¬ 
ible fluids may include or model biological agents which have their own type of 
interactions. In addition, as explained later, our techniques also apply to other 
types of “momentum” equations. The range of possible pressure laws is then 
even wider. 

Based on these examples, the possibilities of pressure laws are many. Most are 
not monotone and several do not satisfy m, proving the need for a theory able to 
handle all sort of behaviors. 

2.4.2 Stress tensors 

One finds a similar variety of stress tensors as for pressure laws. We recall that we 
denote D(u) = (V« + (Vu) T )/2. 

• The isotropic stress tensor with constant coefficients 

V = p A u + (A + p) V div u, 

which is the classical example that can be handled by the Lions-Feireisl 
theory: See for instance E5U, [55] and [59] with 7 > d/2. See also the recent 
interesting work by P.I. Plotnikov and W. Weigant (see [5B]) in the two- 
dimensional in space case with 7 = 1. 

• Isotropic stress tensors with non constant coefficients better represent the 
physics of the fluid however. Those coefficients can be temperature d dependent 

T> = 2 div (/z(d) D(u)) + V (A(d)divu). 

Provided adequate non-degeneracy conditions are made on p and A, this case 
can still be efficiently treated by the Lions-Feireisl theory under some as¬ 
sumptions on the pressure law. See for instance 1551 or [561 . 

• The coefficients of the isotropic stress tensors may also depend on the density 

V = 2 div {p{p, d) D(u)) + V (A (p, d) div u). 

This is a very difficult problem in general. The almost only successful insight 
in this case can be found in 

[m null [Ml HE] with no dependency with respect 
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to the temperature. Those articles require a very special form of p(p) and A (p) 
and without such precise assumptions, almost nothing is known. Note also the 
very nice paper concerning global existence of strong solutions in two-dimension 
by A. Kazhikhov and V.A. Vaigant where p is constant but A = with 
/3 > 3, see [35] . 

• Geophysical flow cannot in general be assumed to be isotropic but instead 
some directions have different behaviors; this can be due to gravity in large 
scale fluids for instance. A nice example is found in the Handbook written by 
R. Temam and M. ZiANE, where the eddy-viscous term T> is given by 

V = p,hA x u + p z d^u + (A + p) Vdiv u, 

with ph 7^ Pz ■ While such an anisotropy only requires minor modifications for 
the incompressible Navier-Stokes system, it is not compatible with the Lions- 
Feireisl approach, see for instance PS- 


3 Equations and main results: The barotropic case 

We will from now on work on the torus n d . This is only for simplicity in order to 
avoid discussing boundary conditions or the behavior at infinity. The proofs would 
easily extend to other cases as mentioned at the end of the paper. 

3.1 Statements of the results: Theorem 13.11 and 13.21 

We present in this section our main existence results. As usual for global existence 
of weak solutions to nonlinear PDEs, one has to prove stability estimates for se¬ 
quences of approximate solutions and construct such approximate sequences. The 
main contribution in this paper and the major part of the proofs concern the stability 
procedure and more precisely the compactness of the density. 

I) Isotropic compressible Navier Stokes equations with general pressure. 

Let us consider the isotropic compressible Navier-Stokes equations 

f d t p + div(pu) = 0, , , 

\ dt{pu ) + di v(pu (g) u) — pAu — (A + p)Vdivu + VP(p) = p/, 

with 2 p/d+ A, a pressure law P which is continuous on [0, +oo), P locally Lipschitz 
on (0, +00) with P(0) = 0 such that there exists C > 0 with 

C" V -C< P(p) <Cp^ +C (3.2) 

and for all s > 0 

|P'(a)| < iV" 1 . (3.3) 


One then has global existence 
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Theorem 3.1 Assume that the initial data uq and po > 0 with f nd Po = -Mo > 0 
satisfies the bound 

Eo= f _|_ p 0 e(po)) dx < +oo. 

Jn d Z P o 

Let the pressure law P satisfies m and USD with 

7 > (max(2, 7 ) + l) (3.4) 

Then there exists a global weak solution of the compressible Navier-Stokes system 
ED in the sense of Definition \2.1[ 

Remark. Let us note that the solution satisfies the explicit regularity estimate 

J^ d ^Pk{x, t )> r )^Pk{y, t )>ri^h{x — y)x{^Pk){t) < ?? l/2 | log| e / 2 ’ 

for some 9 > 0 where Kh is defined in proposition 14.11 5pk and x are defined in 
Section El see (18.11) . 

II) A non-isotropic compressible Navier-Stokes equations. We consider an 
example of non-isotropic compressible Navier-Stokes equations 

f d t p + div(pu) = 0 , , s 

( dt(pu) + div(pu 0 u) — div (A(t) V-u) — (p + A)Vdivu + VP(p) = 0, ' 

with Aft) a given smooth and symmetric matrix, satisfying 

Aft) = /rid + 8A{t), p > 0, — p + A — ||<5^4(i)||> 0. (3.6) 

d 

where 5A will be a perturbation around p Id. We again take P locally Lipschitz on 
[ 0 , + 00 ) with P( 0 ) = 0 but require it to be monotone after a certain point 

C~ l p 7-1 — C < P'{p) < C p~<- 1 + C. (3.7) 

with 7 > d/2. The second main result that we obtain is 

Theorem 3.2 Assume that the initial data uq and po > 0 with f nd po = Mq > 0 
satisfies the bound 

E 0 = f (1 iE/fhlL + p 0 e(p 0 )) dx < + 00 . 

Jn d %Po 

Let the pressure P satisfy (ED with 
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There exists a universal constant C* > 0 such that if 

IIMHoo <C*(2/x + A), 

then there exists a global weak solution of the compressible Navier-Stokes equation in 
the sense of Definition \2.1\ revlacina the isotropic energy inequality by the following 
anisotropic energy 

E(p, u)(t) + f f (\7 x u T Aft) Wu + (p + A) |divw| 2 ) < E 0 . 

Jo Jn 

Remark. Let us note that the constraint on 7 corresponds to the constraint on p: 
p >7 + 7/(7 — 1 ) where p is the extra integrability property on p. 

3.2 Important comments/comparison with previous results 

The choice was made to focus on explaining the new method instead of trying to 
write results as general as possible but at the cost of further burdening the proofs. 
For this reason, Theorems EH and 13.21 are only two examples of what can be done. 

We explain how to apply our new method to the Navier-Stokes-Fourier system 
(with an additional equation for temperature) in section 1101 The Navier-Stokes- 
Fourier system is physically more relevant than the barotropic case and as seen from 
the discussion in subsection m it exhibits even more examples of non-monotone 
pressure laws. 

I) Possible extensions. In section QTJ we also present applications to various other 
important models, in particular in the Bio-Sciences where the range of possible 
pressure laws (or what plays their role such as chemical attraction/repulsion) is 
wide. But there are many other possible extensions; for instance (E3 could be 
replaced with a more general 

C" 1 p 11 - C < P(p) <Cp^ 2 + C, 

with different exponents 71 7 ^ 72 - While the proofs would essentially remain the 
same, the assumption (E3) would then have to be replaced and would involve 71 
and 72 . Similarly, it is possible to consider spatially dependent stress tensor A(t,x) 
in Theorem l3.2l This introduces additional terms in the proof but those can easily be 
handled as long as A is smooth by classical methods for pseudo differential operators. 

II) Comparison with previous results. 

II-l) Non-monotone pressure laws. Theorem 13.11 is the first result to allow for com¬ 
pletely non-monotone pressure laws. Among many important previous contributions, 
we refer to [3011311501 ED and [551 155 ;. 1551 !2T for the Navier-Stokes-Fourier system, 
which are our main point of comparison. All of those require P' > 0 after a certain 
point and in fact typically a condition like (13.71) . The removal of the key assumption 
of monotonicity has important consequences: 
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• From the physical and modeling point of view, it opens the possibility of work¬ 
ing with a wider range of equations of state as discussed in subsection 12.41 
and it makes the current theory on viscous, compressible fluids more robust to 
perturbation of the model. 

• Changes of monotonicity in P can create and develop oscillations in the density 
p (because some “regions” of large density become locally attractive). It was 
a major question whether such oscillations remain under control at least over 
bounded time intervals. This shows that the stability for bounded times is very 
different from uniform in time stability as t —> +oo. Only the latter requires 
assumptions of thermodynamical nature such as the monotonicity of P. 

• Obviously well posedness for non monotone P could not be obtained as is done 
here for the compressible Euler system. As can be seen from the proofs, the 
viscous stress tensor in the compressible Navier-Stokes system has precisely 
the critical scaling to control the oscillations created by the non-monotonicity. 
This implies for instance that in phase transition phenomena, the transition 
occurs smoothly precisely at the scale of the viscosity. 

• Our results could have further consequences for instance to show convergence 
of numerical schemes (or for other approximate systems). Typical numerical 
schemes for compressible Navier-Stokes raise issues of oscillations in the den¬ 
sity which are reminiscent of the ones faced in this article. The question of 
convergence of numerical schemes to compressible Navier-Stokes is an impor¬ 
tant and delicate subject in its own, going well beyond the scope of this short 
comment. We refer for instance to the works by R. Eymard, T. Gallouet, 
R. Herbin, J.-C. Latche and T.K. Karper: See for instance [55] 32 for 
the simpler Stokes case, [31]|38]|40] for Navier-Stokes and more recently to the 
work [221 TRJ . 

Concerning the requirement on the growth of the pressure at oo, that is on the 
coefficient 7 in (lT4lt 

• In the typical case where 7 = 7, (EH 1 ) leads to the same constraint as in P.- 
L. Lions m for a similar reason: The need to have p £ L 2 to make sense of 
pdivu. It is worse than the 7 > d/2 required for instance in 1351 . In 3d, we 
hence need 7 > 9/5 versus only 7 > 3/2 in [35] . 

• It may be possible to improve on (13.41) while still using the method introduced 
here but propagating compactness on appropriate truncation of p; for instance 
by writing an equivalent of Lemma 17.11 on cj)(p{x)) — 4>(p(y )) as in the multi¬ 
dimensional setting by E. Feireisl. This possibility was left to future works. 
Note that the requirement on 7 > d/2 comes from the need to gain integrability 
as per (12.91) along the strategy presented in subsection 12.11 Our new method 
still relies on this estimate and therefore has no hope, on its own, to improve 
on the condition 7 > d/2. 

• In the context of general pressure laws, and even more so for Navier-Stokes- 
Fourier, assumption (13.41) is not a strong limitation. Virial-type pressure laws, 
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where P(p) is a polynomial expansion, automatically satisfy it for instance as 
do many other examples discussed in subsection 12.41 


II-2) Anisotropic stress tensor. Theorem 13.21 is so far the only result of global exis¬ 
tence of weak solutions which is able to handle anisotropy in the stress tensor. It 
applies for instance to eddy-viscous tensor mentioned above for geophysical flows 

T> = phA x u + p z d z u + (A + p) Vdiv u, 
where ph ^ p z and corresponding to 

Aij = ph S Z j for i, j = 1, 2, A 33 = p z , Aij = 0 otherwise. (3.8) 

This satisfies the assumptions of Theorem 13.21 provided \ph — p x \ is not too large, 
which is usually the case in the context of geophysical flows. 

Additional examples of applications are given in section [TT] but we wish to em¬ 
phasize here that it is also possible to have a fully symmetric anisotropy, namely 
div (A Du) with D(u) = Viz + Viz T in the momentum equation. This is the equiv¬ 
alent of the anisotropic case in linear elasticity and it is also an important case for 
compressible fluids. Note that it leads to a different form of the stress tensor. With 
the above choice of A, Eq. 13.81 one would instead obtain 

div (A Du) = phA-xU + p z d z u + p z Vckiz 2 + (A + p) Vdiv it. 

Accordingly we choose Theorem 13.21 with the non-symmetric anisotropy div (A Viz) 
as it corresponds to the eddy-viscous term by R. Temam and M. Ziane mentioned 
above. But the extension to the symmetric anisotropy is possible although it in¬ 
troduces some minor complications. For instance one cannot simply obtain div it 
by solving a scalar elliptic system but one has to solve a vector valued one instead; 
please see the remark just after (15.1111 and at the end of the proof of Theorem l3.2l in 
section [21 

Ideally one would like to obtain an equivalent of Theorem 13.21 assuming only 
uniform elliptic bounds on A(t) and much lower bound on 7 . Theorem 13.21 is a first 
attempt in that direction, which can hopefully later be improved. 

However the reach of Theorem [X2] should not be minimized because non isotropy 
in the stress tensor appears to be a level of difficulty above even non-monotone 
pressure laws. Losing the pointwise relation between div zi and P(p{x)) is a major 
hurdle , as it can also be seen from the proofs later in the article. Instead one has to 
work with 

divzi(f, x) = P(p(t,x)) + LP{p) + effective pressure, 

with L a non local operator of order 0. The difficulty is to control appropriately this 
non local term so that its contribution can eventually be bounded by the dissipation 
due to the local pressure term. 

Notations. For simplicity, in the rest of the article, C will denote a numerical 
constant whose value may change from line to line. It may depend on some uniform 
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estimates on the sequences of functions considered (as per bounds (15.511 or (15.41) for 
instance) but it will never depend on the sequence under consideration (denoted with 
index k) or the scaling parameters h or Hq. 


4 Sketch of the new compactness method 

The standard compactness criteria used in the compressible Navier-Stokes frame¬ 
work is the Aubin- Lions-Simon Lemma to get compactness on the terms pu and 
pu ® u. A more complex trick is used to get the strong convergence of the density. 
More precisely it combines extra integrability estimates on the density and the effec¬ 
tive flux property (a kind of weak compactness) and then a convexity-monotonicity 
tool to conclude. 

Here we present a tool which will be the cornerstone in our study to prove com¬ 
pactness on the density and which will be appropriate to cover more general equation 
of state or stress tensor form. 

In order to give the main idea of the method, we present it first in this section 
for the well known case of linear transport equations, i.e. assuming that u is given. 
We then give a rough sketch of the main ideas we will use in the rest of the article. 
This presents the steps we will follow for proofs in the more general setting. 

4.1 The compactness criterion 

We start by a well known result providing compactness of a sequence 

Proposition 4.1 Let pk be a sequence uniformly bounded in some L p ((0,T) x n rf ) 
with 1 < p < oo. Assume that ICh is a sequence of positive, bounded functions s.t. 


i. V ?7 > 0, sup I ICh(x) dx < oo, suppKh £ B(0,R), 



If dtPk £ L q ([ 0, T) x W 1,9 (n d )) with q> 1 uniformly in k and 





then pk is compact in L p ([ 0, T] x n d ). Conversely if pu is compact in L p ([ 0, T] x n d ) 
then the above quantity converges to 0 with h. 

For the reader’s convenience, we just quickly recall that the compactness in space is 
connected to the classical approximation by convolution. Denote ICh the normalized 
kernel 



22 






Write 


\\pk-K-h*xPk\\ p L p < nir .ip [ ([ IC h (x- y)\p k (t,x)-p k (t,y)\dx) dy 

\\^-h\\ L i Jn dK Jn d ' 

— TmF~I1— / £h(x-y)\pk(t,x)~ p k (t,y)\ p dxdy 

which converges to zero as h , — > 0 uniformly in k by assumption. On the other- 
hand for a fixed h, the sequence Kh* x Uk in k is compact in x. This complete the 
compactness in space. Concerning the compactness in time, we just have to couple 
everything and use the uniform bound on dt.pk as per the usual Aubin-Lions-Simon 
L emma. 


The K-h functions. In this paper we choose 

*‘ w = (inriiiF' forH£1/2 ' 


with some a > d and Kf t non negative, independent of h for \x\ > 2/3, Kh constant 
outside .8(0,3/4) and periodized such as to belong in C°°(II d \ 8(0, 3/4)). For 
convenience we denote 


Kh{x) 


K h {x) 

WKhWl i' 


This is enough for linear transport equations but for compressible Navier-Stokes 
we also need 


/C/,o (a;) 



We will see that 


(j-i (j 

{h 0 + |a;|) d “ ICho ^ ^ “ (h 0 + |a;|) d ' 


4.2 Compactness for linear transport equation 

Consider a sequence of solutions pk, on the torus n d (so as to avoid any discussion 
of boundary conditions or behavior at infinity) to 

d t pk + div {p k u k ) = 0, (4.1) 


where Uk is assumed to satisfy for some 1 < p < oo 


sup||« fc || L p wr i.p < oo, 

k * 


(4.2) 


with div Uk compact in x, i.e. 


lirnsup e k {h) = / / 

k W^iiWl 1 Jo J n 2d 


K h {x—y) |div Uk{t, a:) —div Uk{t, y)\ p —> 0, (4.3) 
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(4.4) 


as h —> 0. The condition on the divergence is replaced by bounds on p k 
77 < inf p k < sup p k <C, V t £ [0, T], 

X x 


One then has the well known 


Proposition 4.2 Assume p k solves (14.11) with the bounds (14.21) . (14.31) and (14.41) : 
assume moreover that the initial data is compact. Then p k is locally compact and 
more precisely 



V ) | Pk(t,x) - Pk{t,y) | dxdy <C 


\\K h \W 

|log(fr + e fe (/i) +£ fe (/i))|’ 


where 

h{h) = I, y I, / K h (x - y) \p k {x) - p° k (y)\ dxdy. 


This type of results for non Lipschitz velocity fields u k was first obtained by R. J. 
Di Perna and P. L. Lions in [53] with the introduction of renormalized solutions 
for u k £ TT 1,1 and appropriate bounds on div u k - This was extended to u k £ BV, 
first by F. Bouchut in [9] in the kinetic context (see also M. Hauray in [42]) and 
then by L. Ambrosio in [4] in the most general case. We also refer to C. Le Bris 
and P.-L. Lions in 53 m], and to the nice lecture notes written by C. De Lellis 
in [2B]. In general u k £ BV is the optimal regularity as shown by N. Depauw in 
m- This can only be improved with specific additional structure, such as provided 
by low dimension, see [3 HH MM 1ST] . Hamiltonian properties [19] [44] , or as a 
singular integral m- 

Of more specific interest for us are the results which do not require bounds on 
divrtfc (which are not available for compressible Navier-Stokes) but replace them by 
bounds on p k , such as (14.41) . The compactness in Prop. I4.2l was first obtained in [5], 

Explicit regularity estimates of p k have first been derived by G. Crippa and 
C. De Lellis in ]55j (see also [45] for the W 1,1 case). These are based on explicit 
control on the characteristics. While it is quite convenient to work on the charac¬ 
teristics in many settings, this is not the case here, in particular due to the coupling 
between div u k and p(p k )- 

In many respect the proof of Prop. 14.21 is an equivalent of the method of 
G. Crippa and C. De Lellis in [25] at the PDE level, instead of the ODE level. Its 
interest will be manifest later in the article when dealing with the full Navier-Stokes 
system. The idea of controlling the compactness of solutions to transport equations 
through estimates such as provided by Prop. I4.1l was first introduced in [7] but relied 
on a very different method. 


Proof One does not try to propagate directly 



V) I Pk(t,x) - p k (t,y)\dxdy. 
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Instead one introduces the weight w(t,x ) solution to the auxiliary equations 

d t w + Uk ■ Vui = -AMtu| t= o = 1, (4.5) 

where M f denotes the maximal function of / (recalled in SectionEI and A is constant 
to be chosen large enough. 

First step: Propagation of a weighted regularity. We propagate 


R(t) = [ K h (x 
Jn 2d 


y) I Pk(t,x) - pk(t,y)\w(t,x)w(t,y)dxdy. 


Using urn we obtain that 


dt\Pk{x) - Pk(y) I + div x ( u k {x ) | p k {x) - pk{y) |) + div y ( u k {y) \p k {x) - p k (y )|) 

< i(di vu k {x) + div u k {y)) \p k {x) - p k (y) \ 

- i(div Uk(x) - di vuk(y)) ( Pk{x ) + p k {y)) s k , 

where s k = sign (p k {x) — p k {y))- We refer to subsection [7J] for the details of this cal¬ 
culation, which is rigorously justified for a fixed k through the theory of renormalized 
solutions in J2BJ- From this equation on | p k (x) — Pk{y) |, we deduce 


d_ 

dt 


R(t) = / VK h (x-y)(u k (x)-u k (y))\pk(t,x)-Pk{t,y)\w(t,x)w(t,y) 
J u 2d 


/ n 2d 


K h (x - y) (div u k {x) - di vu k (y)) ( p k {x) - p k (y)) s k w(x)w(y) 


+ J ^ K h {x - y) \pk{x) - p k {y)\ (d t w + u k ■ X7w + i divu k vj'j w(y) 

+ symmetric. 

Observe that by ( 14.511 , w < 1 and therefore by ( 14.411 and the definition of e k in ( 14 . 01 ) , 
the second term in the right-hand side is easily bounded 


/ K h (x - y) (div u k {x) - di vu k (y)) (, p k {x ) - p k {y)) s k w(x)w(y) < Ce k {h). 
J n 2d 

For the first term, one uses the well known inequality (see |611 62] or section ED 
\uk(x) -U k (y) I < Cd\x-y\ (M\Vu k \{x) + M \Vu k \(y)), 
combined with the remark that from the choice of Kh 


\S7K h (x - y )| \x — y\ < CK h (x - y ). 

Therefore 



S7K h {x - y) ( u k (x) 


u k (y)) | Pk{t, x) - p k (t, y)| w{t, x) w(t, y) 


<C [ K h (x - y)(M\Vu k \{x) + M\Vu k \(y))\pk(x) - pk(y)\w(x)w(y), 
J n 2d 
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and combining everything 
d 


dt 


R(t) < C \\K h \\ L i e k (h) 


+ / K h \p k (x) - pk{y)\ {d t w + u k ■ Vto + C(dxvu k + M\S7u k \)w) w(y) 

Jn 2d 

+ symmetric. 

Since div u k < d |Viifc| < dM |Vufc|, by taking the constant A large enough in (14.51) . 
dtw + Uk • Vw + C (divUk + M |Vufc|) w < 0, 

and hence 

R(t) < [ K h \p° k (x)-p° k (y)\+Ct\\K h \\ L iE k (h) <C\\K h \\ L i(£ k (h)+E k (h)). (4.6) 
J n 2d 

Second step: property of the weight. We need to control the measure of the set where 
the weight w is small. Obviously if w were to vanish everywhere then the control of 
R(t) would be trivial but of very little interest. From Eq. (14.51) 


d t (p k | logw|) + div (p k \ logw|) = A p k M |Vti fc |. 


And thus 


4 / | logw| p k dx = A f p k M\Vu k \dx <\\\p k \\ LP * \\M\\7u k \\\ L p <C, (4.7) 

at J n d J n d 

by ( 14 . 21 ) . ( 14 . 41 ) and the fact that the maximal function is bounded on L p for p > 1 . 
Third step: Conclusion of the proof. Using again m 

C f c 

\{x, w(t,x) <r]}\< - - - \\ogw\p k dx < - - 

I log 7 ? I Jn d I log??| 


Thus 


/ U 2d 


K h (x - y) | p k (t, x) - p k (t , y)\ 


= K h \p k (t,x) - p k (t,y)\ 

J w(x)>r], w(y)>r] 

/ K h (x-y)\p k (t,x)~ p k (t,y)\, 

' w(x)<r] or w(y)<.r] 


and so 


[ K h (x — y) \p k (t, x) — p k (t, y)\ < 4 f K h \p k (t, x) - p k (t, y)\w(t, x) w(t, y) 

Jn 2d T J 


C 


<C\\K h \\ L i 


£k(h) + £ k (h) 


L 1 


+ 


l lo S??l 

which by minimizing in r) finishes the proof. 


□ 
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4.3 A rough sketch of the extension to compressible Navier— 
Stokes 

The aim of this subsection is to provide a rough idea of how to extend the previous 
method. We will only consider the case of general pressure laws and assume that 
the stress tensor is isotropic. When now considering the compressible Navier-Stokes 
system, the divergence divit*, is not given anymore but has to be calculated from p k 
through a relation of the kind 


di vu k = P{pk) + Rk, (4.8) 

where R k includes the force applied on the fluid, the effective stress tensor... To keep 
things as simple as possible here, we temporarily assume that 

sup ||iifc||L~ < 00, 

k 

t (4 

limsupe fc (/i) =-— 4 — [ f K h (x - y)\R k {t,x) - R k {t,y )\ p — > 0. 
k Jo Jn 2d 

We do not assume monotonicity on the pressure P but simply the control 

\P'(p)\ < Cp<~\ (4.10) 

A modification of the previous proof then yields 
Proposition 4.3 Assume p k solves (14.11) and the bounds 

sup||pfc|| L ?= L i < oo, sup ||pfe||z.f < °o with p> j+1. 

k k * ,x 

Assume moreover that sup fe ||ujfe||L|jji < oo, that SHU holds with the bounds SHD 
on R k and (14.101) on P. Then p k is locally compact away from vacuum and more 
precisely 


JPk(x)>V, Pk(v)>V 

where 


K h (x—y) \p k (t,x) - pk(t,y) \ dxdy < C, 


\Ki 


h\\Li 


'' 'log(/i+ (£k(h)) 


£k(h) = 


\\K h \\» Jn 2d 


K h {x - y) | pl(x) - pl(y)\ dxdy. 


Unfortunately Prop. I4.3l is only a rough and unsatisfactory attempt for the following 
reasons 

• The main problem with Prop. 14.31 is that it does not imply compactness on the 
sequence p k because it only controls oscillations of p k for large enough values 
but we do not have any lower bounds on p k . In fact not only can p k vanish 
but for weak solutions, vacuum could even form: That is there may be a set of 
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non vanishing measure where pk = 0. This comes from the fact that the proof 
only give an estimate on 



y) I Pk(t,x) - pk(t,y )| w(x)w{y)dxdy, 


but since there is no lower bound on pk anymore, estimates like 1X7]) only 
control the set where w(x)w(y ) is small and both pk(x) and Pk(y) are small. 
Unfortunately \pk(t,x) — Pk(t,y)\ could be large while only one of pk{x) and 
Pk(y) is small (and hence w(x)w(y) is small as well). 

The solution is to work with w(x) + w(y) instead of w(x) w(y). Now the sum 
w(x) + w{y) can only be small if \pk(t,x) — pk(t,y)\ is small as well, meaning 
that a bound on 



K h {x - y) | pk{t,x) - pk{t,y)\ (w(x) + w{y)) dxdy, 


together with estimates like (SZ 3 ) would control the compactness on pk . Un¬ 
fortunately this leads to various additional difficulties because some terms are 
now not localized at the right point. For instance one has problems estimating 
the commutator term in S7Kh ■ (uk(x) — Uk{y)) or one cannot directly control 
terms like div'Ufc(ir) w(y) by the penalization which would now be of the form 
M |Vu*;|(a;) w(x). Some of these problems are solved by using more elaborate 
harmonic analysis tools, others require a more precise analysis of the structure 
of the equations. Those difficulties are even magnified for anisotropic stress 
tensor which add even trickier non-local terms. 

The integrability assumption on pk, p > 7 + 1 is not very realistic and too 
demanding. If p = 7 (l+ 2 /<j) — 1 as for the compressible Navier-Stokes equations 
with power law P(p) = ap 7 , then this requires 7 > d. Improving it creates 
important difficulty in the interaction with the penalization. It forces us to 
modify the penalization and prevent us from getting an inequality like (14.71) 
and in fact only modified inequalities can be obtained, of the type 


sup / | log w(t, x)\ e Pk{t, x) dx < 00 . 

k J n d 


The bounds (14.91) that we have assumed for simplicity on Rk cannot be deduced 
from the equations. The effective pressure is not bounded in L°° and it is not a 
priori compact (it will only be so at the very end as a consequence of pk being 
compact). Instead we will have to establish regularity bounds on the effective 
pressure when integrated against specific test functions; but in a manner more 
precise than the existing Lions-Feireisl theory, see Lemma T8.31 later. 

This is of course only a stability result, in order to get existence one has to 
work with appropriate approximate system. This will be the subject of Section 

m 


Proof One now works with a different equation for the weight 

d t w + Uk ■ Vw = -A (M |Vufc| + pi ), to|t= 0 = 1, (4.11) 

where M f is again the maximal function of /. 


First step: Propagation of some weighted regularity The beginning of the first step 
essentially remains the same as in the proof of Prop. 14.11 One propagates 


R(t)= / K h (x - y) \pk(t,x) - pk{t,y)\w(t,x)w(t,y)dxdy. 

Jn 2d 

The initial calculations are nearly identical. The only difference is that we do not 
have (14.311 any more so we simply keep the term with div Uk ( x ) — div u k ( y ) for the 
time being. We thus obtain 

d f 

-t: R(t) < / K h (x - y) (divu k (x) - divu k (y)) 
at j U 2 d 

( Pk(t , x) + pk(t, y)) s k w(t , x) w(t, y) (4.12) 

-A / (pl(x) + pl(y))\pk(t,x) - pk(t,y)\w(t,x)w(t,y), 

J n 2d 

by taking the additional term in Eq. (14.1111 into account. This is of course where 
the coupling between Uk and pk comes into play, here only through the simplified 
equation (I4.8D . Thus 


hi 2d 


K h (x - y) (div Uk(x) - di vu k (y)) ( Pk(t,x ) + p k (t,y)) s k w(t,x)w(t, y ) 


ln 2d 


K h (x - y) (P(p k (x)) - P(pk(y))) 


(4-13) 


(Pk(t, x) + p k (t, y)) s k w(t, x ) w(t, y) 

+ / K h (x-y)(Rk(x)-Rk(y))(pk(t,x) + pk(t,y))skw(t,x)w(t,y). 
Jn 2d 

By the uniform L p bound on p k and the estimate (14.91) . one has 

n Kh(x - y) ( Rk(x ) - Rk(y)) ( Pk(t,x ) + Pk(t,y )) s k w(t,x) w(t,y) 

[2d 

<C\\K h \\ Ll (e fc (/i)) 1 - 1/ P. 

Using now (EH, it is possible to bound 


(4-14) 


I P(Pk(x)) - P(pk(y ))| < | Pk(t,x) - p k (t,y )| [ \P'(sp k (x) + (1 - s)p k (y))\ds 

Jo 

< c (pl~ 1 (x) + pl~ 1 (y)) \pk(t,x)~ Pk(t,y)\, 
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leading to 



y) ( P(pk{x )) - P{pk{y ))) {pk(t, x) + p k (t, y)) s fc w(t, x) w(t, y) 


<C [ K h (x — y) (pi l {x) + pl 1 (y))(p k (t,x) + p k {t,y)) 

J n 2d 

I Pk(t, x) - p k (t, y )| w(t, x) w(t, y) 

<C / K h (x - y) (pl(x) + pl{y))\p k (t,x) - p k (t,y)\w(t,x)w(t,y). 
J n 2d 


Using now this estimate, the Equality (14.131) . the compactness (14.141) and by taking 
A large enough one finds from (14.121) 

R(t) < R( 0) + C \\K h \\ L i ( s k (h )) 1 - 1 / p . 


Second step and third steps: Property of the weight and conclusion. The starting 
point is again the same and gives 


iu d 


| log w(t, x)| pk{t, x) dx < C, 


with C independent of k but where we now need pk G L p with p > 7 + 1 because of 
the additional term in Eq. (14.111) . Now 


'pk(x)>V, p k (v)>v 


K h (x - y) | Pk(t,x) - Pk(t,y )| 


JPk(x)>ri, Pk(y)>ri, w(x)>r]\ w(y)>t]' 

+ W K \\l' f (1 + Pk{x)). 

Jpk(x)>V, w(x)<r/' 


K h {x - y) I Pk(t,x) - Pk(t,y)\ 


On the one hand, 


f (1 + Pk(x)) < (- + 1 ) f p k (x) 

JPk( x )>V, w(x)<r)' Jw(x)<rj' 


A .x 1 
< c—f 1)7-:— 

77 I log lp 


- f\logw{x)\p k {x)<(- + l)- 
\ J V 


C 


log 77 '| 


On the other hand 


/ K h (x - y) \pk(t,x) - pk(t,y)\ 

’ Pk(x)>r], p k (y)>ri, w(x)>T]', w(y)>r/' 

-7712 l Kh{x — y) \pk(t, x) — Pk{t, y)\ w(t, x) w(t, y). 
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Therefore 


[ K h (x - y)\p k (t,x) - p k (t,y)\ < (-T 1) . ° 

lp*{x)>n, ph(v)>v I lo S ? 7'l 

C 


+ II^IU 1 (s k (h) + (e k (h)y-^) 


which concludes the proof by optimizing in rf. 


□ 


5 Stability results 

5.1 The equations 

5.1.1 General pressure law. 

Let ( p k ,u k ) solve 


d t pk + div (p k u k ) = a k Ap k . 


and 


(5.1) 


A i k (t,x) div u k = P k (t,x,p k ) + A div (d t (p k u k ) + div (p k u k <g» u k )) + F k . (5.2) 

Important remark. Note that we allow here a possible explicit dependence on t and 
x in P k . This does not really affect our stability results and it will allow later to use 
directly the stability estimates for the case with temperature. We will write some 
comments for reader’s who are only interested by the hypothesis mentioned in the 
barotropic theorem. 

The viscosity of the fluid is assumed to be bounded from below and above 


3 P, - < Pk(t;X) < p 

fl 

We consider the following control on the density (for p > 1) 

SU P [IIPfellx,°°([o, t], l 1 (n d )nLT(n d )) + l|Pfc||L*>([o, r]xn d )] < 00 • 
k 


and the following control for u k 


sup 
k L 


I Pk \^k | 2 || L°° ([0, T], LqiP*)) + ll^ 7 M fe||L 2 (0,T;i 2 (n d )) 


< Too. 


We also need some control on the time derivative of p k u k through 
3p > 1, sup \\d t {p k u k )\\ L 2 W -i.p < oo, 


(5.3) 


(5.4) 


(5.5) 


(5.6) 


and on the time derivative of p k namely 


3q > 1, sup \\d t pk\\ L i w -i,i < oo. 
fc * x 


(5.7) 
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Remark. Note that usually (see for instance [35 -[5TH. 15.61 and 15.71 are consequences 
of the momentum equation and the mass equation using the uniform estimates given 
by the energy estimates and the extra integrability on the density. 


Concerning the equation of state, we will consider that for every y 1 x, Pkit , x, s ) 
continuous in s on [0, +oo) and positive, Pt(s = 0) = 0, Pk locally Lipschitz in s on 
(0, +oo) with one of the two following cases: 

• i) Pressure laws with a quasi-monotone property: 

There exists , P, po independent of t, x such that: 

If s > po, P k (t,x,s) is a function P k (t,x) plus a function independent of t, x 
and 

d s (\P k (t,x, s) - P k (t,x)]/s) > 0 for all s > p 0 , lim P k (t,x,s) = + oo, 

L J s—>•+OO 

I Pk{t,x,r) - P k (t,y,s )I 

< P|r - s| + Q k (t, x, y) for all r, s < p 0 , x, y G IT*, 

lim sup sup [ I ^ (| P k (t,x)~ P k (t,y)\ 

h^O k Jo J n 2d W^hUL 1 V 

+ Qk(t, x, y)^j dx dy dt = 0. 

(5.8) 


• ii) Non-monotone pressure laws (with very general Lipschitz pressure laws); 
There exists P > 0, 7 > 0 and P k (t, x) in L 2 ([ 0 , T] x II d ), Qk and R k in 
L 1 ([0, T] x n 2d ) such that for all t, x, y 


| P k (t,x,p k (t,x)) - P k (t,y,p k (t,y ))| < Q k (t,x,y) 


P ({pk(t,x)y 1 + (p fc (f,x)) 7 : ) +P k (t,x) +P k (t,y) \pk(t,x) - p k (t,y)\, 


Pk(t,x,p k (t,x)) < P{pk(t,x)y +R k {t,x) with Rk > 0, 

sup(||Pfc|| i 2([ 0j Tixn^ + ll-RfclliiQo, T]xn d )) < ofo 
k 

lim sup sup [ [ ^ (| P k {t, x) - P k (t, y)\ 

h-> 0 k Jo Jn 2d \\ K h\\L l v 

+ I R k (t,x) - R k (t,y) I + Q k (t,x,yfj dxdydt = 0. 

(5.9) 


Two important remarks. 

f - It is important to note that the general hypothesis on pressure laws will be used 
to conclude in the heat-conducting compressible Navier-Stokes case: see Section [TUI 
The pressure law will include for instance a radiative part (namely a part depending 
only on the temperature as in [36l 1 and a pressure law in density with coefficients 
depending on temperature (see the comments in the section under consideration). 
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Note that compactness in space for the temperature will be obtained directly from 
the thermal flux: This property will be strongly used to check hypothesis for the 
pressure under consideration especially the last one 


lim sup sup 

h—tO k 




K h {x - y) 
\\Kh\\v 


(| P k {t,x) - P k (t,y )| 


+ | R k (t, x ) - R k (t, y )| + Q k (t, x, y)') = 0. 


(5.10) 


Of course, the first step will be to show that with the pressure under consideration 
estimates on the velocity and the density mentioned in the stability process are 
obtained. 

2 - In the basic case where P k does not depend explicitly on t or x (namely P k = 0, 
Qk = 0 and R k = 0) then (15.91) reduces to the very simple condition 

\Pk{r) - P k {s)\ < Pr^~ l \r — s|. 

Note that this assumption is satisfied if P is locally Lipschitz on (0, +oo) with 

\P'(s)\ < Ps ^- 1 

namely with the hypothesis mentioned in Theorem 3.1. Readers who are interested 
by the barotropic case with this hypothesis are invited to choose P k = 0, Q k = 0 
and Rk = 0 in the proof. Hypothesis on the pressure law are used in the subsection 
named ” The coupling with the pressure law”. 

Remark. Note that i) with the lower bound P{p) > C~ 1 p y — C provides the same 
assumptions than in the article [31] by E. Feireisl. Point i) will be used to construct 
approximate solutions in the non-monotone case. 


5.1.2 A non-isotropic stress tensor. 

In that case, assume ( p k ,u k ) solve (15.11) with a k = 0 

d t pk + div (p k u k ) = 0. 

and 

div u k = v k P k {p k ) + v k a^Afj, P k {p k ) 

+ v k (A- a M E k )~ l div (d t (p k u k ) + div ( p k u k ® u k )) (5.11) 
where A M = (A - Efc) -1 E k . 

Remark. If one considers a symmetric anisotropy, div (A Du) in Theorem 13.21 then 
instead of (15.11[) we have the more complicated formula 

div u k =v k P k (p k ) + v k a^A^ P k (p k ) 

i (O. aA) 

+ v k div (A J - a M E k )~ ( d t (p k u k ) + div (p k u k O u k )), 
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where A ^ = (A I — a^Ek)~ l ■ Ek- But now Ek and Ek may be different and are 
vector-valued operators, so that in particular (A I — a^Ek) -1 means inverting a 
vector valued elliptic system. Except for the formulation there would however be no 
actual difference in the rest of the proof. 

Coming back to ( 15 . 111 ) . we assume ellipticity on vk 


(5.13) 


0 < v_< Vk < v < oo. 


We assume that Ek is a given operator (differential or integral) s.t. 

• (A — Efc ) -1 A is bounded on every L p space, 

• A^ = (A — .Efc ) -1 Ek is bounded of norm less than 1 on every L p space and 

can be represented by a convolution with a singular integral denoted by A^ 


still 



Note here that to make more apparent the smallness of the non isotropic part, we 
explicitly scale it with a M . We consider again the control (15.41) on the density but for 
p > y 2 /(y — 1), and the bound (15.51) for Uk■ We also need the same controls: (15.61) 
on the time derivative of pkUk and (ED on the time of the pk- 

The main idea in this part is to investigate the compactness for an anisotropic 
viscous stress obtained as the perturbation of the usual isotropic viscous stress tensor, 
namely —div {A Vu) + (A + p) Vdivu assuming A = p Id + 5A and = ||<5A|| < e 
for some small enough e. 

5.2 The main stability results: Theorems 15.11 T5.2I and 15.31 
5.2.1 General pressure laws 

The main step in that case is to prove the two compactness results 

Theorem 5.1 Assume that pk solves (O), Uk solves ED with the bounds ED, 
(15.51) . (15.61) . (15.71) . and that pk and Fk are compact in L l . Moreover 

i) If ak > 0, we assume the estimate ED on pk with 7 > 3/2 and p > 2 and 
quasi-monotonicity on Pk through ED- 

ii) if au = 0, it is enough to assume ED with 7 > 3/2 and p > max( 2 , 7 ) and 
only (15.91) on Pk- 

Then the sequence pk is compact in L] oc . 

We also provide a complementary result which is a more precise rate of compact¬ 
ness away from the vacuum namely 
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Theorem 5.2 Assume again that pk solves m with ak = 0, Uk solves m with 
the bounds <53*, (ESD. m, m and that p k and F k are compact in L l . Assume 
that ( 15 . 41 ) holds with 7 > d/2 and p > max( 2 , 7 ) and that Pk satisfies ( 15 . 91 ) . Then 
there exists 6 > 0 and a continuous function e with e( 0 ) = 0 , depending only on p k 
and Fk s.t. 



^■Pk{x)>r] I pk(y)>v k-h(x 


y)x{5pk) < 


C\\K h \\ L i 

t ? 1 / 2 | log(e(/i) + h e )\ 9 / 2 ' 


For instance if Pk, Rk, Pk and Fk are uniformly in W S J for s > 0, then for some 
constant C > 0 


J^^Pp^^Pkly)^ R h(x y)x(dpk) < C ^ 1/2 | Yogh\ 9 / 2 ' 

Since those results depend on the regularity of pk and Fk, we denote £o{h) a contin¬ 
uous function with e 0 ( 0 ) = 0 s.t. 

f f Kh{x - y) (\F k (t,x) - F k (t,y)\ + \pk(t,x) - pk(t,y)\ 

Jo Jn 2d v 

+ I Pk(t,x) - P k {t,y )| + | R k {t,x) - R k it,y)\ + \Q k {t,x,y)\^j < e 0 {h) ||A^|| l i. 

( 5 . 14 ) 


5.2.2 Non isotropic stress tensor 

In that case our result reads 

Theorem 5.3 Assume that p k solves <53*, Uk solves ( 15 . 111 ) with the bounds ( 15 . 51 ) . 
( 15 . 61 ) . ( 15 . 71 ) and ( 15 . 131 ) together with all the assumptions on E k below ( 15 . 111 ) . Assume 
as well that Pk satisfies ( 15 . 81 ) and that ( 15 . 41 ) with'y > d/2 and p > 7 2 /( 7 ~ 1)- There 
exists a universal constant C* > 0 s.t. if 


then pk is compact in L\ oc 


Remarks. Theorems 15.11 15.21 15.31 are really the main contributions of this arti¬ 
cle. For instance, deducing Theorems 13. II and 13.21 follows usual and straightforward 
approximation procedures. 

As such the main improvements with respect to the existing theory can be seen 
in the fact that point ii) in Theorem 15.II does not require monotonicity on P k and in 
the fact that Theorem 15.31 does not require isotropy on the stress tensor. 

Our starting approximate system involves diffusion, a k 7 ^ 0, in the continuity 
equation & As can be seen from point i) of Theorem 15. II our compactness result 
in that case requires an isotropic stress tensor and a pressure P k which is monotone 
after a certain point by (15.81) . This limitation is the reason why we have to consider 
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also approximations P k and E k of the pressure and the stress tensor. While it may 
superficially appear that we did not improve the existing theory in that case with 
diffusion, we want to point out that 

• We could not have used P.-L. Lions’ approach because this requires strict 
nronotonicity: P k > 0 everywhere. Instead any non-monotone pressure P sat¬ 
isfying (15.91) can be approximated by P k satisfying (15.81) simply by considering 
Pk = P + £k as long as 7 > 7 and thus without changing the requirements 
on 7 . 

• E. Feireisl et al can handle “quasi-monotone” pressure laws satisfying (15.81) 
together with diffusion but they require higher integrability on pk for this: 
p > 4 in (15.41) . This in turn leads to a more complex approximation procedure. 


6 Technical lemmas and renormalized solutions 


6.1 Useful technical lemmas 

Let us recall the well known inequality, which we used in subsection 14.21 and will use 
several times in the following (see [51] for instance) 

|$(a:)-$(l/)|<C'|a:-y|(M|V$|(a:)+M|V$|(i/)), ( 6 . 1 ) 


where M is the localized maximal operator 

Mf{x)= sup * f f(x + z)dz. ( 6 . 2 ) 

r<i |-B( 0 , r)\ J B (0,r) 

As it will be seen later, there is a technical difficulty in the proof, which would lead 
us to try (and fail) to control M\X7u k \(y) by M\Vu k \(x). Instead we have to be more 
precise than m in order to avoid this. To deal with such problems, we use more 
sophisticated tools. First one has 

Lemma 6.1 There exists C > 0 s.t. for any u £ II / 1 ’ 1 (II <i ), one has 


Kz) - u(y )| <C\x-y\ ( D\ x _ y \u{x ) + D\ x _ y \u{y)), 


where we denote 


D h u{x) 


1 I |Vu(x + .-)| ^ 

h J\ z \<h \z\ d 


Proof A full proof of such well known result can for instance be found in [15] in 
a more general setting namely u £ BV. The idea is simply to consider trajectories 
7 (t) from x to y which stays within the ball of diameter \x — y\ to control 

\u(x)-u(y)\< [ 7 '(t) • Vu( 7 (t)) dt. 

Jo 
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And then to average over all such trajectories with length of order \x — y\. Similar 
calculations are also present for instance in [33]. □ 

Note that this result implies the estimate m as 
Lemma 6.2 There exists C > 0, for any u £ W 1,p {Jl d ) with p > 1 

D h u(x) < CM\S7u\{x). 


Proof Do a dyadic decomposition and define ig s.t. 2 1 < h < 2 10 


D h u{x) < j 


| Vu(x + z) | 


i>i 0 


2~ i ~ 1 <\z\<2~ 


Id—1 


dz 


^E 

i>i 0 


2(i+l) (d-1) 

h 


/ 2~ i ~ 1 <\z\<2~ i 

2 _i 
h 


|Vu(a: + z )| dz 


< 2 d ~ l |J3(0, n.)| — M\Vu\(x) < CM\Vu\(x). 

i>io 


□ 


The key improvement in using Dh is that small translations of the operator Dh 
are actually easy to control 

Lemma 6.3 For any 1 < p < oo, there exists C > 0 s.t. for any u £ iZ^lP) 


dh 


lh 0 J n<i 


K h {z) ||-D| Z | u(.) - D \ z | u{. + z)\\ L p dz — <C (6.3) 


where the definition and basic properties of the Besov space are recalled in 

section\T3 _ [ As a consequence 


dh 


/ / Kh{z) ||D| Z | u{.) - D\ z \ u(. + z )\\ L 2 dz— < C\ logft. 0 | ' \\u\\ H i. (6.4) 

Jh 0 Jn d n 

It is also possible to disconnect the shift from the radius in D r u and obtain for 
instance 

fl p 7 7 

/ / Kh{z)K h (w)\\D\ z \u(.) - D\ z \u(. + w )\\ L 2 dz dw < C \ log h 0 \ 1/2 ||u|| ff i. 

JhoJn™ h 

(6.5) 


We can in fact write a more general version of Lemma 16.31 for any kernel 
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Lemma 6.4 For any 1 < p < oo, any family N r £ pF s,1 (n d ) for some s > 0 s.t. 


sup sup r s / \z\ s \N r (z) — N r (z + rijj)\dz < oo, 

M<i r Ju d 

sup (IliVrllii + r S ||-/Vr-||w s ' 1 ) < OO, 

r 

there exists C > 0 s.t. for any u £ L p (Jl d ) 

f [ K h (z)\\N h -ku(.) - N h *u{. +z)\\ L p dz'-^ < C\logh 0 \ 1/2 \\u\\ LP . (6.6) 

Jhn Jn d h 


We will mostly use the specific version in Lemma 16.31 but will need the more general 
Lemma 16.41 to handle the anisotropic case in Lemma 18.61 Both lemmas are in fact a 
corollary of a classical result 

Lemma 6.5 For any 1 < p < oo, any family L r of kernels satisfying for some s > 0 

f L r = 0, sup(||L r || L i +r s H-Lrll^a.i) < C L , supr -5 ( \z\ s \L r (z)\dz < C L . 
Jr r J 

(6.7) 

Then there exists C > 0 depending only on Cl above s.t. for any u £ L p (J \ d ) 


C dr 

/ ||L i .*u|| l ,-<C'|H| b o i . 

Jh 0 r 


j4s a consequence for p < 2 

f l 
'h 0 


r i r j r 

/ || L r * u\\ LP — < C | log/i 0 | 1/2 \\u\\ L p. 

Jh n r 


( 6 . 8 ) 


(6.9) 


Note that by a simple change of variables in r, one has for instance for any fixed 
power l 


[ || L r i * u\\ LP — <Ci | log/i 0 | 1/2 |M|lp. 

Jho r 


Remark. The bounds and m could also be obtained by straightforward 
application of the so-called square function, see the book written by E.M. Stein [6l], 
We instead use Besov spaces as this yields the interesting and optimal inequalities 
|E3D-([62D as an intermediary step. 


Proof of Lemma 16.31 and Lemma 16.41 assuming Lemma 16.51 First of all 
observe that Dh u = Nh * u with 


N h 


1 

h |2| d_1 
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which satisfies all the assumptions of Lemma 16.41 Therefore the proofs of Lemmas 
16.31 and linTl are identical, just by replacing Dh by IV/,*. Hence we only give the proof 
of Lemma 16.31 Calculate 


— dh f 

K h(z) — < / 

, h J h 


1 Ch v ~ d dh 

— < 


C 


Iho n Jho (h + \ z \y h (\z\ + ho) 

Note also for future use that the same calculation provides 


d ' 


dh 


1 


L 0 Kh{z) h -c(\z\ + h 0 y 


( 6 . 10 ) 


Therefore, using spherical coordinates 

r l r 


I h 0 Jn d 

< 


K h (z) ||.D|z| u(.) - D |z| u(. + z)\\ LP dzdh 

C [ [ || D r u(.) — D r u(. + r uj)\\lp - —duj. 

Js d -iJh n r + h o 


Denote 


L„{x) = - ir^yfry, L UtJ .(x) = r d L u (x/r ), 


and remark that £ W Stl with a norm uniform in u> and with support in the unit 
ball. Moreover 

D r u(x) — D r u(x + rco) = J |Vit|(:r — r z) L UJ (z) dz = L Utr * |Vit|. 


We hence apply Lemma 16.51 since the family L u<r satisfies the required hypothesis 
and we get 

dr 


'ho 


\Luj,r *V u\\lp — < ^ II^IIbJ p » 


with a constant C independent of oj and so 


I h 0 J n d 

< 


K h (z) ||H| 2 | u(.) - D\,\ u{. + z)\\ L p dzdh 

C [ [ \\L u r -k Wu\\ L p — du < C [ |M| b i dw, 

Js*-lJhn r JS d ~l 1,P 


yielding m - The bound (16.41) is deduced in the same manner. The proof of the 
bound (16.51) follows the same steps; the only difference is that the average over the 
sphere is replaced by a smoother integration against the weight l/(l + |w>|) a . □ 


Proof of Lemma 16.51 First remark that L r is not smooth enough to be used as 
the basic kernels in the classical Littlewood-Paley decomposition (see section fl3l) 
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as in particular the Fourier transform of L r is not necessarily compactly supported. 
We use instead the Littlewood-Paley decomposition of it. Denote 

Uk = 'bk* u. 


The kernel L r has 0 average and so 

L r *U k = / L r {x-y)(U k {y)-U k {x))dy. 

Jn d 

Therefore 

\\L r *U k \\ LP < [ L r (z)\\U k (.)-U k (. + z)\\ LP dz 
Jn d 

< / L r (z)\z\ s \\U k \\w‘.p dz, 

Jn d 

yielding by the assumption on L r , for k < | log 2 r | 

\\L r *U k \\ LP <Cr s 2 ks \\U k \\ LP , (6.11) 

by Prop. 113.21 Note that C only depends on / |^| s L r (z ) dz. 

We now use similarly that L r £ W Sjl and deduce for k > | log 2 r| by Prop. 113.21 

\\L r * U k \\ L p < \\L r \\ w ..i \\U k \\ w - s , P < Cr~ s 2~ ks \\U k \\ LP , (6.12) 

where C only depends on sup r r s ||i r ||vy«,i. From the decomposition of / 


dr 

r 


r 1 A r 00 r 1 

/ \\L r *u||lp — = E / \\L r -kU k \\LP 
Jho k =0 “ 

00 ( r 2 ~ k dr r 1 

<Cj2\\U k \\L P (n fc <|i og2h0 | / r*2 fcs -+ / 
fc=0 V Jh ° 


r~ s 2~ k 


, dr 


max(fao,2 -fc ) 


by using (16.111) and (16.121) . This shows that 


dr 


' ho 


L r *u\\ L P—<C J2 Wk\\LP+C Yj 


~) — k s 


fe<| l°g2 h o\ 


_ /jS 

k>\log 2 h 0 \ 0 




Lp • 


(6.13) 


Now simply bound 


-ks 


E \\Uk\\L>+ E -rr\\Uk\\LP<CY 2 k \\Uk\\LP 


k <| log 2 h 0 \ 


fc>| log 2 h 0 \ 


k—0 


= C\\u\\ b o a , 


which gives 
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Next remark that 


~TT- W U k\\LP < Csup \\U k \\ L p < C || U || B o i0o . 

fe>|log 2 ft, 0 | 0 

Therefore (16.131) combined with Lemma It .3.31 yields 

/ \\L r -ku\\ LP — < C y/ |log 2 ho\ \\u\\lp + C||u|| B o , 

Jho r 

which gives (16.91) by Prop. 113.2I D 

Finally we emphasize that 
Lemma 6.6 The kernel 

>C h o(z)= 

Jho n 

also satisfies i) and ii ) of Prop. 14.11 

Proof This is a straightforward consequence of using (16.101) . □ 


6.2 A brief presentation of renormalized solutions 

Many steps in our proofs manipulate solutions to the transport equation, either 
under the conservative form 


dtp + div (pu) = 0, 

(6.14) 

or under the advective form 


dtw + u ■ Vw = F. 

(6.15) 

We will also consider the mrticular form of (16.151) 


dtw + u ■ Vu> — f w, 

(6.16) 


which can directly be obtained from (16.15[) by taking F = f w. 

However since u is not Lipschitz, we do not have strong solutions to these equa¬ 
tions and one should in principle be careful with using them. Those manipulations 
can be justified using the theory of renormalized solutions as introduced in [5S]. In¬ 
stead of having to justify every time, we briefly explain in this subsection how one 
may proceed. The reader more familiar with the theory of renormalized solutions 
may safely skip most of the presentation below. 

Assume for the purpose of this subsection that u is a given vector field in 
The basic idea behind the renormalized solution is the commutator estimate 
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Lemma 6.7 Assume that p G L 2 x and w G L 2 X . Consider any convolution kernel 
L G C 1 , compactly supported in some B(0,r) with f nd Ldx = 1. Then 

\\div{L e -k x {pu) — uL e -ku)\\ L i —> 0 as e —> 0, 

\\{L e -k x (u ■ V x w) — u ■ \7 x L e *u>|| i i —> 0 as e —> 0. 


The proof of Lemma fo. 71 is straightforward and can be found in [28]. Note however 
that the techniques we introduce here could also be used, a variant of Prop. 14.21 to 
make the estimates even more explicit. From Lemma 16.71 one may simply prove 

Lemma 6.8 Assume that p G L 2 x is a solution in the sense of distribution to 
(16.14[) . Assume w G L 2 X , with F G L 1 , a solution in the sense of distribution to 
(16.1511 . Then for any x £ VF 1,00 (K), one has in the sense of distribution that 

dtx(p) + div{x{p ) u) = (x(p) - px'(p)) divu, 

9 t x( w ) + u ■ Vx(io) = F x'(w). 

Finally if in addition p G L Pl , w G L p2 , u G L P3 with 1/pi + l/p 2 + 1/P3 < 1 and 
F G L q tx with 1/pi + 1/q < 1, then in the sense of distribution for any x £ kF 1,00 (R) 

dt(px(w )) + div(px{w)u) = Fx'(w) p. 


Of course Lemma 16.81 applies to (16.161) in the exact same manner just replacing F 
by fw, provided that f G L p and w G L p with 1/p* + 1/p = 1 (so that F G L 1 ) 
and f w G L q x . 

Lemma 16.81 can be used to justify most of our manipulations later on. Remark 
that all terms in the equation make sense in V: For instance u ■ S7w = div {u w) — 
w divu which is well defined since u , divrt and w belong to L 2 . The proof of Lemma 
16.81 is essentially found in }28j and consists simply in writing approximate equations 
on L e -k p 1 L e * w, perform the required manipulation on those quantities and then 
simply pass to the limit in e. 

As a straightforward consequence, we can easily obtain uniqueness for (16.141) . 
Consider two solutions pi, p 2 G L 2 x to (16.141) with same initial data. Apply the 
previous lemma to p = p\ — p 2 and \(p) = \p\ and simply integrate the equation 
over n d to find 

d n 

— \pi(t, x) - p 2 (t,x)\dx = 0. 
at Ju d 

Thus 

Lemma 6.9 For a given p° G L 2 , there exists at most one solution p G L 2 x to 

m . 
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The uniqueness for the dual problem (16.15[) or (16.16[) is however more delicate and 
in particular the previous strategy cannot work unless divzt £ Lff. The estimates 
are now slightly different from (16.151) or (16.161) and we present them for (16.161) as we 
use more this form later on. 

If one considers two solutions w\ and u >2 to (16.161) and a solution p to (16.14[) , one 
has 

37 / p(t,x)\w!{t,x) - w 2 (t,x)\dx = [ f p(t,x)\wi{t,x) - w 2 (t,x)\dx, 
at J n d J U d 

leading to 

Lemma 6.10 Assume that 

• p £ Lj x solves (16.141) . 

• w\ and W 2 are two solutions in L/ x to (16.161) with w\(t = 0 ) — W 2 (t = 0 ) for 
a given f £ L/ x and Wi £ L p , i = 1, 2, with 1/p* + 1/p = 1. 

• p £ L Pl , w £ L P2 , u £ L P3 with 1/pi + l/p 2 + l/ps < 1 and fw£ L\ x with 
1/px + 1 /q < 1. 

• Finally either f £ L°° or f < 0 and p > 0. 

Then w\ = W 2 p a.e. 

Of course if p > 0 everywhere then Lemma 16.101 provides the uniqueness of the 
solution to (16.151) . But in general p could vanish on a set of non zero measure (this 
is the difficult vacuum problem for compressible Navier-Stokes). In that case one 
cannot expect in general uniqueness for (16.151) . 

We will use the same strategy of integrating against a solution p to the conser¬ 
vative equation (16.141) to obtain some bounds on logw for w a solution to (16.151) . 


Lemma 6.11 Assume that 

• p > 0 in L^ x solves (16.141) . 

• w\ is a solution in x to (|6.16l) with 0 < w± < 1 , wi(t = 0 ) = w° for a given 
f £ L/ x and w £ L p , i = 1, 2, with 1/p* + 1/p = 1. 

• p £ L Pl , w £ L p2 , u £ L P3 with 1/pi + l/p 2 + 1/P3 < 1 and fw£ L\ x with 

1/pi + 1 /q < 1. 

Then one has for any 0 < 9 < 1 


/ (1 + | log w(t,x)\) e p(t, x) dx < / | log w°\ e p(t,x)dx 

Ju d J n d 

+0 f f \f(s,x)\(l + \logw(s,x)\) 6 ~ 1 p(s,x)dxds. 

Jo Jn d 
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The lemma is proved simply by applying Lemma 16.81 (the last point) to a sequence 
Xe(w) = (1 + | log(e + w)\) e , as for a fixed e > 0, Xe is Lipschitz. One then integrates 
in t and x and finally passes to the limit e —> 0 by the monotone convergence theorem. 

Note that the log transform allows to derive (16.151) from (16.161) but requires in 
addition log w £ L 2 while Lemma 16.111 does not require any a priori estimates on 
log w. 

Let us finish this subsection by briefly mentioning the existence question. This 
does not use renormalized per se, although as we saw using the solutions once they 
are obtained requires the theory. 

For uniqueness, the conservative form was well behaved and the advective form 
delicate. Hence for existence, things are reversed. Unless diva £ L°° it is not 
possible to have a general existence result for (16.141) . In general a solution to (16.141) 
with only divw, £ L 2 may concentrate, forming Dirac masses for instance. 

But it is quite simple to obtain a general existence result for (16.151) 

Lemma 6.12 Assume that w° £ L°°(n d ) and either that f £ L°°(n d ) or that f < 0, 
/ £ L\ x and w° > 0. Then there exists w £ L°°([0, T\ x n d ) for any T > 0 solution 
to (16.151) in the sense of distributions. 


Proof Consider a sequence u n £ C°° s.t. u n converges to u in L 2 Hf. Define the 
solution w n to 

d t w n + u n ■ V x w n = f w n , w n (t = 0) = w°. 

This solution w n is easy to construct by using the characteristics flow based on u n . 
Now if f £ L°° then 

IK(i,.)lkr<lk°IUr e m ^. 

In the other case if w° > 0 then w n > 0. Furthermore if / < 0 then 

IK(t,OIU- <lk°IU-. 

So in both cases w n is uniformly bounded in L°°([0, T] x n d ) for any T > 0. 
Extracting a subsequence, still denoted by w n for simplicity, w n converges to w in 
the weak-* topology of L°°([0, T] x n d ). 

It only remains to pass to the limit in u n ■ V x w n = div ( u n w n ) — w n divu„ which 
follows from the strong convergence in L 2 of u n and divw n . Similarly one may pass 
to the limit in / w n . □ 


7 Renormalized equation and weights 

We explain here the various renormalizations of the transport equation satisfied by 
Pk, We then define the weights we will consider and give their properties. 
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7.1 Renormalized equation 

We explain in this subsection how to obtain the equation satisfied by various quan¬ 
tities that we will need and of the form Z k ’ v x{p k ~ P V k ) where Z k ’ v is chosen as 
Zf’ v = Kh(x - y)W k ’ k with W k ’ y = W kt h(t,x,y) appropriate weights. 

Lemma 7.1 Assume pk solves m with (IQl l for p > 2 and that Uk satisfies (1531) . 
Denote for convenience 

Pk=Pk(t, x), p v k = p(t, y), u% = u k (t,x), u v k = u k (t, y), 
5p k =Pl-p v k , Pk = Pk + Pi- 

Then for any x £ W 2, °° 


[ [K h (x - y)W k $ x($Pk)\ (*) - [ [K h {x - y)W%$ x{$Pk)\ (0) 

Jn 2d Jn 2d 

+ [ f (x'(Spk)Spk-x{fipk))(div x u% +div y u y k )K h (x-y)W k f 

Jo Jn 2d 

- f f x(dPk) [u% ■ V x K h (x -y) + u y k ■ V y K h (x - y) 

Jo Jn 2d 

+ a k ( A x + A y )K h {x - y)]W k $ 

ot k S7I<h(x - y) xifipk) 


- 2 


Jo Jn 2d 

ft 


V W x,v — V W x,v 
v x yy k,h v V yy k,h 


2 cx k 


o Jn 2d 


K h (x - y) x(Spk) [A* W k f + A y W k fi 



n 2d 


x(5pk)[d t W x ; y + u x k -\/. 


W x,v 4- v y 
• yv k,h w a k 


v y w x : y 


- a k (A x + A y )W k ' y \K h (x - y) 

~ \ f [ x'($Pk)Kh{x - y)W x ’ y (div x u% - div y u v k )p k 
1 Jo Jn 2d 

+ \ f f x'{3pk)Kh(x ~ y)W kk (div x u k + div y u k )Sp k 
1 Jo Jn 2d 

-a k [ [ x"(Pk)K h (x~y)W k f(\X/ x Sp k \ 2 + \X/ y Sp k | 2 ). 
Jo Jn 2d 


Proof The result essentially relies on a doubling of variable argument and straight 
forward algebraic calculations. 

Since p k solves m, one has that 

d t p x k + di v x {p%u%) = a k A x p k , 

d t p v k + di v y {p v k u v k ) = a k A y p v k . 
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Recalling 5pk = p% — p k , and using that pk £ Tf x with p > 2 and hence p^div Uk is 
well defined, one can check that 


d t 5pk + div x (u% Spk ) + div v {u v k 5p k ) = a k (A x + A y )Sp k - p v k div x u k + pldiv y u v k . 
Then, recalling the notation p k = p k + pf, we observe that 

- p v k div x u x k + p x k div y u v k = ^ (div x u x k p x k - div x u% pf div y u\ p x k - div y u\p\ 

- div x u% p% - div x u x k p v k + div y u\ p x k + div y u\ p v ^j 
= ^ (div x u x k + div y u v k ) Sp k - ^ (div u x k - div u y k ) p k . 

Consequently, we can write 

d t 5p k + div x (u k 5p k ) + div y (u v k 5p k ) = a k (A x + A y )Sp k 

1 1 (7-1) 

+ -(di v x u x k + di v y u v k ) Sp k - -(di vu x k - divuf) p k , 

We now turn to the renormalized equation that means the equation satisfied by 
x(8 k ) for a nonlinear function s H > x(s). Formally the equation can be obtained 
by multiplying ED by x'{5p k ). If a k = 0 and p k is not smooth then the formal 
calculation can be justified following Di Perna-Lions techniques using regularizing 
by convolution and the estimate (15.51) . i.e. u k £ L^H X . Then 


dtx(Spk) + di v x {u x k x(5p k )) + div y(u v k x(Sp k )) 

+ (x'(Spk)8p k ~ x(8p k )) (di v x u k + di v y u v k ) 

= a k {A x + A y )x(8p k ) - ^ x!{8Pk)(div x u k - div y u y k )p k 
+ \x'{8p k ){div x ul + div y u v k )Sp k - a k x" (pfc)(|V x 8p k \ 2 + |V y dp fc | 2 ). 

For any V k , V k and smooth enough Z k ’ k , one has in the sense of distributions 
Z x : y A x V k x = di v x (Z x JX7xV k x ) - div x {V k A7 x Z k ’ v ) + V£A x Z x k y 




div 


y(Z k .h v y 


VyV y ) 


div y (V k V 


yZ k 'h) + V k A y Z 


%,y 

k,h' 


Consequently, we get the following equation for Z k ’ k x(8p k )' 


dt[Z x k % X {8p k )] + di v x (ulx(8p k )ZZ;V) + di v v (ulx{8p k )Zl%) 

+ ( x'{8p k )Sp k - x(8pk)) (div x u^ + di v y u v k )Z x ' y 
X(8p k )[d t z x : y + U% ■ S7 x z x k % + u\ ■ V y Z x : y - a k (A x + A y )Z^ h ] = r.h.s , 
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with 


r.h.s. = ~x {$Pk)Z%%(div x u x k - di v y u v k )p k + ^x' (^Pk)Z k ’J(dW x u% + di v y u v k )6p k 
- a kX "(p k )Z^ y h {\V x Sp k \ 2 + |V y <5p fc | 2 ) 

+ 2 a k Z k ’J (A x + A y ) X (6p k ) + OL k [div x (Z k ’JV xX (Sp k )) - di v x ( X (Sp k )V x Z k ^) 

+ div y(Z k 'JVyx(Sp k )) - div y (x(Spk)VyZ x ’J)]. 

Integrating in time and space and performing the required integration by parts, we 
get the desired equality writing Z k ’ y = Kh{x — y)W k ' y . □ 


7.2 The weights: Choice and properties 

In this subsection, we choose the PDEs satisfied by the weights, we state and then 
prove some of their properties. 

7.2.1 Basic considerations 

We define weights w(t,x) which satisfy 


dtiv + u k (t,x) ■ X7 x w =—D w + cx k A x w, w(t = 0, x) = w°(x) (7-2) 

for some appropriate penalization D depending on the case under consideration. 
Note that w depends on k. but we do not precise the index to keep notations simple. 
The choice of D will be based on the need to control “bad” terms when looking at 
the propagation of the weighted quantity. The choice will also have to ensure that 
the weights are not too small, too often. 

7.2.2 Isotropic viscosity, general pressure laws. 

The case with a k > 0 and monotone pressure. The simplest choice for the penaliza¬ 
tion D to define wq is 

£>o = AM|V« fc |, (7.3) 

with A a fixed constant (chosen later on) and Ad the localized maximal operator as 
defined by ( 16 . 21 ) . We choose accordingly in that case 

w 0 |t=o = Wo = 1- 


The case a k = 0 and non-monotone pressure. In the absence of diffusion in (15.11) 
(i a k = 0) and when the pressure term P k is non-monotone for instance, one needs to 
add a term in the penalization. This would lead to very strong assumptions, in 
particular on the exponent p in (15.41) (and hence 7 ) as explained after Prop. 14.31 It 


47 


is possible to obtain better results using that p £ L p for some p > 2, by taking the 
more refined 

-y- =pk |divM fe | + |divu fc | + M |Vn fc | + p\ + Pk Pk + Rk, ( 7 -4) 

for the general compactness result. We take for simplicity 

wi| t =o = w° = exp (-Asupp°) . 

The reason for the first term in D\ compared to Dq is to ensure that w\ < e~ Xf>k 
which helps compensates the penalization in pj. to get the property on p^ \ logWi| e 
for some 9 > 0. The three last terms are needed to respectively counterbalance: 
Additional divergence terms in the propagation quantity compared to wo, the same 
M |Vztfc| as for wq and the p'J, for terms coming from the pressure. 

7.2.3 Anisotropic stress tensor. 

The choice for the penalization, denoted D a in this case and leading to the weight 
w a , is now 

^ = M\X7u k \ + K h * (|divufc| + |A M p^|). (7.5) 

Note that the second term in the penalization is used to control the non local part 
of the pressure terms. As initial condition, we choose accordingly 

w a \t=o =w° a = l. 

7.2.4 The forms of the weights. 

Two types of weights W are used 

W(t,x,y) =w(t,x) +w{t,y), or W(t,x,y) = w(t,x)w{t,y). 

The first one will provide compactness and will be used with (17731) or G3D- The 
second, used with (1741) gives better explicit regularity estimates but far from vacuum 
and is considered for the sake of completeness. Therefore one defines 

Wo{t,x,y)=w 0 {t,x) + w 0 {t,y), Wi(t,x,y) = wi(t,x)+wi(t,y), 

W 2 {t,x,y) = wi(t,x)wi(t, y), W a {t,x,y) = w a (t,x) + w a (t,y). 

As for the penalization, we use the notation W when the particular choice is not 
relevant and W», i = 0, 1, 2 or a otherwise. For all choices, one has 

d t W + u k (x) ■ V X W + u k (y) ■ V y W = -Q + a k A x , y W. (7.7) 

The term Q depends on the choices of penalizations and weights with the four pos¬ 
sibilities 

Qo = D 0 w 0 (t,x) + D 0 wo(t,y), Q i = D 1 wi(t,x) + £>i wi(t, y), 

Q 2 = (Di(t,x) + D 1 (t,y))wi(t,x)wi(t,y), Q a = D a w a (t,x) + D a w a (t, y). 
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7.2.5 The weight properties 

We summarize the main estimates on the weights previously defined 

Proposition 7.2 Assume that pk solves (15.ip with the bounds (15.51) on Ukand (15.41) 
with p > max( 2 , 7 ). Assume that Pk, Rk are given by (15.91) and in particular Pk 
is uniformly bounded in L^ x and Rk in L\ x . Then there exists weights wq, w w a 
which satisfy Equation m with initial data respectively 


'u>°|t=o = 1 , wi|t=o = exp(—A sup p 0 k ), w a | t= o = 1 

and Do, D\, D a respectively given by (ESP , El and E3 such that 
i) For any t,x: 


0 < w 0 (t,,x) < 1, 0 < w a {t,x) < 1, 0<wi(t,x)<e ( 7 . 8 ) 


ii) One has 



if oik = 0 and p > max( 2 , 7 ) then similarly there exists 9 > 0 s.t. 



while finally if p > 7 + 1 



(7.9) 


Hi) For any ip we have the two estimates 



and if p > 7 + 1 



iv) Denoting w a ,h = Kh * vu a , if p > 7 , we have for some 0 < 9 < 1 



K h * ( K h * {\divu k \ + lA^p'lDwa)- 

(K h * {\divu k \ + I^Pfcl)) w a ,h Lg dt ~Jf- ^ CI l°g h 0 \ e , 


with q = min( 2 ,p/ 7 ). 
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Remark 7.3 Part i) tells us that Wi is small at the right points (in particular when 
Pk is large). On the other hand we want Wi to be small only on a set of small mass 
otherwise one obviously does not control much. This is the role of part ii). We use 
part iv) to regularize weights in the anisotropic case. Part Hi) is also used to get 
control under the form given in \4.1\ from the estimates with weights. 

Remark 7.4 Even when ctk > 0, it would be possible to define D\ in order to have 
a bound like w\ < e ~ x Pk(t,x) q _ p or instance take 

=p q k ~ 1 (?- 1) \divu k \ + |V log w\\ 2 

+ | divuk | + M | V'ltfc | + L pf + Pk Pk + Rk 1 

with 

(\ a\ \ _ q- 2 

V ! + V fc J (9- 1 ) 

However one needs q < p/2 and q > 2 which already forces p > 4. Moreover the 
main difficulty when ctk > 0 comes from the proof of Lemma \8.1\ which forces us to 
work with I\ ^ * wq and not Wq. Because of that any pointwise inequality between 
w(x) and pk{x) is mostly useless. 


Proof Point i). This point focuses on the construction of the weights satisfying 
the bounds CLHD - 

Construction of wq . The construction of the weights wq is classical since it satisfies 
a parabolic equation; moreover since Dq is positive then one has 0 < wq < 1. 
Construction of w a : Choosing 'tr a |t=o = 1 and / = — D a < 0 and noticing that 
D a G L' 2 and we hence easily construct (see Lemmas 16.12116.8l and l6.91) w a such that 
0 < w a < 1 solution of 

d t w a + Uk ■ S/w a + Da Wa = 0, W a |t =0 = l. 

Note that D a > |divzifc| because M |Vufc| > |Vitfc| > |divufc| and w a also solves 
d t w a + dW(u k w a ) + {D a - divzi fc ) w a = 0, 

so that by the maximum principle, we also have w a < pk where pk > 0. This means 
that we can actually uniquely define w a by imposing w a = 0 if pk = 0. 

Construction of w\: Choosing wi|t = o = exp(—Asupp°) and / = —D± < 0 and 
noticing that D i G L 1 and D\ > |divufc| we again construct w\ just as w a such that 
0 < w\ < 1 with w±{t, x) = 0 where pk{t, x) = 0 and solution to 

d t wi + div(zz fc uii) + (Pi - divztfc) wi = 0, uq| t=0 = exp(-Asupp°). 

Note that using renormalization technique on the mass equation 

d t [exp(—Ap fe ] + div(u fc [exp(-Ap fe )]) + [-Xpkdivu k - divu fc ] exp(-A p k ) = 0. 
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Subtracting the two equations we get the following equation on g = w\ — exp(— Xp k ): 

d t g + div(u fc p) + (D 1 - di vu k )g = -(£>1 + Xp k divu k ) exp(-A p k ) 

Recall now that D\ > |divitfc| and D\ > — Xp k d\vu k , thus using the maximum 
principle, we get 

w\ < e~ xpk 

recalling that we have w\ = 0 where p k = 0. 

Point ii). By point i), w < 1, hence |logic| = — logic and from (17721) . denoting 
| log Wi | = Ai, 


d t (Ai) + u k ■ \7 x (Ai) - Di = -—Awi 

Wi 

= a k AAi - a k \VAi\ 2 . 


(7.10) 


For Aq, we directly apply 

dt{p k Aq) + V, ■ ( p k Aq) = Dq p k + a k A(p k Aq ) — 2 a k X7p k • VHo — a k p k \X7 Aq\ 2 , 
and integrate to find 

J p k A 0 dx=C + j J D 0 (t,x) p k dxds - a k j J p k \S7 A 0 \ 2 dx ds 

— 2a k / V p k • V Aq dx ds. 


Simply bound 


J Vp/c • VAq < — J Pfc|VA 0 | 2 + J 


IVpfcl 2 

Pk 


On the other hand, using renormalization techniques, 


^■ J Pk log p k dx = - J p k div u k dx - a k J ^ Pk ~ dx. 

However as p > 2 then p k div u k is bounded uniformly in L] x and p k log p k in 
This implies, from the previous equality, that 

ai f [)Ail dxds <c, 


0 J Pk 


and consequently 


-2 a k J J Vpfe ■'VA i dxds<^- J J p k \\7 Ai\ 2 dx ds + C. 
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Using this in the equality on J p k Aq given previously, we get 

J p k A 0 dx = C + J j D 0 (t,x)p k dxds-^ J j p k \S7A 0 \ 2 dxds. 

In that case, we know that ||D 0 ||l 2 < C A and since p > 2 we get 

J p k | logicol dx + a k J J p k |VA 0 | 2 dxds < C (1 + A). (7.11) 


Concerning ie a , the estimate is similar, even simpler as a k = 0, to get 

J p k | log w a | dx < C( 1 + A) 


assuming p > 7 + 1. Indeed M\S7u k \ and Kh * |divufc| are bounded in L 2 by (15.51) . 
Finally 


0 J n d 


Pk Kh * (| A^ f? k |) 



7 /( 7 +!) 


0 J n d 


\A,Pl\ 


,71 (7+l)/7 


since A^ is continuous on any L q space for 1 < q < + 00 . The right-hand side is 
bounded assuming p > 7 + 1 . 


For w 1 , the estimate is as before a bit different. We now assume that a k = 0, 
define A\ = (1 + A\) 9 and obtain from (17.101) 


d t A 1 + u k ■ Vli - e ^ 1 + ^y_ e 


Integrating and recalling A\ > A p k , M |Vufe| and P k by (15.91) are uniformly bounded 
in L 2 and R k in L 1 


Pk A\ dx ‘P ( ' (‘ 


1 + Pk 

1 


|divufc| dxds 


+ C [ [ - Pk , . (M |V u k \ + P k p k + R k ) dx ds 

Jo J l + p ( k 9) 

+ c ft , 

Jo J 1 + p k 


for some 9 > 0 depending on p — max(2, 7 ). This gives the desired control regarding 
p\ logwi | 0 for an exponent 6 small enough. 
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Point iii) Estimate (17.111) will not be enough in the proof and we will need to control 
the mass of pk where Kh*wo is small. Denote 

tth,ri = {x £ H d , K h *w 0 (x ) < 77 }, £l h>v = {x £ Slh, v , w 0 {x) > y/rj}. 


The time t is fixed during this argument and for simplicity we omit it. 

One cannot easily estimate | flh,r /1 directly but it is straightforward to bound 
Assume x £ £lh,r] Kh *wq(x) < 77 . From the expression of Kh, if |<5| < h 


Kh(z + 5) < 


\\K h \\~ L l 

(/i+|^ + <5|) q 


<CK h (z), 


we deduce that for any y £ B(x, h) 


Kh * wo{y) < Ci 7 . 

Now cover £l Vt h by C d with C d disjoint hyper-cubes of diameter h/C. For any i, 
denote f V v h = Cl v> h D Cf. 

If Q l h 7 ^ 0 then K h * wo(x) < C i] on the whole C d . In that case 


Cr) h d 


1 * wq > 



[ K h (x 


y) w 0 (y)dydx 


We conclude that |n* h \ < C y/rjh d . Summing over the cubes, we deduce that one 
has 

\& v ,h\ < Cy/rj. 

Finally 



p k dx < 



Pk dx + 


2 

llogr/l 


f Pk | log W 0 \dx < C ?? V 2 -l /27 + 

J n<* I l°g 7 71 


since p £ for some 7 > 1. This is the desired bound. 

The same bound may be obtained on the quantity pk^K hi , Wa < r] in a similar way 
when p > 7 + 1 because of bound ( 17 . 91 ) on pk \ logw a |. 


Point 


iv) Denote to simplify / = |divttfc| + |A M pj, |, then by the definition 
r 1 ft _ _ _ dh 

/ / K h *(K h * fw a ) - (K h * f)w a ,h dt — 

J ho Jo L “ d 

< [ f Kh(z ) (K h */(. + z) - K h */(.)) w Q (. + z) 

J ho JO 

< [ f K h (z) (K h *f(t,.+z)-K h *f{t,.)) 

Jo Jho 

< C | log * 0 | 1/2 I \\f(t,.)\\ L «dt<C\\ogh 0 \ 1/2 , 

Jo 


dh 

Li h 


of Wa,h 


dh 

h 
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by a direct application of Lemma l6.4l with Nf l: = I\h provided / is uniformly bounded 
in L\ L| which is guaranteed by q < min(2, p/ 7 ). 

□ 


8 Proof of Theorems 15.1L 15.21 and 15.3 

We start with the propagation of regularity on the transport equation in terms of 
the regularity of divzifc, more precisely div x it/ — di v y u v k . We prove in the second 
subsection some estimates on the effective pressure. This allows to write a Lemma 
in the third subsection controlling div^u^: — div y u v k and then to close the loop in the 
fourth subsection thus concluding the proof. 


8.1 The propagation of regularity on the transport equation 

This subsection uses only Eq. m on pk without yet specifying the coupling between 
div Uk and pk (for instance through (15.21) 1. Recall that we denote 


Spk = Pk(t,x ) - pk(t,y), 

Choose any C 2 function y s.t. 


Pk = pk(t,x) + p k (t,y). 


x(0 - ^x'(£K 


<2^(0^ x'(£)£<cx( 5 )<^ia 


( 8 . 1 ) 


It is for instance possible to take y(£) = £ 2 for |£| < 1/2 and x(0 = |£l f° r ICI ^ 1- 
Note however that the first inequality on the right forces y(£) > |^|/C. 

Similarly for the anisotropic viscous term, for some £ > 0, choose any y„ S C 1 


\-£ 


Xa( 0 - ^ < 

(r+ 1) (-xait - m - i )+2 xa(e - o) 


X'a(C)^ X , „(0^<C , Xa(0<C'|^| 1+/ s 


( 8 . 2 ) 


1 


> -(r-r)-j-xa(e- 0 (^ + 0 - 


Note that it is possible to simply choose Xa = But to unify the notations and 

the calculations with the other terms involving y, we use the abstract y a . 

The properties on this non-linear function will be strongly used to characterize 
the effect of the pressure law in the contribution of div X iik(x) — di v y Uk(y) in the 
anisotropic case. The form of y a and the choice of £ will have to be determined very 
precisely so that the corresponding which will be exactly counterbalanced by the A 
terms in Lemma 18.21 

We can write two distinct Lemmas respectively concerning the non monotone 
pressure law case and the anisotropic tensor case. 
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Lemma 8.1 Assume that p k solves ED with estimates ED and ED on Uk- 
i) With diffusion, ak > 0, if p > 2, 3Sh 0 {k) —> 0 as k —> oo for a fixed ho 

/* l /* dh 

/ / Kh{x-z)Kh(y-w)Wo(t,z,w)Kh(x-y)x(Spk)dxdydzdw — 

Jh 0 J n 4d h 

< C (e ho (k) + | log /i 0 | 1/2 ) 

1 ft f dh 

-77 / / / Kh(x-y) (divu k (x)-divu k (y))x'oPkW 0 K h (x-z) I\ h (y-w)— 

1 Jho Jo Jn* d ", 

1 ^ '* 


m Kh(x - y) ( divu k (x) + divu k (y)) 

i 4d 


A 

2 


f l /.t 


(x'(dpk) 5p k -2x(Sp k )) W 0 K h (x - z) Kh(y - w )~f 
K h {x - y)x(Spk) K h (x - z) M\Vv,k\(z) w 0 {z) 


Ih 0 Jo Jn 3d 

where we recall that Wq = Wo(t,x) + wo(t,y). 
ii) Without diffusion, a k = 0, if p > 2, 


n i /» dh 

/ / K h (x-y)x(Spk)(wi(t,x)+wi(t,y))dxdy — < C | log/i 0 | 1/2 ||u fc || ff i 

Jh 0 Jn 2d h 


— 2 A 


/ho Jo Jn 2d 

rt /*! 


-Kb (a; - 2/) (pZ(z) + -Pfc(z) Pfc(z) + Rk) wffx) xi^Pk) 


dh 


- 2 


m K h (x - y) ( divu k (x ) - divu k {y )) (\x!($Pk)Pk 
i 2d 

+ x{6pk) - \x!{S Pk) Spk) Wi(x) y. 

For the derivation of explicit regularity estimates, we have as well the version with 
the product weight, namely 


/ n 2d 


K h (x - y)x{Spk) wi(t,x) w± (t, y) dx dy < C 


n K h (x - y) (divu k {x) - divu k (y )) x'(Sp k ) Pk wffx) wffy) 

[2d 

-A f f K h (x — y) (p1(x) + P k (x) p k {x) + R k (x) + p k (y) 

Jo Jn 2d v 

+ Pk{y)Pk(y) + Rk(y)) x^Pk)wi{x)w\(y). 


For convenience, we write separately the result we will use in the anisotropic case 
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Lemma 8.2 Assume that p k solves (ED with estimates (15.5I) - (15.4I) . Without diffu¬ 
sion, a k = 0 , assume (18.21) on \a with p > 7 + £ + 1, and denote w a ,h = Kh*w a , 




K h (x - y) 


lh 0 Jn 2d 
< 


(w a ,h(x) + Wa t h(y)) Xa(Spk)(t) 


Kh{x - y) 


h 0 J n 2d 


( W a ,h(x ) +W a ,h(y)) Xa(,Spk.)\t=0 + C \ log h 0 \ 6 + I + 11 - II a , 


with the dissipation term 


n„ = a 


0 Jh 0 Jn 2d 


w a ,hi x )Xa{Spk) Kh* {\divu k \ + \A^pJ\){x)K h ^-, 


while 


I = -- 


Kh{x - y) 


h 0 Jo Jn 2d 


( divu k (x) - divu k {y))x'a(fiPk) Pk 


(w a ,h(x) + W a ,h(y )), 


and 


1 rt r is 


11 = — 
2 


h 0 Jo J n 2d 


Kh(x - y) 

h 


( divuk(x ) + divuk(y)) 

( Xa( S Pk)SPk - ^Xa(Spk)) (W a ,h{x) + W a>h (y)). 


Remark. We emphasize that the A terms in relations i) and ii) of Lemma ED come 
from the penalization in the definition of the weights wq and w\. They will help 
to counterbalance terms coming from the contribution by &\v X u k (x) — div y u k (y). 
Similarly the nonlocal term II a follows from the definition of the weights w a ■ bad 
terms coming from the viscosity anisotropy. 


Proof 

Case i). Denote 

W 0> h(t,x,y)= / K h (x - z) K h (y - w)W 0 (t, z,w) dzdw, 


/ n 2d 


and let us use \ in the renormalized equation from Lemma 17.11 We get 
/ W 0: h{t,x,y)K h (x - y) x(5pk) dxdy = A + B + D + E 


hi d 



K h (x - y) (divu fc (x) - divu fe (y)) x!(5pk) Pk W 0 ,hdxdy, 

K h (x- y) (div u k {x) +divu fc (y)) {x'{Sp k )6p k - 2x(5p k ))W 0ih dxdy, 
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with, by the symmetry of K k , Kh and Wo,h, and in particular since V x Wo,h(x,y) = 
V v W 0 , h (y,x) 

A = / (uk(t,x) - u k (t,y))-\7K h (x - y) x{5 Pk)W 0 , h dxdydt, 

Jo Jn 2d 

n K h (x - y) (d t W 0 ,h(x) + u k ■ V x W 0:h - a k A x W 0: h) x($ Pk) dxdydt, 

[2d 

! OL k / / x(<5 p k ) [A x K h (x - y) W 0th + 2 K h (x - y) A* W 0 ,h] dx dy dt, 

Jo Jn M 

E = -2a k [ [ K h W 0yh x"(S Pk)\^ x fip k \ 2 dxdy dt. 

Jo Jn 2d 

Now note that by the convexity of y, E < 0. Then simply bound using (18.11) . 

D <8a k h~ 2 H-K^llii ||p||li < C a k h~ 2 ||A"/j|| l i, 
leading us to choose £h 0 (k) = ak h~ 2 


B = 2 


D = 2 < 


As for B , using Eq. ( 17 . 21 ) 

B = Bi - 2 J K h (x - y) x($Pk) K h * XtV R 0 dx dy dt, 

with 

B\ = 2 J K h (x - y) x($Pk) ( u k {x ) - w fc (^)) • V x K h (x - z) 

K h (y — w) Wo(t, z, w) dx dy dz dw dt. 

We recall that Ro(x, y) = Dq(x) wq(x) + Do(y) wo(y) with Dq = XKh * (M |Vufc|) 
and we thus only have to bound B\. By Lemma 16.11 we have 

Bi < C J K h (x - y) x(fipk) ( D\ x _ y \u k {x ) + D\ x _ y \u k {z)) 

\x - z\ \S7K h (x - z)\K h (y - w) W 0 (t,z,w) 

< C J K h (x - y) x(S p k ) ( D\ x _ y \u k {x ) + D\ x _ y \u k (zj) I< h (x - z)K h (y - w) W 0 , 

as |z| \VK h \ < CK h . 

Next recalling that Wo(t, z, w) = wo(t, z) + wo{t, ui), by symmetry 

Bi<C J K h (x - y)x{Spk) D\x- y \u k {z) K h (x - z)I< h (y ~ w) w 0 (t,z) 

+ C J K h (x — y) x(dpk) (D\ x _ y \Uk(x) + D\ x _ y \u k {w) — 2D\ x _ y \u k (z)) 

K h (x - z)K h (y - w) w 0 (t, z). 
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Since D\ x _y\u k {z) < C M\S7uk\(z), for A large enough, the first term may be bounded 
by 

-- J K h (x-y)x(6pk)K h *x (M\X7u k \w 0 )dxdydt. 

Use the uniform bound on ||pfc||LP with p > 2, to find 

J Kh{ x - y ) x($Pk) ( D\ x _ y \u k {x ) + D\ x _ y \u k {w) - 2D\ x _ y \u k {z))K h w 0 (t,z) 

< C J ||.D| r |Ufc(.) + -D| r |Ufc(. + r + u) - 2D\ r \Uk(. + v)|| i2 K h (r) I\ h {u) K h (v), 

where we used that w = x + (y — x) + (w — y) . We now use Lemma 16.31 and more 
precisely the inequality (16.51) to obtain 

/ / K h {x - y) x($Pk) ( D\ x _ y \u k {x ) + D\ x _ y \u k {w) - 2D\ x _ y] u k (z))K h w 0 — 


< C'|log/i 0 | 1/2 / \\u k (t,.)\\ H i dt < C\\ogh 0 \ 1/2 . 
Jo 


/o 

Therefore we have that 
r 1 


' ho 


B I|A '" IIl1 dh < C e ho (k) + C | log /to | 1/2 


3A 


h 0 JO 


Kh(x - y) x{$Pk) M\Vu k \(z) Kh(x - z) wq(z) dzdxdydt 


dh 

h 


The computations are similar for A and we only give the main steps. We have, using 
again Lemma 16. II 


l n 2d 


VK h (x - y) ■ ( u k (t,x ) - u k {t,y))x(5pk) W 0 ,h 


<C K h {x-y)(D\ x _y\Uk{x) + D\ x _ y \u k {y))x{&Pk)Wo, h . 

J n 2d 

By decomposing M-'o,/,,, we can write just as for Bi 

/ VK h {x - y) ■ ( u k (t,x ) - u k (t,y))x(Spk ) W 0 ,h 
Jn 2d 

<C [ K h (x - y) M\Vuk\(z)x(S Pk) K h (x - z)w 0 (z) 

J n 2d 

~b C j K k (x y) (D| a ._y| u k {^x') -{- -D| j; _y| u k {y) T D\ x _y\ a k {w) 3 D |x_y| 

Jn 2d 

x(S p k ) I<h{x - z) K h (y - w) w 0 (z). 

The first term in the r.h.s. can again be bounded by 

fK h {x-y)x(Spk)K h *x (• M\S7u k \w 0 )dxdydt . 
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The second term in the r.h.s. is now integrated in h and controlled as before thanks 
to the bound (1531) in Lemma 16.31 and the uniform L p bound on pk and H 1 on uk- 
This leads to 


[ A -j- < C \ log 1 1/2 

J K h 

A f ^ f dh 

+ T / / / K h {x-y)M\S7u k \(z)x{5p k )K h (x-z)w 0 {z)dxdydzdt—. 

4 Jho Jo Jn 2d h 


Now summing all the contributions we get 



{A + B + D + E)\\K h \\~ L l 


dh 

h 


<Cs ho (k) + C\logh 0 \ 1/2 


A 

2 


' ho J 0 


K h {x - y)x(Spk)I<h * (M|Vu fc | w 0 )■ 


Note that indeed £h 0 {k) —> 0 as k —>■ oo for a fixed h 0 . This concludes the proof in 
that first case. 


Case ii). In this part, we assume a*, = 0. We may not assume that pk is smooth 
anymore. However by [28j since pk £ L 2 and so is Vufc, one may use the renormalized 
relation with (p = X and choose ITfyJ = IT,;. We then can use the identity given in 
Lemma 17T1 Denoting 

*(£) = ( x(0 - \ x'(£) c) , 


we get for i = 1,2 


/ K h (x - y)x(<)Pk) Wi(t, x, y) dxdy < A t + B t + D iy 
Jn 2d 

where by the symmetry in x and y 


/ / (uk(t,x)-u k (t,y))-\7Kh{x-y)x(Spk)Widxdydt 

Jo Jn 2d 

A [ [ K h (x-y)(M\Vuk\{x)wi(x) + M\Vu k \(y)wi(y))x(5pk), 
Jo Jn 2d 


while 


Ao = 


n (u k (t,x) - Uk(t,y )) • V K h {x - y) x($Pk) W 2 dxdydt 

[2d 

A [ [ K h (x-y)(M\Wu k \(x) + M\Wu k \(y))w 1 (x)w 1 {y))x(Spk)- 
Jo Jn 2d 
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And 


B, = 2 


n K h {x - y) ( d t wi(x) + u k (x) ■ V^wi + 2div X u k (x) x{8p k ) w\ 

[2d 

+ A M\S7u k \{x) wi) x(S pk) dxdydt, 


while 


B , = 


2/ / K h (x - y) (ptw^x) + u k (x) • V^wi + div x it fc (x) x(^Pfc) w>i 

Jo Jn 2d 

+ XM\Vu k \(x) W\ ) wi(y) x(<5 Pk) dxdydt. 


Finally 


D i = -2 


- 2 /) (div u k {x) - div u k (y)) (^x'i^Pk) Pk 
[ 2 d \z 

x($Pk) - -jx!{dPk)dpkj Wi(x) dxdy dt, 


and 


D 2 = - 


n K h (x - y) (div u k (x) - di vu k (y)) x’(8Pk) Pk 

[2d 

w\ (x) w i (y) dx dy dt. 


Note that we have split in several parts a null contribution in terms of the maximal 
function namely the ones with M\X/u k \ 2 ■ Notice also the additional terms in D\ that 
come from cross products such as divrtfc(y) x.X w i( x ) which would pose problems in 
Bi. 

The contributions D\ and are already under the right form. Using Eq. m 
with m , one may directly bound 

Bi<- 2 \ [ [ K h (x - y) (pl(t,x) + P k {x) pk(x) + Rk(x))wi(t,x)x(dPk), 

J o Jn 2d 

and 

B 2 <- 2 \ [ [ K h (x-y)(pl(t,x) + Pk(x)pk{x)+Rk{x))w 1 (x)wi(y)x(dpk), 
Jo Jn 2d 

giving the desired result by symmetry of the expression in x and y. 

The term A 2 is straightforward to handle as well. Use m to get 


Ao < 


[ [ \VK h (x-y)\\x-y\(M\Vuk\{x) + M\Vuk\(y))x(f> Pk)Wz 
Jo Jn 2d 

A [[ K h (x-y){M\Vuk\(x) + M\Vuk\{y))wi(x)wi(y))x(5pk)- 
Jo Jn 2d 


60 


Since |x| \S7Kh\ < CKh, by taking A large enough, one obtains 


A 2 < 0 . 


The term A± is more complex because it has no symmetry. By Lemma 16. II 


Ai < C 


n \X/K h {x -y)\\x-y\ (. D\ x _ y \ u k (x) + D\ x _ y \ u k {y)) \(8 p k ) 

[2d 

wi(t, x) dx dy dt 

X / Kh{x — y) AI |Vufc|(x) w\(t, x) x{& Pk) + similar terms in w\(t, y). 
Jo Jn 2d 


The key problem here is the D k u(y) wi(x) term which one has to control by the term 
M | Vit| (x) w± (x) . This is where integration over h and the use of Lemma RPl is needed 
(the other term in w\(y) is dealt with in a symmetric manner). For that we will add 
and subtract an appropriate quantity to see the quantity D\ x _ y \u k (x) — D\ x _ y \u k (y). 


By the definition of K h , \z\ \X/K k (z)\ < C K h (z) and by (18.11) . x(^Pk) < C (p k (x)+ 
p k {y)) with p k € L 2 uniformly moreover uq in uniformly bounded. Hence using 
Cauchy-Schwartz and denoting z = x — y 

i P5ITT Sit L •*« - + *>■*• T 

+ 4 C / K.h 0 {x-y)D\ x _ y \u k {x)x{8 Pk)Wi 

Jo Jn 2d 

-2X [ [ 1C ho {x - y) M \S7u k \(x) w\{x) x{8pk) 

Jo Jn 2d 

pi P 7 7 

<2 C / K h (z) ||D| Z | u k {.) - D\ z \u k {- + z )\\ L 2 —, 

Jh 0 Jn d n 

by taking A large enough since Lemma [71721 bounds D\ x _ y \ u k (x) by M |Vufc|(:r). 


Finally using Lemma 16.31 


f 1 A 1 dh 
ho ll^lli 1 h 


< C | log /t-o | 1/2 \\u k \\ H i. 


Summing up Ai + Bi + Di 1 and integrating against 
the proof. 


dh 

\\Kh\\ L ih 


for i = 1, concludes 


□ 
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Proof of Lemma 18.21 In this part, we again assume a k = 0 and still use [2Bj to 
obtain the renormalized relation of Lemma EE] with ip = Xa and = K h -kW a = 
w a ,h + Wa,h- With this exception the proof follows the lines of point i) in the previous 
Lemma ISTT1 so we only sketch it here. 

From Lemma 1 7. II we get 


(w a ,h{x) + Wa,h{y)) Xa($Pk)(t) 

( Wa,h{x) + Wa,h{y)) Xa(Spk)\t=0 + A + B + D + I + II + IL a , 


Kh(x - y) 


Jh 0 Jn 2d h 

< f 1 f Kh{x-y) 

Jh 0 Jn 2d h 

with the terms 


A = / (uk(t,x)-u k (t,y))-X7K h (x-y)Xa{5pk)(wa,h(x)+Wa,h(y)), 

Jo Jn 2d 


while 


B = 2 


L J o Jn 3d 


K h (x - y)xa(Spk)(u k {x) - Uk{z)) ■ V x K h (x - z)w a (t , z) 


- A 


o Jn 2d 


Kh{x - y)xa{5pk)I<h * (M|Vu fc |w a ) 


and 


D = - A 


[ Xa(Spk) I<h * ((|div u k \ + \A^p 1 \)w a )(a;) K h (x - y 
o Jh 0 Jn 2d n 

t n r _ _ d/r 

w a ,h(x ) Xa(Spk) Kh * (|div Uk\ + \A IJ/ p' y \)(x) Kh(x - y) 

o Jh 0 J n 2d 


h . 


The dissipation term is under the right form 

n i /* dh 

/ ix a ,h{x) Xa{Spk) Kh * (|div Uk\ + \A IJ ,p' y \)(x) K h (x - y )—, 
o Jn 2d h 


and so are by symmetry 


l r t r 17 


I = - 


K h {x - y) 


h 0 Jo Jn 2d 


(div u k (x) - di vu k {y))x'a{fiPk) Pk w a ,h{x ), 


and 


II = — 


h 0 Jo Jn 2d 


—y) ^[ VUk ^ + di vu k (y)) 


(Xa( S Pk)Spk ~ 2 Xa(Spk))w aih (x). 
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The terms A and B are treated exactly as in case i) of Lemma 18. II they only require 
the higher integrability p > 7 + 1 + t. 

The only additional term is hence D which is required in order to write the 
dissipation term II a in the right form. D is bounded directly by point iv) in Lemma 
17.21 Thus 

A + B + D < C | log h 0 \ e , 

for some 0 < 9 < 1 which concludes the proof. □ 


8.2 The control on the effective flux 


Before coupling the previous estimate with the equation on divrtfc, we start with a 
lemma which will be used in every situation as it controls the regularity properties 
of 

Dpu k = A -1 div (d t (pk Uk) + div (p k u k <S> u k )), 

per 

Lemma 8.3 Assume that p k solves (ED, that ED-ED hold and ED with 7 > 
d/2. Assume moreover that <f> £ _L°°([0, T\ x II 2d ) and that 


:= 


ln d 


Kh(x - y) §{t,x,y) dy 


wy Wx 1,1 


+ 


/ n d 


Kh(x - y) <&{t,x,y) dx 


1.1 TI7-1.1 


< OO, 


WB L W, 


then there exists 6 > 0 s.t. 


n l<h{x - y) $(t,x,y) ( Dpu k (t,y ) - Dpu k (t, x)) dx dy dt 

[d 

<C\\K h \\ L1 h e (||$||l- +<7*). 


Proof This proof is divided in four steps: The first one concerns a control on 
Pk\u k — u k , v | where u k , v is a regularization of u k defined later-on; The second step 
concerns the proof of an estimate for 3> x = ( I\h * $)(£, a;) and = (K k * $)(£, y) 
in LjLP with p' = p/{p — 1); The third step concerns a control with respect to h 
when $ 3 ,, <f> y in ffi W t 1 ’ +oo W _1,0O— 0 with p' = p/{p — 1); The last term is the 

end of the proof obtained by interpolation. 

i) A control on p k \u k — u ktV \ where u k<r) is a regularization in space and time defined 
later-on. Choose a kernel /C^ € C'^°(R + x R + ) s.t. /C v (t,s) = 0 if \t — s| > 77 , for 
smoothing in time. We still denote, with a slight abuse of notation 

^C »7 * t u k {t)= / JC v (t,s)uk(s)ds. 

•j M _|_ 
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Now since u k is uniformly bounded in L/H/. C L/L® in dimension d = 3 and in 
LjrL| for any q < oo in dimension d = 2, one has 


{u k (t,x)-u k (s,x)) 2 

Pk(t,X) - - 77 -ry—j- 7 - 77 AC,(i, s) dt ds dx 

l+\Uk{t,x)\ + \u k (S,x)\ 

— / Pk id i ^0 iV'k idt *^) ^k (s 5 a?) ) AC, (t, s) I\ jj' ~k x 


Uk(t , •) - Ufc(s, ■) 


1 + |ztfc(f, .)| + |wfc(s, -)l 


+ C ( ? /) e l|Pfc||.L“L2i 

with 9 > 0 if 7 > d/2 for d = 2, 3. Note that 

11 pk 11 Wx ~ 1 ’ 1 — 

and by interpolation, as j > d/2 and thus 7 > 2d/{d + 2), there exists 9 > 0 s.t. 


\\d 6 t Pk\ 


L^H- 1 ^ C - 


Thus 


[ \ (Uk(t,x) - u k (s,x)) 2 

/ Pk(t,x)—- — , — 7 -rr AC,(i, s) dtdsdx 

J 1 + \u k (t,x)\ + \u k {s,x)\ 

< / ( Pk(t ) Wfe(t, x) - p k (s ) Mfc(s, x)) s) K v > * 


.) - n fc (s,.) 


" " 1 + Mt,.)| + Ms,.)l 


+ C (77') 61 + C \\u k \\ L 2 H i . 

Using (15.611 . one deduces that 

f u ^ (u k (t,x) - u k {s,x)) 2 

/ Pk[t,x)—— :— — 7 -y )C v {t,s)dtdsdx 

J l + \u k (t,x)\ + \u k {s,x)\ 

<cv ,e + c^- d +c^- d . 

Optimizing in 7 / and interpolating (using again 7 > d/2), one finally gets that for 
some 9 > 0 

J p k {t, x) ( 7 tfc(t , x) — u k {s, x)) 2 s) dtds dx < Cb/. 

Define for some v > 0 

= K'lj 1 ' *x AC^ *£ U k . 

Using the regularity in a: of u k one has that 

J Pk(t,x) I u k (t,x) - U k , v (t,x )| 2 dxdt < Crf, 
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and finally by m , there exists 9 > 0 s.t. 


||Pk [uk — Uk, v ) || L i+o L i+o < Crf. 


(8.3) 


ii) The case where $ 2 :, is only in with p' = p/(p — 1). We recall that p is 

the exponent in (15.611 . Denote 

= f [ K h (x - y)<&(t,x,y) (Dpu k (t,y) - Dpu k (t,x))dxdydt 
Jo Jn 2d 

which can be seen as a linear form on $. Recall as well 

$x= K h (x - y)$(t,x,y)dy, $ y = K h (x - y) $(i, x, y) dx. 

Ju d Jn d 

By (15.61) . Dpuk is uniformly bounded in Therefore 


|J*|<C'||J0,|| L 1 (||$x|| i?ir +\\*v\\J 


')• MM>°- < 8 - 4 ) 


P P 

iii) The case <& x , in L^L^. n W t 1 ’ +oo W _1,oo_0 with p ' = p/(p — 1). Denote 

and 


— || ‘hxll T 2 T p' "b II ^s/ll r2 T p' II^MI w. 1 ’ 00 w-Uoo-o + ll^ > yllry t 1 ’ oo W /_1 ’ oo_0 > 


Ri=A x di vpk(u k -u k , v )- 
Observe that by (18.31) and integration by part in time 

r t r 

$ x d t Ridxdt < C^y 6 \\K h \\ L i. 


0 Ju d 

The same procedure can be performed with div ( p k u k <8>u k )- Denoting 
Dpuk, v = A -1 div (d t (p k u ktV ) + div (p k u k ® u ktV )), 

one then has 

I*<C*\\K h \\» rf 

+ K h (x-y)$(t, x, y)(Dpu k ,r,{t,y) - Dpu k ,r,(t, x)) dxdydt. 

Jo Jn 2d 

However using m 

d t (p k u ktV ) + div (p k u k ® u k , v ) = Pk{d t u k:V + u k ■ Vu*,,,) + a k u k , v A p k . 
For some exponent n 

||A _1 div (p k {d t u ktTI + u k -Vu ktV ))\\ L i w i,i <Crj~ K , 
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and 


Therefore 


O-k || A x div (a k u k , v ^Pk)\\ L 2 L 2 <Cr] K y/ok- 



y)®{t,x,y) (Dpuk,r,(t,y) 


<C®r] h \\KhWL 1 (h + y/ak) 1 °. 


Dpu k , v {t, x)) dx dy dt 


Finally 


I$<C*\\K h \\ L i ( V e + rj- K (h + ^)), 
and by optimizing in ip there exists 9 > 0 s.t. 

i*<c\\K h \\ L ih e (\\$ x \\w?-w-i.~-° + \\*v\\ w ;-°°w-i"°-°) ■ ( 8 - 5 ) 

iv) Interpolation between the two inequalities (EH) and (EUD- For any s £ (0,1) 
there exists 9 > 0 s.t. 

I § < C H-fOillii h + ll^ > ylli2 L p' 

+ ll^xH^.-J + O^-^r+O + H^Jw-l.S+O^-l.r+o), 

with 

1 _ 1-s 1 _ 1-s 

q 2 ’ r p' 

On the other hand, if for example <P X belongs to Lf x and to W t 1 ’ 1 W'“ 1 ’ 1 then by 

interpolation is in W - /’ 1 /* -0 Hence C<f, controls the W t s ’ 9+0 W'“ s,r+0 

norm provided 

, , 1-s 1-s 

s < 1 q = —-—, s < 1/r = —. 

2 p' 

This is always possible by taking s small enough (but strictly positive) as p > 1 and 
hence p' < oo. This concludes the proof. 

Note finally that the interpolations between Sobolev spaces are not exact which 
is the reason for the 1/s — 0 or q + 0 and r + 0 (one would have to use Besov spaces 
instead, see for instance [49]). □ 


8.3 The coupling with the pressure law 

We handle all weights at the same time. For convenience, we denote 

Xi(t,x,y) = ^x'ifiPk) Pk + x(Pk) ~ \x!{5pk)Spk- 
In the case without viscosity, one has 
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Lemma 8.4 Assume that p k solves (ED with a k = 0 , that m , (i53i) . m with 
7 > d/2 and p > 2 hold. Assume moreover that u solves (15.21) with p k compact in 
L 1 and satisfying (& F k compact in L 1 , Pk satisfying m- 
i) Then there exists a continuous function e(.) with e(0) = 0 s.t. 


n K h (x - y) (divuk(x) - divu k (y )) Xi wi(x) < C \\K h \\ L i e(h) 

[2d 

+ C / K h (x - y) (l + p/{x) + P k {x) p k {x) + R k (x)\ x(fipk)wi(x). 

Jo J v ' 

ii) There exist 6 > 0 and a continuous function e with e(0) = 0, depending only on 
p and the smoothness of p k and F k , s.t. 

n K h (x - y) ( divu k (x) - divu k (y )) x'($Pk) Pk wi(x) wi(y) 

[2d 

<C\\K h \\ L i (e(h) + h e )+C [ [ K h (x - y) (l + pl(x) + P k (x) p k (x) 

Jo J n 4d v 

+ R k {x) + p/(y) + P k (y) p k (y) + Rk{y)) x(^Pk)wi{x) w k {y). 


For instance if p k and F k belong to W S}1 for some s > 0 then one may take e(h) = h e 
for some 9 > 0 . 

While in the case with viscosity 

Lemma 8.5 Assume that p k solves m, that m , 33i) . m with 7 > d/2 and 
p > 2 hold. Assume moreover that u solves m with p k compact in L 1 and satisfying 
(15.31) . F k compact in L 1 , P k satisfying (EHD- Then there exists a continuous function 
e(.) with e(0) = 0 and depending only on the smoothness of p k and F k s.t. 


1 

2 

1 

2 

A 

2 



y) ( divu k (x ) - divu k (y))x'{bpk) Pk W 0 K h {x - z) K h (y - 

y)(divu k (x) + divu k (y))(x'{Sp k )Spk - 2x($Pk)) 
W 0 (z,w)K h {x - z)K h {y - w) 
y)x(dpk) Kh(x - z) M|Vufc| w 0 (z ) 

nt r 


<C\\K h \\ L i£(h) + C 


/ / K h (a 

Jo J n 4d 


y) x($Pk) Wo K h (x 


z ) K h (y — w). 


w) 


Proof The computations are very similar for i) and ii) in Lemma 18.41 and for 
Lemma [83] For simplicity, in order to treat the proofs together as much as possible, 
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we denote 

Gi(t, x, y) = xi{t,x,y)w l (x), G 2 (t,x,y ) = x!(&Pk) PkWi(t,x)w 1 (t,y), 
Go(t,x,y) = l;Xo( s Pk)Pk f W 0 (t,z,w) K h (x - z) K h (y - w)dzdw. 

2 J n 2d 


The first step is to truncate: Denote I k {t,x,y) = <f>(pk{t,x)/L)(f)(pk(t,y)/L) for 
some smooth and compactly supported <p, 


n K h (x - y) (div u k (x) - di vu k (y)) Gi 

[2d 

<C\\K h \\ L i L~ e ° - [ [ K h (x - y) (divM fc (:r) - divu k (y)) G 
Jo Jn 2d 


T l 
i 1 k ■ 


Here for i = 0,1, 2, Gi < C (p k {t,x) + p k (t,y)) (even G 2 < 2) and consequently, as 
divrtfc £ L 2 uniformly, only p > 2 is required with 9q = (p — 2)/2 > 0. 

Introduce an approximation pk :V of p k , satisfying m and s.t. 


Il^fc.^llw, 2 ^ 00 — GV i ^ £o(v)i 

/ / K h (x-y)\p k ,r,{t,x) - p k ,v(t,y)\dxdydt< \\K h \\ L ie 0 {h), 

Jo Jn 2d 

from (15.141) . Use (15.21) to decompose 


( 8 . 6 ) 


J I< h (x - y ) (di yu k (x) - di vu k {y)) Gi I k dxdy 

= 2 Ai + 2 Bi + 2 Ei, 


with 


and 


r t l 

Ai = - K h ( x - y) (. Pk{pk(x )) - Pk(pk(y))) Gi -r dxdy, 

J Pk,rt{X) 


Bi= K h (x - y) F k (t , x, y) Gi 


J L 

1 k 


Pk,r](x) 


dx dy, 


where 

Fk(t,x,y) =F k (t,x) - F k (t,y) + p k (y) Pk, v (x) div u k (t, y) 
- p k (x) Pk, v (®) div u k (t, x) 


Pk, v (x ) pk{y) 


l 


Pk,ri(x) Pk{x) J 


Finally 


f I L 

Ei= K h {x - y) (Dpukit, y) - Dpu k (t, x)) Gi - dx dy, 

J Pk,r,(x) 
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with as before 


Dpuk = A 1 div (d t {pk Uk) + div (p k Uk ® Uk))- 
For Bi by the compactness of F k , Pk, estimates (18.611 and (15.141) . and by (15.31) 


Bi < CL 


o J n 2d 


K h {x-y) \F k (t, x, y)\ < C L (so(h) + £o(v)) ll-MU 1 - ( 8 - 7 ) 


Note that again \Gi\ < C ( p k {x) + p k (y )) for i = 0, 1, 2. 
For Ei , we use Lemma 18.31 bv defining simply 

1 


®i{t,x,y) = Gilfc(t,x,y ) 


Pk,Tl(t,X) ' 


By (ED, 


1 $ 


i L° 


< CL. 


As for the time derivative of <f>, for i = 1, 2, Gi is a combination of functions of 
Pk(t, x ), pk(t, y) and Wi which all satisfy the same transport equation (with different 
right-hand sides). By (15.11) . 


d t Gi + div x ( u k (x ) Gi) + div y (u k (y) Gi) = fi,idiv x Uk(x) + f 2 ,i div y u k {y) 

+ h,iDi(x) + U,iDi(y), 

where the Di are the penalizations introduced in Section 17.21 and the are again 
combinations of functions of p k {t, x), p k (t, y) and Wi. Every /„^ contains as a factor 
4>{pk (x) /L) or a derivative of <j> and thus 


\\fn,i\\L°° < CL, Vn, i. 


Finally by the smoothness of 

d t ®i + div x ( u k {x ) $i) + div y (u k (y) $i) = fipdiv X u k (x) + h,i div y u k (y) 

4” fz,iDi(x) -f- f<i,i Di(y) ~i~ Qrii 

and it is easy to check that the constant C'$ i as defined in Lemma 18.31 is bounded 
by C L ? 7 -1 . 

The case i = 0 is slightly more complicated as Wq is integrated against Kh so the 
equation on $o involves non local terms and we have to take into account extra terms 
as mentioned in the statement of Lemma 8.4. By (17.21) . denoting ino,/i = K] x -kw o 

dtw 0ih +u k {t,x) ■ S7w 0 , h - a k A x w 0 ,h = ~K h * ( D 0 w 0 ) + R h - K h * (divu fc w) 0 ), 


with 


Rh 



VK h {x 


z) ■ ( Uk(t,x) 


Uk(t, z)) Wo(t, z) dz. 
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Remark that Rh is uniformly bounded in L 2 x by usual commutator estimates. 
Finally as p k , v is smooth in time, one has 

<9 t $ 0 + div x ( u k (x ) $ 0 ) + div y ( u k (y ) $ 0 ) - ot k (A x + A y ) $ 0 
= f lj0 div X u k {t,x) + f2,0 div v u k (t,y) + a k (/3,o \^xPk(x )\ 2 + / 4 ,o \^xPk(y)\ 2 ) 

- — '■— (K h * ( D 0 w 0 ) + Rh - K h * (div ii fc w 0 )) - 2 • V x w 0 ,h 

+ ^0 9ri, 

where <b p = 5p k p k I k (t,x, y)> 9v a function involving first and second derivatives 
of pk,y in t and x and Viifc. The fj ,q are combinations of functions of p k (t,x) and 
Pk{t,y ), multiplied by w 0 ,h, and involving </>(p k (x)/L), (j)' (p k (x) / L), or (/>"(p k (x)/L) 
and the corresponding term with p k (y). 

By the L°° bounds on $ p , wo, each fj ,o and by m , one obtains 


d t / K h (x - y)$ 0 (t,x,y)dy 


l\ w x ~ 


< CLy- 1 . 


Therefore C<t> < C Ly Thus for all three cases, Lemma 18.31 yields 


EiKCLrj-'WKhWLih 0 , 


( 8 . 8 ) 


for some 8 > 0. 


Proof of Lemma Iff. 51 The term ^4o • The terms Aj are where lies the main 
difference between Lemmas 18.41 and 18.51 as P k is not monotone in the first case and 
monotone after a certain threshold in the second. For this reason we now proceed 
separately for Lemma 15751 and Lemma l8. ll In the case with diffusion for Lemma [5751 
there also exist extra terms to handle, namely J +1 with 

J = “ [ [ K h (x — y)x(5pk)Kh(x — z) M\Vu k \(z) w 0 (z), 

1 Jo Jn 3d 

and 

I = ~\f [ K h (x - y) (div u k (x) + div u k (y)) 
z Jo J n 4d 

(x'(Spk)Sp k -2x(Spk)) W 0 I<h(x - z) K h (y - w). 


We decompose this last term in a manner similar to what we have just done, first of 
all by introducing the truncation of p k 


I<C\\K h \\ L xL- e ° 



y) (div u k (x) + div u k (y)) 

(x($Pk)Spk -2x(Spk))I k Wq K h (x 


z) K h(y ~ w), 
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with again 9q = (p — 2)/2. Now introduce the p k 

I < CHA^Hii L e - i [ f K h (x-y)(p k (t,x)divu k (x)+p k (t,y)divu k (y)) 
z Jo J n 4d 


1 k 


p k<v (t,x) 


{x($Pk)Spk -2x(Spk)) W 0 K h (x - z)K h (y - w) 


I K h (x-y)H k (x,y)(x'Spk-2x(Spk))IkW 0 K h (x-z)K h (y-w) 
z Jo Jn id 


where 


H k (t,x,y) = p k (x)dWu k (x)(- 


1 


K p k , v {x) p k (x) 
By the compactness of p k , one has that 


) - Pk(y)divuk(y)( 


Pk{y) Pk, v {x) 


)■ 


K h {x - y)H k (x, y)(x! $Pk ~ 2 x{&Pk)) W 0 K h (x - z) K h (y - w) 

/o dn 4d 

< \\K h \\ L ieo(h)\\u k \\ L?H i\\p k \\ L 2 <Ce 0 {h)\\K h \\ L i. 


This implies that 


I<- 


77 [ [ K h(x - y) (p k (t,x)divu k (x) + p k (t, y)div u k (y)) 
z Jo J n 4d 


pk,r) (ti * 


(x'(dpk) Spk - 2 xi s Pk)) Wo K h (x - z) K h {y - w) 


+ C\\K h \\ L i(L-» +e 0 (h)). 

Using (15.21) or namely that p k d\vu k = Dpu k + F k + P k (p k ), the quantity Aq + I + J 
may be written 


Aq + I + J < C \\K h \\ L i(L + £o(h)) + I\ + I 2 


with 


h = A 0 — 


[ f K h (x-y){P k (x,p k (x))+P k (y,p k (y))) -- 

Jo J n 4d Pk,ri\t,' 


(8.9) 


0 


ix'idpk) Spk -2x(Spk)) Wo Kh{x - z) K h (y - w), 


and 


h = ~ 


7 ; [ ( K h(x - y) (Dpu k (x) + F k (x) + Dpu k {y) + F k (y)) -—- 

^ JO J Il 4d f J 'k,rjy'^ x ) 

(x'{fipk) Spk -2x{Spk)) W 0 I\ h (x - z) K h (y - w ) 
K h (x - y)x(Sp k )K h (x - z) M\Vu k \(z) w 0 (z). 


0 Jn 3d 
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In this case with diffusion, because P k is essentially monotone, the term Aq is mostly 
dissipative and helps control the rest. More precisely 


h = --z / / K h {x-y) [{P k {x,p k (x)) - Pk(y,pk(y)))x'(8p k )p k 
& Jo Jn 2d L 

+(Pk{x, pk{x)) + P fc (y, pk{y))){x'($Pk)Spk - 2 x($Pk))\ 

I L r _ _ -I 

- J- —r- / W 0 K h (x - z) K h (y - w)dzdw dxdy. 

Mfc ,%) JYl 2d 

As P k > 0 and by (EH), x'i^PkjSpk - 2 x($Pk) > -x'{$Pk)Spk, thus 

C P k (x,p k (x )) - Pk(y, Pk{y)))x!Pk + (Pk(x,p k {x)) + Pk{y,pk{y)))(x'6pk - 2x(6p k )) 
> x'{$Pk) [{Pk{x,pk{x)) - Pk{y,Pk{y))) Pk - (Pk(x,p k (x)) + P k {y,Pk(y))) Sp k \ ■ 


Without loss of generality, we may assume that pk(x) > p k {y) and hence x'(^Pk) > 0. 
Developing 

(P k (x,p k (x)) - Pk(y,pk(y))) Pk - (■ Pk(x,Pk(x )) + Pk{y,pk(y)))5p k 
= 2 P k {x, Pk(x)) Pk(y) -2 p k {y,Pk(y)) Pk(x). 


We now use the quasi-monotonicity (15.811 of Pk(x,s)/s. First of all if p 0 < p k (y) < 
Pk{x) then necessarily Pk depends only on pk(x) or pk(y) plus Pk . Thus 

P k {x, Pk(x)) p k (y) - Pk{y, Pk(y)) Pk(x) > -\P k (t,x) - P k (t,y)\. (8.10) 


If pk(y) < po, by 05.81) Pk(x,s) —> oo as s —> oo while P k (y, p k (y)) is bounded. 
Hence there exists p large enough with respect to po, s.t. if p k {x) > p, then again 


P k {x, p k {x)) pk(y) - Pk(y, p k (y)) Pk(x) > 0. 


The only case where one does not have the right sign is hence where both p k {x) 
and Pk(y) are bounded by p and po- Therefore using the local regularity of P k given 
by (EED' 


(P k (x, p k (x)) - P k (y, Pk{y))) Pk - (. Pk(x,p k {x )) + P k (y,Pk(y)))Sp k 

> -P\6p k \ - Qk(t,x,y). (8.11) 

Introducing these estimates (18.101) and (18.111) in Ii yields 


h<P 


K h (x - y ) (| P k (x) - P k {y)\ +Q k + \S Pk \) 
o Jn 2d %) 


I u 2d 


Wo K h (x— z) K h {y- w) 


( 8 . 12 ) 


< £q (h) || A/,11^1 + P 


o J n 4d 


Kn(x - y)x(Spk) W 0 K h (x - z) K h (y - w). 
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Turning now to I 2 , we observe that p k M\Vu k \ > p k di vu k > Dpu k + F k and that 
X(4) > (2x(4) - x'(Sk)6 k )/C. Therefore for A large enough, using W 0 (z,w) = 
wq (z) + w>o (w) and the symmetry 


I 2 < — 


7 ; [ [ K h (x- y) {Dpu k {x) - Dpu k (z) + F k {x) - F k (z) 

z Jo Jn 4d 


+Dpu k (y ) - Dpu k (z) + F k (y) - F k (z)) 


;(x'(8p k )6pk -2 X (6p k )) 


W 0 K h (x - z) K h (y - w). 


The differences in the F k are controlled by the compactness of F k and the differences 
in the Dpu k by Lemma 18.31 as for the terms Ei. Hence, finally 

h <C\\K h \\ L i (so{h) + Lp- 1 h e ). (8.13) 


Conclusion for Lemma Iff.51 We sum up the contributions from Bq in (18.71) , Eq 
in (18.81) , A 0 + I + J in (18.91) together with the bounds on I\ in (18.121) and I2 in (18.131) 
to obtain 


7 :f I K h (x — y) (div u k (x) — div u k (y)) G 0 
z Jo J n 2d 

I K h(x-y) (divu k (x)+divu k (y))(x'(Sp k )Sp k -2x(Sp k )) 
z Jo Jn id 

W 0 (z, w) I\ h (x - z) K h (y - w) 

K h (x - y) x(dpk) K h (x - z) M |Vu fc |( 2 ;) w 0 (z) 


0 Jn 3d 

< C \\K h \\ L i ( L~ e ° + L (eo(fc) + e 0 (v)) + L r,- 1 h e ) 

+ C [ [ K h (x - y) X (Spk)Wo(z,w) K h (x - z) K h (y - w). 

Jo Jn id 

Just optimizing in L and p leads to the desired e(h) and concludes the proof of 
Lemma 18.51 


Proof of Lemma \ 8.f\ The term Aj with i = 1, 2. It now remains to analyze 
more precisely the terms (P k (x, p k (x)) — P k (y, p k (y))) Gi for i=l,2 concerning the 
case without diffusion but with non monotone pressure. We will split the study in 
three cases but remark that now the possible dependence of P k in terms of x affects 
the estimates. For this reason, we carefully write explicitly this dependence. 

Case 1). The case ( P k (x , p k (x)) — P k (y , p k (y)))5p k > 0. Since G 2 obviously have the 
same sign as Sp k , one simply has 

( Pk(x,p k (x )) - Pk(y,Pk(y))) Gi > 0 , 
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for i = 2. In the same case, for G\ 


( Pk{x,p k (x )) - Pk(y,Pk(y))) Gi 
= (Pk(x,p k (x)) - P k {y,Pk{y))) 

> \Pk{.x, Pk(x)) - Pkiy, p k (y))\ 


\x!{Spk) Pk + x(Spk) - ]pi{5pk) Spk 


\x'{tipk)\p k - 


x(Spk) - ^x\Sp k )Sp k 


m(x) 


wi(x) 


> 0 , 


by m as pk > \5pk\- Therefore in that case the terms have the right sign and can 
be dropped. 


Case 2). The case (P k (x,p k (x)) - P k (y, p k (y)))5p k < 0 and p k (y) < Pk{x)/2 or 
Pk{y) > 2 pk{x) for some constant C. 

- For i = 1, first assume that P k {x,p k (x)) > P k {y,p k (y)) while p k (y) > 2 p k {x). 

( P k {x,Pk{x )) - P k (y,Pk(y)))G 1 > -P k (p k (x)) {\x!{&Pk)\ Pk + x(fipk)) VJi(x). 
Now observe that since p k (y) > 2 p k ( x) then 

\x\fipk)\Pk < ^ \x\6pk)\p k {y) < 3 \x'{dpk)\ \Spk\ < Cx(Spk), 

by (E3J. Therefore in that case by m 

(P k (x, p k (x)) - P k (y, Pk{y))) G 1 > -C (p k (xy + R k (x)) x($Pk) «h(x), 

Note that the result is not symmetric in x and y: We also have to check also 
P k {x,p k {x)) < P k (y,Pk(y)) and p k {y) < p k (x)/2. Then simply bound since now 
Pk{y) < Pk(x) 

( P k (x,p k (x )) - P k {y,p k {y)))Gi > -C ((pkiy)) 7 + R k (y)) x{fipk)wi(x) 

> -C((p k (x)y + Rk(y)) x{$Pk) wi(x). 


In both cases, one finally obtains 

(P k {x,Pk(x)) - Pk{y,Pk(y)))Gi >- C {(p k (x)y + Rk(x)) x(^Pk)wi{x) 

- |i?(x) - R(y)\pk wi(x). 


- For i = 2. The calculations are similar (simpler in fact) for G 2 and this lets us 
deduce that if P k (x, Pk{x)) — P k (y, p k {y)) and p k {x) — p k {y) have different signs but 
Pk(y) < Pk(x)/2 or p k (y) > 2 p k (x) then 

(P k [x,p k {x)) - Pk(y,Pk{y)))G 2 > -C((p k (x)y + P k (x)p k {x) 

+ (Pk(y)y + P k (y)Pk(y)) x(^Pk)w 1 (x)w 1 (y). 
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Case 3). For i = 1, 2, the situation where P k (x, Pk(x)) — Pk(y, Pk(y)) and pk(x) — 
Pk{y) have different signs but Pk(x)/ 2 < pk(y) < 2 pk(x). Then one bluntly estimates 
using the Lipschitz bound on Pk given by m 

I (Pk(x,pk(x)) - Pk{y,Pk{y)))\ < C((p k (x)y ~ 1 + P k (x) 

+ (pk(y)y~ 1 + Pk{y)) \Spk\ + Qk- 

Bounding now the Gi by m, 

( P k {x,p k (x )) - Pk{y,Pk{y))) G 2 <C (( Pk(x)y + Pk(x)p k (x) + (Pk(y)y 

+ Pk(y) Pk(y )) x($Pk ) wi(x) wi(y) + Q k Pk wi(x), 


and 

(Pk{x,p k (x)) - Pk{y,Pk(y)))Gi < C((p k (x)y + P k (x)p k (x) + (Pk{y)y 

+ P k {y) Pk(y)) x($Pk) uii{x) + Q k p k 
<C( 1 + p1(x) + Pk(x) p k (x))x(Spk)wi(x) 

+ (Qk + | Pk(x) - Pk(y)\p k ) p k wi{x), 

as pk(x) and Pk(y) are of the same order. From Prop. 17.21 point i. we know that 
Wi(x) < e~ Xpk ^ v . One the other hand, we are precisely in the case where pk(x) 
and Pk(y) are of the same order. Hence p l k w i is uniformly bounded for any l > 0. 
Hence we finally obtain in this case 

(Pk(x,Pk(x)) - Pk(y, Pk(y))) G 2 < C(( Pk (x)y +Pk(x)p k (x) + (pk(y)y 

+ Pk(y) Pk(y )) x($Pk) Wi(x) w\(y) + Q k , 


and 

(Pk(x,pk(x)) - P k (y,Pk(y)))G i < C(l + pl(x) + P k (x) pk(x)) x(Spk) wi(x) 

+ Qk + | Pk(x) — Pk(y) I- 


From the analysis of these three cases, one has that 

M <C J K h (x - y) (1 + pl(x) + P k (x)pk(x) + R k (x))x(Spk)w 1 (x)dxdydt 

+ C J K h (x - y) (Qk + \P k (x) - Pk(y) | + | Rk(x) - Rk(y) |) dxdydt. 

Therefore by the compactness properties of P k and the estimate on Q k in the as¬ 
sumption (15.81) 

M <C J I\ h (x - y) (1 + p1(x) + P k (x)pk(x) + Rk(x)) x(Spk) wi(x) dxdy dt 
+ II^IIl 1 £o(h), 

(8.14) 
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and 


A 2 <C J K h (x - y) (1 + pl(x) + P k (x) p k (x) + R k (x) + pl(y) 

+ Pk(y ) Pk(y) + Rk(y))x(Spk)wi(x) w k (y) dxdy dt + ||^||i,i e 0 (h). 


(8.15) 


Conclusion of the proof of Lemma \ 8.f\ Summing up every term, namely (E3>-(E 
and (I8.14I1 - (I8.15I) . we eventually find that 


n K h {x - y) (div u k (x) - divw fc (j/)) Gi 

i 2 d 

< C \\K h \\ Ll ( L~ e + L (e 0 (h) + e 0 (r,)) pLrf 1 h e ) 


C / K h (x - y) (1 + pkix) 1 + P k {x) p k {x) + R k )x(Spk)wi(x), 
Jo J n 2d 


while 


n K h {x - y) (div u k (x) - divu k (y)) G 2 

[2d 

< C \\K h \\ L i ( L~ e + L (e 0 (h) + £ 0 {y)) + Lr,' 1 h e ) 

+ C [ [ K h (x - y) (l + p k (xy + P k (x) p k (x) + R k (x) + p k (xy 
Jo J n 2d v 

+ Pk(y ) pk(y) + R k {y)^ x{$Pk)w i(x) wi(y). 

To conclude the proof of Lemma 18.41 one optimizes in ?? and L. Just remark that 
since the inequalities depend polynomially in L and rj then the result depends on e s Q 
for some 6. □ 


8.4 Conclusion of the proofs of Theorems 15.11 and 15.21 

We combine Lemma IBTT1 with Lemma [5~5l or [H~TI and we finally use Prop. 17.21 Let us 
summarize the required assumptions. In all cases one assumes that p k solves (EB 
and that divi^ is coupled with pk through (|5.2[) ; bounds are assumed on the viscosity 
as per m , on the time derivative of p k u k per m , on u k per (15.51) . Finally the 
viscosity p k and the force term F k are assumed to be compact in L 1 . Moreover 

• In the case with diffusion, a k > 0, one assumes that the pressure term P k 
satisfies (15.81) and the bounds (15.41) on p k with 7 > d/2 and p > 2. 

• In the case without diffusion, a k = 0, one needs only (15.91) on the pressure P k 
and the bounds (15.41) on p k with 7 > d/2 and p > 2. Moreover for Prop. 17.21 
it is necessary that p > 7 (in general 7 = 7 < p so this is not a big issue). 
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Then one obtains by taking A large enough and using a simple Gronwall lemma 

ih 
h 


/* 1 n dfi, 

/ / K h (x - z) K h (y - w) (w 0 (t, z) + w 0 (t, w)) K h (x - y) xi^Pk) dx dy dz dw— 

Ih 0 Jn id 


I n id 


K h (x - z) K h (y - w) ( w 0 (t, z ) + w 0 (t, w)) JC ho (x - y) x(<5 p k ) dx dy dz dw 


< C (|logh 0 | 1/2 +£h 0 (k) + e{h)^pj, 


(8.16) 


and 


dh 


x) + wi(t, y)) K h (x - y) x{Sp k )dx dy 
lh 0 Jn 2d n 

/ (wi(t,x) + wi{t,y)) K, ho {x - y) x{Spk) dx dy 
J n 2d 

< C log h 0 1 1/2 + e{h) , 


(8.17) 


with finally 


>n 2d 


wi(t,x)wi(t,y) K h (x - y)x(dpk) dxdy 


<C\\K h \\ L i (h e + e(h)), 

where e depends only on the smoothness of p k and F k and of p > 2. 


(8.18) 


The key point in all three cases is to be able to remove the weights from those 
estimates. For that, one uses point ii) of Prop. 17.21 

The case with iuq(x) + wo(y)- Denote u} v = {K k *Wo(t,x) < y} C [0, T] x II d and 
remark that 


/ n 2d 


£h 0 (x - y)x(dpk) = 


I h 0 J n 2d 


Kh(x - y)xi s Pk) 


fj 

J ho J x€ 


Now 


— r dh 

K h {x - y) x(°Pk) -r 

Ih 0 J igu' or y€uic 

[ [ K h {x - y) x(Spk) Y' 

'ho JxEujn and yGtJri 


dh 


'ho Jx£gjc or 


K h {x - y)x(dpk) ~r 


< 


1/ . 

V Jho J n 2d 


K h (x - y)(K h -kw 0 (x) + K h *w 0 (y )) x(dpk) 


dh 
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while by point iii) in Prop. 17.21 using that p G L p ((0,T) x II d ) with p > 2 and 
recalling that %(£) < C |£| 


Jho JvGuJn and yGuJ-q 


< C | log h 0 1 
| log 7 ?| 1/2 ’ 


Therefore combining this with (18.161) . one obtains 



y)x($Pk)dxdy < C 


' £ hoik) + I log ho 1 1/2 + fk 0 e(h) 


\ dh 
h 


I log ?y | 1/2 J 


recalling that 


and denoting 


E ho (h) — 




e(h 0 ) 


1 

I log h 0 | 



Remark that e(ho) 0 since e(h) —> 0: For instance if e (h) 
| logh 0 | _1 . The estimate then reads 


h e then e(ho) ~ 


f K-h 0 {x-y)x{dPk)dxdy <C (\\ogh 0 \ akh ° + I log I - + £ i h o ) 

J n 2d V P 

llogTyl 1 ^;- 

As HAChJii ~ | log ho|, by optimizing in i], the following estimate is obtained 



y) x{dpk) dx dy < 


_ c\\ic ho \\ L , _ 

I log(ofcho 2 + | log h 0 1 : / 2 + ^(ho ))! 1 / 2 


Per Prop. S3 this gives the compactness of 6pk as the r.h.s. is negligible against 
\\IC ho \\ L i as h 0 —> 0. And it proves the case i) of Theorem 15.II 


The case with w\{x) +W\{y). Similarly, from (18.171) . one then proves that in the 
corresponding case 


I K. ho (x-y)x(Spk)dxdy<C\\ogh 0 \ - +£ ^°^ 

J n 2d V V 

< C\\lC h0 \\ L , 

- | log(| log ho) -1 / 2 + e(h 0 ))| e ’ 

again using part ii) of Prop. 17.21 


1 

I log 771 e 
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In both cases, using Prop. [42] together with Lemma [6.61 in the second case, one 
concludes that pk is locally compact in x and then in t, x. Thus we’ve shown the 
case ii) and concluded the proof of Theorem 15.11 

The case with w\{y)w\(x). The situation is more complicated for (18.181) and the 
product wi(y)wi(x). Indeed wq(x) + wq( y) or w\(x) + w\(y) are small only if both 
wq{x ) and wo(y ) are small (or the corresponding terms for w\). But w±(x) w\(y) can 
be small if either w\{x) or w\(y) is small. This was previously an advantage with 
then simpler computations but not here and (18.181) does not provide compactness. 

This is due to the fact that one does not control the size of {w w < ??} but only 
the mass of pk over that set. The difference between the two is the famous vacuum 
problem for compressible fluid dynamics which is still unsolved. 

The best that can be done by part ii) of Prop. 17.21 is for any rj, 7 / 

/ l Pk(x)>r,\ k (y)> v K h (x - y)x{bPk) < ~fo [ w 1 (x)w 1 {y)K h (x-y)x(Spk) 

J n 2d '/ J 

, r ll^ftlU 1 

rj 1 / 2 | \ogr ]'\ e / 2 ’ 

using that pk G L 2 uniformly. Using (18.181) and optimizing in rf, one finds for some 
6 > 0 


L Ipk{x) - VIpfe(w) -’ 7 Kh{ - X V) X ^ Pk) - ° 77V2 | l 0 g(e(/l) i + h»)\°/ 2 • 

If pk and Fk are uniformly in IV s ’ 1 for s > 0, then 

J^ 2d \ k (y)>v K h{x - y) x(Spk) < C ^i/2 | j 0 g^|s/2 ’ 

which concludes the proof of Th. 15.21 Note however that in many senses (18.181) is 
more precise than the final result. 

8.5 The coupling with the pressure in the anisotropic case 


In that case we need the weight w a and its regularization w a ,h, defined by (17.21) with 
m in order to compensate some terms coming from the anisotropic non-local part 
of the stress tensor. 

Lemma 8.6 There exists U* > 0 s.t. assuming that pk solves (EH) with ctk = 0, 
that m , & eh with 7 > d/2 and p > 7 + 1 + £ = 7 2 /(y — 1 ) hold where 
t = 1 /( 7 — 1); assuming moreover that Pk satisfying (15.81) . that u solves (15.111) with 

a M < C*. (8.19) 
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Then there exists 0 < 9 < 1 s.t. for \ a verifying (18.2J1 for this choice of £ 



[ ^-^7 -—( W a ,h( x ) + W a ,h{y )) Xa(Sp k )(t ) 

Jj2d fl 

< f [ Rh ^ X - —{w a , h { X )+W at h(y))Xa{$Pk)\t=O+C{l+£)\\ogho\ 0 . 

Jh 0 Jn 2d n 


Proof To simplify the estimate, we assume in this proof that P k {pk) — P k , the 
extension when P k satisfies instead [5751 being straightforward. We also recall that \ a 
satisfies (f8~2l) . meaning that for all practical purposes Xa(f) ~ |£| 1+£ . 

We use the point iii) in Lemma HOI 


f f —— [w a ,h{x) + w at h(y))Xa(8pk)(t) 

h 0 Jn 2d n 

-f [ ^"^7- —(Wa,h( X ) +Wa,h{y))Xa{Sp k )\t=0 

Jh 0 J n 2d n 

< C\ log ho\ 6 + I + II — n a 


where 0 < 9 < 1 and with the dissipation term by symmetry 


n a = A 

)o Jh 0 Jn 2d 
while still by symmetry 


Wa,h(x)xa(Sp k ) K h -k (|divufe| + \Afj,Pk(pk)\){x) Kh-j- 


1 = -- 
2 


1 r 17 


Kh{x - y ) 


h 0 do Jn 2d 


(div u k (x) - di vu k {y))x'a(.fipk) Pk w a ,h(x), 


and 


II=~[ I [ ———i ,—— (div u k (x) + div u k (y)) 

1 Jh 0 Jo Jn 2d ", 


{Xa( S Pk)Spk ~ 2Xa(Spk))Wa,h(x). 


I) The quantity I. We recall that in this case one has the formula (15.111) on divzifc 
div Uk =r , kPk{pk) + r'k a^A^Pk(pk) 

1 (8-20) 

+ v k {Afj, - a^Ek) di v(d t {pkU k ) + div (p k u k <8> u k )), 

leading to the notation 

Dpuk = Vk {- a^£(fc) _1 div (d t (p k u k ) + div(p fc u k ® u k )). 
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Therefore, one may decompose 


with 


while 


and 


In = — 


I = J 0 + I D + I R , 
K h (x - y) 


(Dpu k (x) — Dpu k {y)) 

Ih 0 Jo Jn 2d n 

Xa(6pk) Pk ( Wa,h{x) + Ul a ,h{y ))> 


i D _ 


r*l 


IC h {x - y) 


Ih 0 Jo Jn 2d h 

Xa( S Pk) PkWa,h(x), 


( Pk(Pk{x )) - P k {pk{y )) 




r*l rt 


Kh{x - y) 


Ih 0 Jo Jn 2d h 

Xa( S Pk)pkW a ,h(x). 


(A/j,P k (p k (x)) - A^P k (p k (y))) 


1-1) The term 7 0 . This term is handled just as in the proof of Lemmas 18.4118.51 bv 
using Lemma 18.51 and for this reason we do not fully detail all the steps here. First 
note that Lemma 18.31 applies to Dp k u k as well as for D p k u k as 

Dp k u k = [y k { - a^Ek)- 1 A) Dp k u k . 

Then as before, we first truncate by using some smooth function lfc(t,x,y) = 
</>(p k (t,x)/L)<j)(p k (t,y) / L) with some smooth and compactly supported function <j> 
leading to I 0 = Ifr + I^ L with 

I o = “b / f [ Kh ^ h —— {Dpu k (x) - Dpu k (y)) 

J Jho Jo Jn 2d " 

x'a(dpk) Ik Pk ( Wa,h{x) + W a , h (y)), 

and 

Io L = ~l[ [ [ —^^ (Dpu k {x) - Dpu k (y)) 

z Jh 0 Jo Jn 2d n 

Xa( S Pk ) (1 ~I k )Pk (U’a,h(x) + W ajh (y)). 

Remark that divu*, £ P k (p k ) £ L P J2 and since A M is an operator of 0 order, 
AftPkiPk) £ Lpx ■ Therefore by Eq. (18.201) 

sup \\D Pk U k \\ r mi„( 2 ,p/ T ) < 00 . 

b ^t,x 
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On the other hand \x'ai^Pk)\ <C{1+1) (\pk(x)\ e + \Pk{y)\ f ') and this lets us bound 
very simply I RL by Holder estimates 

Iq L <C{l+t)\ lOg/lol \\Dpk Wfe|| L mm( 2 ,p/ 7 ) ||(1 - lt)pl +i H L ma x ( 2 , 9 ) 

< c (1 + l) I logftol 11(1 - It) p\ +l \\ lTT ™, 

with l/g + 7 /p=l. But g (1 + i!) < p by the assumption p > 7 + 1 + £ and similarly 
2 (1 + 1) < p. As a consequence for some exponent 9\ > 0 

lt L <C(l+e)\\ogh o \L- 0 \ (8.21) 

We now use Lemma I5U1 for DpkUk and $ = x'a^Pk) it Pk w a,h{x). We note that 
||<f>|| L =o < C (1 + £) L 1+t . Moreover just as in the proof of Lemmas 18.4118.51 we can 
show that $ satisfy a transport equation giving that 


(7$ = 


/ n d 


I\ h (x - y)$(t,x,y)dy 


1,1 .T.-l , 1 


wr w. 


ln d 


I\ h (x - y)$(t,x,y) dx 




<C{l + £)L 


1 +e 


By Lemma I5T51 we obtain that for some 9 2 > 0 


It <C(l + i) L 1+e [ h 92 ^ < C (1 + I) L 1+e . 

Jh 0 h 


( 8 . 22 ) 


By optimizing in L, this lets us conclude that again for some 0 < 9 < 1 and provided 
that p > 7 + 1 +1 

I 0 <C(l + £)\logh 0 \ e . (8.23) 


1-2) The term I D . This term has the right sign as 

(y)) 

> C [ [ [ - —Xa( S Pk)tiPkpl(Wa,h(x)+W a!h (y)). 

Jh 0 do dn 2d n 

We will actually give a more precise control on I D + II D later on when the corre¬ 
sponding decomposition of II = IIq + II D + II R will be introduced. 


Kh(x - y) 


h 0 do dn 2d 


( pU X ) - Pk(y)) Xa( S Pk) Pk ( W a ,h(x ) + W a ,h 


1-3) The term I R . The difficulty is thus in this quantity. From its definition, A^ 
is a convolution operator. With a slight abuse of notation, we denote by A M as well 
its kernel or 


A^f = / A^x - y) f(y) dy, 

Jn d 
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and note that A^ corresponds to an operator of 0 order, i.e. it for instance satisfies 
the property J A A fx = 0 for any annulus A centered at the origin, |^4 M (a;)| < C \x\~ d . 
Decompose 

A tl = L h + R h , suppL/j C {|x| < d/j}, 

such that both Lh and Rh remain bounded on any L p space, 1 < p < oo, and 
moreover Rh is a regularization of A^ that is Rh = A M * N$ h for some smooth kernel 
N $ h . The scale Sh has to satisfy that 

5 h «h, log « | log/i|. 

Oh 

For simplicity we choose here 6h = h/\ log ft.|. 

Contribution of the Rh part. 

The first step is to decompose Rh into dyadic blocks in Fourier. Introduce a 
decomposition of identity \Eq as in sections 151 and fl3l s.t. 1 = ^ ’ll; and write 

Ilog 2 s h \ 

Rh = ^ * Rh + Rh, Rh — 'ki* Rh = Nh* Ng h * A^. (8.24) 

;=|log 2 h| 2<|log 2 h| 

Note that we of course require of the \Eq to satisfy all the assumptions specified in 
section m for the definition of Besov spaces. Define now Nh = Nh*Ns h , this kernel 
Nh therefore satisfies that for any s > 0 

WNhWw* <Ch~ s , (8.25) 

and moreover by the localization property of the T*,, one has that for s > 0 and any 

M < i 

[ \z\ s \N h {z) + N h (z + cor)\dr<Ch s . (8.26) 

Jn d 

Fix t for the moment and decompose accordingly 

f f K> h ^ II R h*pl(t,.) - R h *pl{t,. + z)\\ L i 
Jho Jn d n 

< f 1 [ ^^\\Rh*Pl(t,.)-R h *pl(t,. + z)\\ Li 

Jho Jn d h 

rl I l°g2 <5^. | . JT , . 

+ / / —1—||T ; * R h * pl(t, ■) - * Rh * Pk(ti • + z )\\lI- 

J ho j |, , | J n d ^ 

u 1 =I log 2 h I 

By (18.2511 and (18.2611 . the kernel Nh satisfies the assumptions of Lemma RTT1 and thus 
with Uk = Afj, * pf, applying Lemma HTTl for any q > 1 

f 1 [ ^±\\Rh*pl(t,.)-R h *pl(t,. + z)\\ Ll 

Jh 0 Jn d " 

= T f ^7T^ \\Nh * Uk(t, .)-N h * u k {. + z)\\li 

Jho Jn d h 

<C\logh 0 \ 1/2 \\U k (t,.)\\ Li . 
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Recalling that A p is continuous on every L p space, one has that ||£4 ||lp < C ||yOfc|| 
hence 


/ / II Rh*pl(t, ■) ~ Rh*pZ(t, ■ + z)\\ L i < C | logh 0 1 1/2 || p k {t, 

Jh 0 Jn d 11 

On the other hand simply by bounding 

\^>i*R h i,pl{x)-^ii,R h *pl(y)\ q <C\^ii,R h i,pl(x)\ (1 + \^ii,R h i,pl{y)\ q 
we write 

1 | log 2 Sh | 


i E 

ho Hlog 2 h\ Jnd 


K h (z) 


l * Rh * P 7 k {t, ^ * Rh* Pfc(L ■ + 2 )ll Lg 


<c E 

i<l log 2 ^ol+log2 I lo g2 M 


\$!i* R h * pl(t,.)\\ L * f 
J 2-1 


dh 


recalling that Sh = h/| log 2 h|. This leads to 


I log 2 i5h| 

E 


K h {z) 


I'lh *R h *pl(t,.) - 'L i * Rh* pl(t, ■ + z)\\li 


h ° i=| log 2 h\ n 

<C £ logl\\yi*Rh*pl{t, OIU®, 

l<2 | log 2 ho | 

and can in turn be directly bounded by 

< Clog | log h 0 1 \\^l* R h*pl{t,-)\\Lg 

l<2 | log 2 h 0 | 

< C log | logho| I log^o| 1/2 ||-ff/i *-)IIl| < C|logho| e \\pk{t,.)\\l„, 

with 0 < 9 < 1 by Lemma Tl3.31 Combining with the previous estimate, we deduce 
that 


rl r ~T 7 


K h (z) 


Ih 0 Jn d 


\\Rh*pl(t,-) -R h *pl(t,- + z)\\ L g 

< C | logh 0 | e \\pk(t, .)||2«7- 


(8.27) 


with 0 < 6 < 1. And therefore since x' a (£) — (1 + ■£) |£|^, by Holder’s inequality with 
the relation 1 /q + (1 + £)/(! + 7 + £) = 1, that is q = (1 + t + r y)/'y 


Kh(x - y) 


1 h 0 J 0 J n 2d 


(■ Rh*Pk(x ) - Rh*Pk(y)) Xa(fipk) Pk (w a ,hi x ) + W a,h(y)) 


>~C( 1+e) /V(*,-)llUE, f 1 f ^-\\R h *pl(t,.)-R h *pl{t,. + z)\\ Li . 

J 0 Lx Jho Jn d n 
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Finally by (18.2711 there exists 0 < 9 < 1 s.t. 


Kh{X , V) (R h *pl(x)-R h *pl(y)) 

Ih 0 Jo Jn 2d n 

Xa( S Pk) Pk ( W a ,h(x ) + W a , h (y)) 

>-C(l + £)\logho\ 6 \\p k \\ltlt+f 


(8.28) 


Contribution of the Lh part. It remains the term with Lh where we symmetrize the 
position of the weight with respect to the convolution with Lh by 

Iff — —( L h*pl(x)-L h -kpl(y))Xa(Sp k )p k Wa,h(x) 

Jh 0 Jo Jn 2d n 

= lh - Diff, 

with 

T h = f [ [ Kh ^ h ^ L h(z)(pl(x - z) - pl(y - z)) 

Jh 0 Jo Jn 3d n 

Xo (Spk) Pk wl~ n (x) wl{x - z), 

for n = 1 — l/y. Recall that since w a ,h = Kh*w a , \w a ,h{x) — w a ,h{x — z)\ < h _1 \z\ 
while \z\ ~ Sh on the support of Lh ■ Thus using that |xl(£)l — f rom <18.21) and 
that \Lh(z)\ < C\z\~ d , one has 


Diff = 



Xa^Pk) Pkwl h 


—^ L h (z)(pl(x -z)- pl{y- z)) 
(x) (w™ h {x) - w a t h(x - z)) 


< C{l+£) 

< c(i + e) 



-t+1 1 -n 


Using \z\ < Sh = j lo £ /t | , we obtain on the other hand that 

f 1 f <c f 1 Kdh = r 1 dh 

Jho hl+n J\z\<8 h \z\ d ~ n ~ Jho h±+n Jho h\logh\ 

<C I log hoi 1 -". 


As n = 1 — I/ 7 , this leads to 

Diff <C{ l+t) \\p k \\ltlt+i I Iog/i 0 | 1/7 . (8.29) 
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As for the first term by Holder inequality, using again that |x„(OI — C £\£,\ l 


Ih <C(l+£) 



\6pk\ e+1 % +1)/e 


a-e){e+i)/e 

W a,h 



W+ 1 ) 



( L h*((pk(-)~ Pk(-+ W )) 


<,(•))) 


(7+t+l)/7 

\x 



7/(7+f+1) K h (w) 

h 


dw, 


provided that l is chosen s.t. 


7 


j+£+l £+1 


= 1, 


or 


implying for instance that 


e + i 
i 


= 7, 


7 + ^ + 1 
7 


7 

7 - 1 ’ 


Given those algebraic relations and recalling that Lh * is continuous on every L p for 
any 1 < p < oo 



Since, using the definition of £, 

I Pk(x) - p\(x + w)\^ +e+1)/ ' 1 <7P 7 |<5pfc| (7_M+1)/7 =7 P 1 \5pkf +1 , 

one has 

Ih<C{ l + £)j [ [ ^0—1\Sp k \ e+1 p 7 k w a ,h{x)dxdw , (8.30) 

J h 0 J n 2d n 

which multiplied by —a^v^/2 will be bounded by I D + II D provided |a.^,| is small 
enough. 

II) The quantity II. Let us turn to II and decompose it as for / 

II = II 0 + II D + II R , 


where 


IIo = ~^>[ f f — —0pu k (x) + Dpu k {y)) 

1 Jh 0 Jo Jn 2d n 

(x'a( S Pk)5pk - 2Xa(8pk)) ( W a ,h ( x ) + w a,h(y)), 
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while 


n ° = r 1 r r K ^ x y) ( Pk ( Pk ( X )) + Pk ( Pk (y))) 

2 Jh 0 Jo Jn 2d h 

(Xa( b Pk) b Pk ~ 2 Xa(Sp k )) (W ath {x) + W a , h (y)), 

and 

nR = _a R Vk r 1 r r K h (x — y) {A p k{pk{x)) + A ,P k (p k (y))) 

2 Jh 0 Jo Jn 2d n 

( Xa( S Pk)Spk - 2 Xa(Spk)) ( Wa,h(x ) + W a , h {y)). 

II-1) Term IIq. For the term Ho, using Lemma 15751 in a manner identical to Iq 

IIo<~ f j [ —^ 7 — —K h *DpUh(x)(x' a {Spk)8pk-2xa(Spk))w ai h(x) 

J ho Jo J n 2d n 

+ C(l + £)\\p k \\ 1 +i e+1 \logh 0 \ e , 

for some 0 < 6 < 1. Using formula (15.111) or (18.201) . one has that 

div u k - a^A fl P k (p k ) > Dpu k , 
and hence since £ + 2 \a > -C (1 + €) \a 

no<Ci\\p k \\ l +i t+1 \\ogh 0 \ e 

+ C( 1+i) [[[ Kh ^ x - —(Kh*(\dWu k \+a tl .\A f j,Pk(pk)\)Xa{^Pk)wa,h(x), 

Jh 0 Jo Jn 2d n 

(8.31) 

and the first integral will be bounded by II a /2 for A large enough. 

II-2) Term II D . The term II D is controlled by I D : For a > b, by (18.21) 

( a 7 + 6 7 ) (~Xa( a - b )i a ~ b ) + 2 Xa{a - b)) > -(a 7 - b 1 ) ^-j^x' a ( a ~ b )( a + b )• 
Therefore 

1° + H D < —C 7 ^ f f — — —— \5p k \ e+1 pl(w h (x) +w h (y)), (8.32) 

2 Jo J n 2d P 

for some C independent of l and 7 . 

II-3) Term II R . The control on the last term, II R , requires the use of the penal¬ 
ization n a 

IjR + 7 n a < - a/Hfc / / / —^— -{AvPl(x)-A IJ K h i'pl(x)) 

2 Jh 0 Jo Ju 2d n 

(Xa( b Pk)bPk ~ 2 Xa(Spk))w a>h (x). 
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We use the same decomposition of A^ = Lh + Rh as for I R . 

Contribution of the Rh part. 

Note that, as Xa(0 < C |£| 1+f and |x(J < C (1 + £) |£| £ , for q = (1 + i + X )h or 
1/q + ( 1 + £)/(! +£ + 7 ) = 1 


- an Vk 


Khi - X y) {R h *pl(x)-K h *R h *pl(x)) 


Ih 0 JO JU 2d 


h 


(Xa( S Pk)Spk - 2 Xa(Sp k )) W a ,h{x) 


<c(i + o fwpkic^+U^ f j ^M\\R h ^ P i( t> .)~R h ^pi(t,. + z)\\ Lg . 
Jo Lx Jho Jn d h 


Now by estimate (18.271) . we have that 
rl rt r K h (x — y) 


- an is k 


h 0 Jo Jn 2d 


(R h * pl(x) - K h * R h * pl( x )) 

(x'a( S Pk)Sp k - 2 Xa(6p k ))w a ,h(x) 


< C(l + t) |log/i 0 | 3/4 ^ \\p k (t, ■WjtU+y \\Pk(t,.)\\ 
<C(l + t)\ log/i 0 | 3/4 Wpkf+ltZ- 


(8.33) 


7 

Ll q 


Contribution of the Lh part. Similarly as for I R 1 we symmetrize the weights leading 
to the following decomposition 

-a^Vk f f f — h ~, —— (L h * pl(x) - K h * L h * pl(x)) 

Jh 0 Jo Jn 2d n 

(. Xa( S Pk)Sp k - 2Xa(Spk)) Wa t h(x) 

= /4 L + Diff 2 , 


where 


IIh=a^Uk[ f f — ^— —Lh*{wl h (x)(pl(x)-K h *pl(x))) 

J ho Jo J n 2d n 

(-Xa($Pk)Spk + 2 X a( 5 pk)) wlfhix), 

with still n = 1 — I/ 7 . The term Diff 2 is controlled as the term Diff in I R using the 
regularity of w a ,h and yielding 

Diff 2 < C(1 + £) llog/jol 7 ||Pfc||l++^. (8.34) 


We handle lift with Holder estimates quite similar to the ones used for the 









term 1^, recalling that L k * is bounded on any L q space for 1 < q < oo 


Hh <Ca tl v k { 1 + t) 


I<h{w ) 


ho J o J n d 
\\\Spk\ e+1 wl-^(x)\\ Li t+iyi dw 
-1 r t r 17 


\< h (Pi(.)-K h *pim\L< + i 


<Ca^v k { l + t) [ [ [ w a,h(x) \p k (x) - K h * p k (x)\ e+1 dxdw 

J ho Jo Jn d " J n d 

f 1 r t 


+ Ca ll v k (l+t) f f I Kh j w i f Wah (x) \6p k \ (t+1)2/e dxdw. 
Jh 0 Jo Jn d n Jn d 


' Kq J 0 

One has immediately that 

r*l r t 


f f f Kh ^ w a , h (x)\5p k \^ +1) ! l dxdw 

Jh 0 Jo Jn 2d n 

K h {w) 


1 r t r 17 


< 

Ih 0 Jo Jn 2d 
as |<5pfc| < p k and again (£ + l)/£ = 7 . 
Finally as (£ + 1) (7 — 1) = 7 

rl '* ' K h (w) 


■ w a ,h(x ) |<5 p k \ f+1 pi dxdw 


Wa,h(x) I pl(x) - K h *pl(x)\ + dx dw 

lh 0 Jo Jn 2d 11 

< 7 [ [ [ Kh( u ) K h (z)w a ,h(x) \pk{x) ~ Pk{x + Z)\ e+1 ( pk(x ) + pk(x + z)) 1 
Jh 0 Jo J n 3d 


< 


7 [ [ [ Kh w a ,h{x) | pk(x) - Pk{x + z)\ 1+t ( p k (x ) + Ph{x + z )) 1 , 
Jhn Jo j U 2d n 


/ h 0 Jo Jn 2d 

as Kh(w ) is the only term depending on w and is of integral 1. Therefore 

Hh < Ca ll ,v k 'y(l + £) f f f w a ,h{x) \Sp k \ 1+t p\. (8.35) 

Jh 0 Jo Jn 2d P 


To conclude, we sum all the contributions, more precisely (18.231) . (18.281) . (18.291) . 
(I8.30D . (18.311) . (18.32L (18.331) . (18.341) . and (18.351) . to find for some 0 < 9 < 1 and 
provided p > 7 + 1 + £ 




K h {x - y) 


I h 0 J n 2d 

r-i 

< 

I h 0 Jn 2d 


h 

Kh{x - y) 


(w a ,h(x) + Wa,h(y)) Xa(Sp k ){t) 

(w a ,h( X ) +Wa,h{y)) Xa(dp k )\t =0 + C (1 + £) \ loghof 


+ C (a M u k (1 + i) - C 7 f f f Kh ^ w a , h {x) \5p k \ l+l p\. 
v 1 7 Jh 0 Jo Jn 2d h 


h 
































This finishes the proof of the lemma: As 1 + l = 7/(7 — 1) is bounded (we recall 
that 7 > d/2), if a M < C* for C* > 0 well chosen, the last term in the rdr.s. is non 
positive. □ 


8.6 Conclusion of the proof of Theorem 15.31 

We combine Lemmas 18.21 and 18.61 to get the following estimate 


n 1 n clJ~l 

/ / K h (x - y)(w a ,h{x) +w a , h {y))Xa(5pk)(t)dxdy— 

Jh 0 Jn 2d n 

< C | log h 0 \ 6 + initial value, 

with 0 < 9 < 1. We now follow the same steps as in the proof of Theorem 15.11 with 
the weight wo(t,x) + wo(t,y). We define u} v = {w a< h(t,x) < rf\ and note that 

r rl r 


/ n 2d 


£h 0 (x - y)xa{6p k ) = 


lh 0 J n 2d 

r l 


K h {x - y)Xa(Spk ) 


dh 


'ho Jx£u)C or 


Kh(x-y) Xa(Spk ) 


+ 


Now 


' ho J xGuJrj and y£uJn 


dh 


Kh{x - y)xa( 5 p k ) 


dh 

h 

dh 


ho J x£gj° or y^uj?. 


Kh{x - y) Xa(Spk) — 


1 

< - 


1 /» cUx 

/ K h (x — y)(w a ,h{x) +Wa, h {y))Xa{.dp k )— < C | log h 0 \ e , 
V J h 0 Jn 2d h 

while by point iii) in Prop. 17.21 and using the L p bound on p , for some 6 > 0 


/ 

1 J X 


ho Jx€(jj v and y^uj v 


K h (x - y) Xa{5p k ) ^ < 2 


< 


h 0 dn 2d 
C | log h 0 \ 

I log v\ e 


K h (x — y) pl +e l-K hirWa < v Y 


Hence we have 


/ n 2d 


fch. 0 (x - y)Xa(5p k )(t) < C | log h 0 \ 


log ho] 8 - 1 


I log 


< 


C\\ICh 0 \\ 


L 1 


I log I log /loll 6 


by optimizing in p and recalling that ||/Ch, 0 ||ii = | log/io|. 

Using Prop. 14.11 together with Lemma 16.61 one concludes that pk is locally 
compact in t, x. Thus we conclude the proof of Theorem 15.31 
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9 Proof of Theorems 13.11 and 13.21 : Approximate se¬ 
quences 

In this section, we construct approximate systems that allow to use Theorems 15.11 
and 15.31 to prove Theorems 13.11 and m 

We do not need to use here pressure laws P which depend explicitly on t or x , 
which simplifies the form of the assumptions on the behavior of P, either (15.81) or 

m • 

9.1 From regularized systems with added viscosity to no vis¬ 
cosity 

Our starting point for global existence is the following regularized system 

! d t pk + div(pkUk) = a k Apk, 

dt(pkUk) + di v(p k Uk < 8 > u k ) - pAuk - (A + p)Vdivu fc - A e *u k (9.1) 
+VP £ (pfc) + a k Vp k ■ Vuk = pkf , 

with the fixed initial data 

Pk\t—0 = P°, Pk u k \t=o = P° u°. (9.2) 

The pressure P e satisfies the bound (13.21) with 7 > 3d,/(d + 2) uniformly in e, that is 

C " 1 p 7 - C < P E {p) < Cp 7 + (7, 

implying that e(p) > C _1 p 7-1 — C. In addition we ask that P e satisfies the quasi¬ 
monotone property (15.81) but possibly depending on e, i.e. there exists po )£ s.t. 

(P £ (s)/s)' > 0 for all s > p 0 , £ . 

And finally we impose an e dependent bound (13.21) on P £ for some y E > d 

C- 1 P 7e -C< p £ (p) < C £ p 7s + C, (9.3) 

Similarly A e is assumed to be a given smooth function, possibly depending on e but 
such that the operator defined by 

T>ef = -pAf - (A + p)Vdiv/ -A e *f 

satisfies (12.21) and (12.31) uniformly in e. 

As usual the equation of continuity is regularized by means of an artificial viscos¬ 
ity term and the momentum balance is replaced by a Faedo-Galerkin approximation 
to eventually reduce the problem on X n , a finite-dimensional vector space of func¬ 
tions. 

This approximate system can then be solved by a standard procedure: The ve¬ 
locity u of the approximate momentum equation is looked as a fixed point of a 
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suitable integral operator. Then given u, the approximate continuity equation is 
solved directly by means of the standard theory of linear parabolic equations. This 
methodology concerning the compressible Navier-Stokes equations is well explained 
and described in the reference books [35], [3B], [53] . We omit the rest of this classical 
(but tedious) procedure and we assume that we have well posed and global weak 
solutions to m- 

We now use the classical energy estimates detailed in subsection 12.11 Note that 
they remain the same in spite of the added viscosity in the continuity equation 
because in particular of the added term a^V p k ■ Vuj in the momentum equation. 
Let us summarize the a priori estimates that are obtained 

sup sup / (pfc \u k \ 2 + pi) dx < oo, sup / / |Vitfc | 2 dx dt < oo, 

k,e t Jn d k,s Jo Ju d 

from m • The estimate (12.91) may actually require the e dependent bound from 
m with 7 e > d to control the additional term Vp*, • Vu k . It provides 



with p £ = 'j s + 2 r y e /d— 1 oip = "f + 2 r y/d— 1 which means p > 2 as 7 > 3d/(d+2). 
This bound may not be independent of e because it requires & However since 
a k vanishes at the limit, it still implies that any weak limit p of pk satisfies 



for p = 7 + 27 /d — 1 . 

From these a priori estimates, we obtain (15.41) and (15.51) . And from those bounds 
it is straightforward to deduce that pk Uk and pk \uk\ 2 belong to L\ x for some q > 1 , 
uniformly in k. Therefore using the continuity equation in (19.11) . we deduce (EH). 
Using the momentum equation, we obtain (15.61) but this bound (and only this bound) 
is not independent of e because of A e . 

Finally taking the divergence of the momentum equation and inverting A 

(A + 2 p) div Uk =P e (Pk) + A _1 div {dt{pk u k ) + div (p k u k < 8 > u k )) 

- A -1 div (p k f + A e * u k ) + «fcA _1 div (Vp fc • Vu k ), 

which is exactly (15.21) with p k = A + 2 p satisfying (1 5.3 p and compact, while 

F k = -A -1 div (p k f + A e *u k ) + afcA ^ 1 div(Vp fc ■ Vu fc ). 

The first term in F k is also compact in L 1 since A e is smooth for a fixed s. On the 
other hand 

afcA -1 div (Vp k ■ Viifc) 

converges to 0 in L 1 since y/ukA7p k is uniformly bounded in L 2 and Vu k is as well 
in L 2 . Therefore F k is compact in L 1 . We may hence apply point i) of Theorem 15. II 
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to obtain the compactness of pk in L 1 . Then extracting converging subsequences, 
we can pass to the limit in every term (see subsection 12.21 for instance) and obtain 
the existence of weak solutions to 

dtp + div(pu) = 0, 

dt(pu) + div(pit Cg u) — pAu — (A + /i)Vdivit — A e * u + VP e (p) = pf. 

(9.4) 


9.2 General pressure laws: End of proof (Theorem 13.lh 

Consider now a nonmonotone pressure P satisfying m and (13.31) . Let us fix co, e = 
1/e and define 

Pe{p)=P(p) if P<C 0 ,e, P e {p) = P(c 0 ,e) + C(p-C 0tE y, if p > Cq, s - 

If 7 < d then we also add to P e a term in £ p/ e to satisfy (19.31) . Note that P e is 
Lipschitz, converges uniformly to P on any compact interval. Due to (13.21) there 
exists po,e with po,e —> +oo as e —> +oo such that for p > po.ej 

(P e (s)/s)' = (P/(s)s - P(s))/s 2 > (C (7 - 1 ){p - c 0 , e ) 7 - P(c 0 iE ))/s 2 > 0. 

The approximate pressure P e still satisfies (13.21) with 7 and due to the previous 
inequality it satisfies (15.81) for p > p 0 e and (15.91) in the following sense: For all s > 0 

|P e '(s)| < Psi-%<co.' + C(7 - l)s 7_1 I s >c 0 . e - 

As a consequence, we have existence of weak solutions (p, u ) to (19.11) for this P e 
(assuming A e = 0 ) for any e > 0 . Consider a sequence £k —> 0 and the corresponding 
sequence ( pk , Wfc) of weak solutions to (19.11) . 

Because the previous a priori estimates were uniform in e for the limit p and u , 
(including (15.61) since A e = 0 ), then the sequence (pk, Uk ) satisfies all the bounds 
(15.41) . (15.51) . (15.61) . (15.71) and (15.91) . 

Moreover the representation (15.21) still holds with pk = 2p+X, compact in L 1 and 
satisfying (15.31) . Finally the exponent p in (15.41) can be chosen up to 7 + 27 /d — 1. 
Since 7 > 3 d/{d + 2) then p > 2 and since 7 > (7 + 1) d/{d + 2) then one has p > 7 
as well. 

Therefore all the assumptions of point ii) of Theorem 15.11 are satisfied and one 
has the compactness of pk- Extracting converging subsequences of pk and Uk, one 
passes to the limit in every term. Note in particular that P ek (pk ) converges in L 1 
to P{p) by the compactness of Pk- the uniform convergence of Peu to P on compact 
intervals and by truncating Peu (Pk ) for large values of pk since the exponent p in 
(El is strictly larger than 7 . 

This proves the global existence in Theorem 13.11 The regularity of p follows from 
Theorem 15.21 which concludes the proof of Theorem 13.11 
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9.3 Anisotropic viscosities: End of proof (Theorem 13.21) 

For simplicity, we take / = 0. Consider now a “quasi-monotone” pressure P satis¬ 
fying (EZD- Observe that P then automatically satisfies (13.21) since P(0) = 0. To 
satisfy (15751) . we have to modify P on an interval (co j£ , +oo) with co i£ —>• +oo when 
£ —> +oo. More precisely we consider P e as defined in the previous subsection 

Pe(p) = P{p) ifp<co, e , P £ W = P(co, E ) + C(p-c 0 , E ) 7 . ifp>c 0 , £ . 


Remark that since 7 > d here, we never need to add a term with 7 £ . 
Now given any smooth kernel, for instance K , we define 


A E k u = div (6A(t) V K e * u). 


Because of the smallness assumption on SA(t), the operator T> e satisfies (12.21) and 
m uniformly in e. Therefore we have existence of global weak solutions to (liHl) 
with this choice of P E and A E . We again consider a sequence of such solutions 
( pk , Mfe) corresponding to some sequence e k —»• 0. Because the estimates are uniform 
in e for ( 19 . 41 ) . we have again that this sequence satisfies the bounds ( 15 . 41 ) . ( 15 . 51 ) . ( 15 . 71 ) . 
We now assume that 




implying that p in (15.41) is strictly larger than 7 2 /(q — 1). Moreover observe that 


ll-4 Efc Mfc || L 2 h -i < C \\Vu k \\ Llx , 

such that HMD is also satisfied. 

For simplicity, we assume that 6A has a vanishing trace (otherwise just add the 
corresponding trace to p). Denote = 2 ||< 5 A||/(2 p + X) and the operator E k 


E k u = —div 


SA(t) 

2\\6A\\ l - 


V K e -k u 



SA i:j (t) 

2 ||<W|| L - 


(1/7 k 'll . 


For small enough, A — a^Ek is a uniform elliptic operator so that (A — E k )~ x A 
is bounded on every L p space, uniformly in k. For the same reason, A^ = (A — 
a^Ek )- 1 E k is bounded on every L p space with norm less than 1 and can be repre¬ 
sented by a convolution with a singular integral. 

Taking the divergence of the momentum equation in (19.11) , one has 

(2/r + A) (A div Uk - a M E k divu) = A P(p k ) + div ( d t (p k u k ) + div ( p k u k ® u k ))■ 

Just write A P E {p k ) = (A — a M Ek) P E {pk) + E k P E , take the inverse of A — E k 
to obtain 

(2 p + A) div u =P(pk ) + (A - a p P fe ) _1 E k P{p k ) 

+ (A - E k ) _1 div ( d t {pk u k ) + div ( p k u k <8> u k )), 


94 









which is exactly (15.111) with v k = (2 p + A) -1 . As a consequence, if < C*, which 
is implied by ||<L4 ||lc*> small enough, then Theorem 15.31 applies and pk is compact. 
Passing to the limit again in every term proves Theorem 13.21 Note that P Sk (pk) 
converges in L 1 to P(p) for the same reason than in the previous subsection. 

The case with D(u ) instead of Vrt. Let us finish this proof by remarking on the 
different structure in the case with symmetric stress tensor div (A Dm). In that 
case, one cannot find divu^, by taking the divergence of the momentum equation but 
instead we have to consider the whole momentum equation. Let us write it as 

£ u k = VP(pfc) + d t (pk u k ) + div (p k u k <8> u k ), 

with £ the elliptic vector-valued operator 

£ u = p A u + (p + A) Vdiv u + div ( 8A D u). 

The operator £ is invertible for 8A small enough as one can readily check in Fourier 
for instance where —£ becomes a perturbation of p |£| 2 I+(p+ A) Its inverse has 

most of the usual properties of inverses of scalar elliptic operator (with the exception 
of the maximum principle for instance). Therefore, one may still write 

div u k = div £~ 1 VP(p k ) + div £~ 1 ( d t (p k u k ) + div ( p k u k ® u k )), 

leading to the variant (15.121) of the simpler formula (15.111) . 


10 Extension to the Navier—Stokes-Fourier system 

We present here a direct application of our compactness results to compressible 
Navier-Stokes systems with temperature. Those systems are considered to be more 
physically realistic. They also offer many examples of non monotone pressure laws 
thus being especially relevant to our approach without any thermodynamic stability 
assumption. 

This chapter is in many respects a preliminary work, with non optimal assump¬ 
tions. It is likely that further works would require estimates specifically tailored to 
the case under consideration, and for the same reason would require to specifically 
tailor the compactness Theorem 15. II But we believe that the result shown here can 
provide a good general basis or first step for such future works. 

10.1 Model and estimates 

The heat-conducting compressible Navier-Stokes equations read 
( d t p + div (pit) = 0, 

< d t (pu) + div(pu (g) u) — div<S + VP(p,i?) = 0, (10.1) 

[ d t (pE) + di v(puE) + div(P(p, d)u) = div(<S-u) + div(ft($)V$), 
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where E = |w| 2 /2+e is the total energy with P = P(p, d) and e = e(p, d) respectively 
stand for the pressure and the (specific) internal energy. As usual the system is 
supplemented with initial conditions. 

We denote the stress tensor T> = div S and as in the barotropic case, we always 
assume that for some function p 


Vu : S dx 


p |Vu| 2 dx. 


Instead of the temperature d, one could also choose as a third unknown the internal 
energy (or the entropy as defined below). Some formulas are easier when using the 
set of variables (p,d) and some with the set (p, e). For this reason and simplicity 
we here follow the classical notations of Thermodynamics for partial derivatives, 
denoting for instance d$f\ p if / is the function f(p,d) and d e f\ P if instead one 
considers the composition f{p,d(p, e)). 

The total mass as well as the total energy of the system are constants of motion 
namely 

J p(t, •) dx = J po dx 

J{\p\ u ? + pe{p,d)){t,-)dx = J (^poKI 2 + p 0 e{po,d 0 ))dx. 


and 


From the mathematical viewpoint, the heat-conducting compressible Navier-Stokes 
system suffers the deficiency of strong a priori bounds. For instance in comparison 
with the previous barotropic Navier-Stokes, the energy bound does not yield any 
L 2 H], bound on the velocity. 

In order to be consistent with the second principle of Thermodynamics which 
implies the existence of the entropy as a closed differential form in the energy balance, 
the following compatibility condition, called ’’Maxwell equation” between P and e 
has to be satisfied 


p=y|2 

op 


a dP 
■& od 


1 de 


and 

ds 


1 dP 

P = ddd 

P 

dp 


p 2 dd 


The specific entropy s = s(p,d) is defined up to an additive constant by 

ds 

dd 

An other important assumption on the entropy function is made, 

the entropy s is a concave function of (p -1 ,e). 

This property ensures in particular the non negativity of the so-called C v coefficient 
given by 


r - — 

v dd 


1 d 2 s 
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If (p, d) are smooth and bounded from below away from zero and if the velocity 
field is smooth, then the total energy balance can be replaced by the thermal energy 
balance 


C v p(d t d + u • V0) - div(/e(0)V0) =S :\7u- ^ dP { p ^ divu. 

ov 

Furthermore, dividing by we arrive at the entropy equation 

dt(ps) + div(psu) -div^--- J = ^{S : Vu + —-— J. 


( 10 . 2 ) 


10.2 The entropy estimate through thermodynamical stabil¬ 
ity 

The first difficulty in the Navier-Stokes-Fourier system is to obtain a L 2 estimate 
on Vm, since, contrary to the barotropic case, it does not follow from the dissipation 
of energy. 

Instead a first approach is to use the entropy. Coupling the equation on the total 
energy and the equation on s, we get 


J(^P\u \ 2 + p[e(p,d)-d ir s])(t,-)dx +J j -y(S:Vu 
= J(\po\ u o? + Po[e(po,d 0 ) - d*so]) dx, 


+ 


c|Vi?| s 


(10.3) 


for any constant temperature $*. This is exactly the type of modified energy estimate 
that we are looking for as it provides a control on n(d) \X7d\ 2 /d 2 and with any 
reasonable choice of stress tensor, a control on |Vu| 2 /i9. 

However it is useless unless one can bound 


# 0 * (p, d) = p[e(p, d) - $*s] 

Observe that (using the mass conservation) 


- / = - / \H^(p,d)- 


dt 


dt , 


dH^(p*>^*) 

dp 


(p- P*) - HoAp*,#*) 


Remark that, using the Maxwell equation and the definition of the entropy, we have 

d 2 H^{p,d*) l&P(p,i?*) 


dH^{p,d) _ p .deippd) 

dd / dd ’ 


dp 2 


dp 


Those estimates were developed in the works by E. Feireisl and collaborators, 
under the assumption of thermodynamic stability 


de{p,d) 

dd 


> 0 , 


dP(p,d) 

dp 


> 0 . 


(10.4) 
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The meaning of such condition is that both the specific heat at constant volume C v 
and the compressibility of the fluid are positive. However as we have already seen, 
the latter condition is violated by several physical state law, such as the standard 
Van der Waals equation of state. 

Under (110.4p . one obtains that is increasing in d for d > 0* and decreasing 
for d < that is that has a minimum at d = $*. The second part of (110.41) 
implies that H#* is convex in p. 

One chooses accordingly p* as the minimum of ff,? t (p, d+) and decompose the 
quantity linked to H-s (p, d) into two parts 


HvAp,#) - P+) ~ #tf,(p*,0*) 

= [HvAPi®) - h dAp^*)\ 

+ [H«. (p, 0*) - ^)] > 0 


(10.5) 


by the thermodynamic stability assumptions on de(p,d)/dd and dP(p,d)/dp. As a 
consequence under (110.41) . one deduces that 




0 


d 2 


dx dt < constant. 


( 10 . 6 ) 


10.3 The pressure laws covered by previous works. 

In every previous work, the viscous stress tensor is assumed to be isotropic 

S = p,(Vu + Vii T ) + A div u Id, 

with coefficients p, A either constant or depending only on d. 

1) The pressure law as a perturbation of the barotropic case. It is due to E. Feireisl 
who considered pressure laws under the form 

P(p,d) = P c (p) + dPo(p), 

where 

P c { 0) = 0, P' c {p) > oip 7_1 - b for p > 0, 

P c (p) < 02 p 7 + b for all p > 0, 

P#( 0) = 0, P'${p) > 0 for all p > 0, ^ 

P#(p) < c(l+p^), 

and 

7 > d/2, ft < ^ for d = 2, ft = ^ for d = 3 

Z O 
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with constants ai > 0, 02 , b and P c , P$ in C[0, + 00 ) flC 1 (0, + 00 ). In agreement with 
Maxwell law and the entropy definition, it implies the following form on the internal 
energy 

e (P,d) = J ^^-ds + Q(d), 

where Q'(d) — C v (d) (specific heat at constant volume). The entropy is given by 

s(p,d) = f Cv ^ ds - S#(p), 
dp* s 

where S$ is the thermal pressure potential given through S#(p) = / P$(s)/s 2 ds. 

dp* 

The heat conductivity coefficient k is assumed to satisfy 


+ 1) < k(i?) < K 2 (d a + 1) for all d > 0, 

with constants K\ > 0 and a > 2. The thermal energy Q = Q(d) = C v (z)dz 
has not yet been determined and is assumed to satisfy / 2g [ 0 +OD ) C v (z) > 0 and 

C v (d) < c(l + d a / 2 ~ 1 ). Because the energy and pressure satisfy (110.41) . the estimate 
on Hi) gives a control on p 7 in L°°{ 0, T ; T 1 (fl)) and through (110.61) a control on d in 
L 2 (0, T; L 6 (fl)) in dimension 3 and in L 2 (0, T; L p (fl)) for all p < +00 in dimension 
2 . 

Because (110.61) does not provide an H] : bound on u, E. Feireisl combines it with 
a direct energy estimate (see below). Therefore one obtains the exact equivalent of 
estimates (12.61) in this case. Using similar techniques, one then proves (12.91) for 
0 < a < mm(l/d,2d/ r y — 1). The proof of compactness and existence follows the 
main steps shown in subsection 12.21 


2) Pressure laws with large radiative contribution. It is due to E. Feireisl and 
A. Novotny who consider pressure laws exhibiting both coercivity of type p 7 and 
d 4 for large densities and temperatures namely 


P(p, d) = ^/(7-l)Q ( _g_ ) + with > 0) 7 > 3/2, 


with 

Q G C^fO.+oo)), Q(0) = 0, Q'(Z) > 0 for all Z > 0, 


and 


lim 

Z —^-|-oo 


Q(z) 

Z T 


— Qoo ^ 0 . 


In agreement to Maxwell law and the definition of entropy, it implies the following 
form on the internal energy and the entropy 


e(p,d) 


1 # 7 /( 7 -!) ^ p ^ ^4 

(7 — 1 ) p l? 1 /! T -1 ) p 
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and 




They impose 


4a d 3 
3 p ' 


0 < -S'(Z) 


1 7 Q(Z)-Q'(Z)Z 

7-1 Z 


< c < +oo for all Z > 0, 


with limz^+oo S(Z) = 0 so that thermodynamical stability (110.41) holds. Therefore 
the energy provides uniform bounds in Lf*L\ for t? 4 and p 1 . One assumes in this 
case that the viscosities and heat conductivity satisfy 


fi,Ae C 1 ([0, +oo)) are Lipschitz and p (1 T d) < p(d), 0 < A($), po > 0, 

and 

k € C 1 ([0, Too), kq(1 T'd 3 ) < /«($) < «i(l T i? 3 ), 0 < kq < K\. 


Almost everywhere convergence of the temperature is obtained using the radiation 
term. Extra integrability on P(p , d) can be derived just as in the barotropic case. 
Finally the same procedure as in the barotropic case is followed to have compactness 
on the density, relying heavily on the monotonicity of the pressure dP(p, , d)/dp > 0. 
This gives global existence of weak solutions (in a sense that we precise later). 


With respect to these previous works, we focus here, as in the barotropic case, 
in removing the assumption of monotonicity on the pressure law, which is consistent 
with the discussion of such laws in subsection 12.41 


10.4 The direct entropy estimate 

We explain here the general framework for our result on the Navier-Stokes-Fourier 
system. The estimates here closely follow the ones pioneered by P.-L. Lions, and 
E. Feireisl and A. Novotny; just as in the barotropic case, our contribution is the 
new compactness argument not the energy estimates. With respect to the previous 
discussion, we only present them here in a more general context as in particular we 
will not need the monotonicity of P. 

If one removes the monotonicity assumption on P then thermodynamic stability 
(IK).41) does not hold anymore. Following P. L. Lions, it is however possible to 
obtain the estimate (jlO.fil) directly by integrating the entropy equation (110.21) 

T Ac('d) ^ ) dxdt < C J ps(t,x) dx. 

Therefore the bound (110.61) holds under the general assumption that there exists C 

s.t. 

s(p, i?) < C e(p, •d) + —. (10.8) 

P 
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Recall that 


and 


e = m(d) + / 
J p 


p ' 2 

d$P 


dp, 


d P s = — 

P 

We also have that d$s = d^e/d, therefore as long as m(d) > 0 with 


m'(s) 




ds < (7(1 + m(d)), 


and 


r^ iC+c 

Jo- P Jo- P 


/ p* p' J p* p 

for some C > 0 then (110.81) is automatically satisfied and one obtains the bound 
(110.61) . Moreover if e(d, p) > p 1 ~ 1 /C then one also has that p £ L/° L/. Assuming 
now that 

/ci (d a + 1) < k(0) < n 2 (iT + 1), 

with a > 2, one deduces from (110.61) that log d £ L d H], and d a ^ 2 £ L\ H], or by 
Sobolev embedding d £ Lf ^ 2 ^ d \ By a Holder estimate, it is also possible to 
obtain a Sobolev-like, L Pl W^ P2 , bound on u 


|Vu| P2 dx 


Pl/P2 


dt < 


|Vu| 5 


Pi/2 


dx dt 


$P2/(2- P2 ) dx 


Pi(2-p 2 )/(P2(2-pi)) 


(2—Pl)/2 


dt 


< 00 , 


provided that j? 2/(2 — P 2 ) = cc/(l — 2/d) and pi/(2 — pi) = a or 


2a 2 ad 

1 + a’ d(a + l) —2 


(10.9) 


Unfortunately this Sobolev estimate does not allow to derive an equivalent of (12.91) . 
Instead one requires a L 2 H^ estimate on u (the critical point is in fact the L\ with 
value in some Sobolev in x ). Instead one can easily extend the argument by E. 
Feireisl and A. Novotny: For any one can write 


1 d 

2 dt 


P M 2 + j t J <t>(p) + j $ :Vu = J(P(d, p) - <j>'{p) p + <j){p)) divu. 


This leads to the assumption that there exists some <f> s.t. 

C^p 1 -C< <t>{p) < Cp 1 + C, 

\P(d,p) - ct>’{p) p + <f>(p) I < c(> + 0* + y/rffrp )), 
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with 


( 10 . 11 ) 




Indeed, with (jlO.lOjl . one has 


5 


: S7udx dt < E(p°,u°,d°) + C J J (^p^ 1 + + y/pe/d, p)^ |divu| dx dt. 

By (12.21) . this leads to 

f f \S/u\ 2 dxdt < CE(p 0 ,u 0 ,d°) + C\\\7u\\ L 2 H/ 1 + 0 /Ja || i 2 (10.12) 

Jo J t,x 

and the desired H 1 bound follows from (110.111) . It is now possible to follow the same 
steps to obtain an equivalent of (12.91) if 7 > d/2 

f ( p^ +a dxdt<C{T,E{p°,vP,d 0 )), for all a < 1/d (10.13) 

Jo Jn 

Note here that the assumptions (110.10D - (|10.11I) are likely not optimal and will be im¬ 
proved in future works. They nevertheless already cover the two examples presented 
above. 


10.5 Main result in the heat-conducting case 

For convenience, we repeat here all the assumptions presented above and the con¬ 
cluding mathematical result which may be obtained. Let P{p,d) a positive pressure 
law such that that for some C > 0 and 7 > d 


P(p, d) such that - f P dp' <C + C f - — dp', 
Jp- p Jp * p 


e(p,d) = m/d) + j 
with m/d) > 0 


(p(p',d)-mp( P ',d)) 


dp 1 > 


p- 1 

~c~ 


~c~p' 


and 


[ ds < c(l + m/d)) <C{l + d a (7+«- i )/2(7+o)) ) (10-14) 

J#* S 


for 


s 

Ki (d a + 1) < n(d) < «2 (d a + 1), n, A constant and a > 4, 

c - V - c < (t>{ P ) <c P ^ + c, 

\P(p, d) - 4>'(p) p + tip) I KCpP'+Cd^+C y/M^dj, 

. \doP(p,d)\ < Cp^ + Cd^ 




^3 < 

2 


7 + a + 1 


< 


a 


2 ’ 


(10.15) 


d 


p + A > 0, 7 ^ > 0 
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where we recall that a < min (1/d, 27 /d — 1) = 1/d since 7 > d here. We also assume 
that the specific heat is positive (as is necessary for the physics) i.e. 


C v = d#e(p, d) > 0, Vp, d, 
and that the pressure contains a radiative part 

d$P(p = 0,d) > 0. 

We do not need to impose monotony on P and it is enough that 
\d p P(p,0)\ < Cp 1 - 1 + Cd b with 0 < 6 < 
Finally the initial data has to satisfy 


po e i 7 (n d ), d 0 ei w ( u d ) 

with po > 0, do > 0 in II d 


and 



> 0 , 


and 

E 0 = [ +p 0 e(p 0 ,d 0 )) <+ 00 . 

Ju d ^ 2 po 2 


(10.16) 

(10.17) 

(10.18) 

(10.19) 

( 10 . 20 ) 


Then, on the condition of approximate solutions construction, we can prove the 
following result. 

Theorem 10.1 Assume that (110.141) (110.181) are satisfied with 7 > d and that the 
initial data satisfy (110.191) and (110.201) . Then there exists 

p e L°°([ 0 , T], L 7 (n d )) nL 7 +a ([ 0 , T] x n d ), for all a < 1/d, 
ugl 2 {[ 0, t], 7 d 1 (n d ))nL°°([o, t], x 2 (n d )), (10.21) 

d e L a ([o, T], ^“/C 1 - 2 /^ (n rf )), iogdeL 2 ([o, T], ^(ir*)), 

with d > 0 a. a. on (0,T) x II d that is solution in the sense of distributions to 

d t p + div(pu) = 0 , , 

( 10 . 22 ) 

dt(pu) + div(pn ® u) — p A u — (A + p)V divu + VP(p, d) = 0, 


while fde equation on the temperature is satisfied in the following sense: In the sense 
of distribution, the entropy solves 


dt{p s(p, d)) + div (p s(p, d) u) — div 


c(d) Vd 

_ d 


>i(5:Vn- 


|Vd | 5 


, (10.23) 


with S = p, (Vn + Vti T ) + A drew Id and recalling that 


J p* P J 0 u 
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This is supplemented by the total energy property 

J(p^f + pe(p,d))dx = j (p°^y- + p° e(p°, d 0 )) dx. (10.24) 

and the initial data conditions satisfied by (p, pit, ps) in a weak sense 

p|t=o = Po, Ht=o = (P u ) o, Ht=o+ > Pos(Po,^o)- 

Remark 10.2 An example of pressure law included in Theorem, \10.1\ is a perturba¬ 
tion of the truncated virial expansion as described in subsection \2.f\ plus a radiative 
part namely 

N 

P(p,i?) = p^+tf^R„(i?)p" 

n—0 

with 7 > 2N > 4. Choosing m = constant for simplicity in this example, this leads 
to 

7—1 N n —1 1 

e(p, d) =m +-—- - ^d 2 B' n { d) -—- - d 2 B[(d) logp + d 2 B' 0 (d)~, 

7 -i ^ «-l P 

For simplicity, let us assume that B\ = constant, which is the normal virial assump¬ 
tion, see [35J for example, so that this term vanishes. The entropy reads 

' r.n—1 1 

s(p, 0) = - £> B' n fd) + B n (d)) - + B x logp + (dB' 0 (d) + B 0 (d))~ 

^2 n P 


1) We assume that the pressure contains a radiative part, namely that Bq is convex in 
d with C~ l fib- 1 < B 0 (d) < i? 6 " 1 and C” 1 d b ~ 2 < B' 0 { d) < d b ~ 2 where 2 <b< a/2. 
This already satisfies (110.171) . 


2) For n > 2 , the coefficients B n can have any sign but we require a concavity 
assumption: 




This ensures, with the assumption on Bg, that the specific heat C v = d$e(p,d) 
satisfies (110.161) . This is again a classical assumption for the virial, see still [35] . 
Note that it would be enough to ask this concavity of some of them and moreover 
that this is automatically satisfied if B n ~ —d v , that is precisely for the coefficients 
contributing to the non-monotony of P in p. 


3) We also require some specific bounds on the B n namely that there exist B n and 
£ > 0 S.t. 


IdB’M + lB^^Cd^b- 1 -^ 

|B„(0)| + \dB n {d) - B n | < (10.25) 
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First of all this gives us a bound from below on e 


e(p, d) >m 


> m + 


P 


.7 — 1 


7- 1 
, 7-1 


,b * 

+ c -^^-ydd z 


p 


n" 

■b-l-e P 


n =2 

b N 


+c- i — -cy[ p'l- 1 - 5 ' 

7 — 1 p ' 


n — 1 

0 6 ~ e '' 


n=2 


P 


by Young’s inequality, so that this implies (110.141) 9. Assumption (110.181) is proved 
with an identical calculation. The same calculation also proves Assumption (110.14l) i 
by showing that s < C (p 7_1 + + 1). 


4) Then choosing 


4>{p) = 


P’ 


_ _„ n 11 

y B n -^~ 

^ n -1 


7-1 


Bn 


2<n<N 

and using again the second part of (110.251) we have that 


\p-<mp)p+m\<c y \vB n {i))-B n \p n <c y d^ i -^- E P n 

n<N n<N 

< ClV(p 7/2 + ^“ /2 - e ), 


still by Young’s inequality. This yields (110. 141) ^7 with the right inequalities on /3\ and 
P 2 in (110.151) . The same calculation also proves that |3^P| < C (p 7,/2+e + ■0 a / 2 ~ e ) 
thus proving (110.14I) k with bounds (110.151) on ^ 3 , / 34 . 

We will treat more general virial expansions in a future work. Note for a fixed d 
then P(p,’d) is indeed increasing with respect to p after a critical p$ which depends 
on d and can be arbitrarily large where d » 1. This is the reason why P does not 
satisfy any of the classical monotonicity assumption such as m or (EH) and why 
our new approach is needed. Our pressure law has two important parts: the radiative 
term (corresponding to n = 0) to get compactness on the temperature and the 4>{p) 
term to get compactness on the density in time and space. 

Remark 10.3 The assumption (110.181) is more demanding in the Navier-Stokes- 
Fourier setting (though satisfied by the example just above). It would be more natural 
to have instead 

\d P P(p, 0)| < P(d) p 7 " 1 + C{d), (10.26) 

for some unbounded function P of d. This should be possible and require some 
modifications in the proof of Theorem 15.11 combined with appropriate estimates on 
d. This will be the subject of the forthcoming work HU. In general obtaining more 
optimal results for the Navier-Stokes-Fourier system seems to depend on adapting 
precisely the estimates and Theorem 15.11 to the specific model under consideration. 

Remark 10.4 With appropriate additional growth assumption, it is possible to show 
that the combination of (110.231) and (110.241) implies the usual energy equation (and 
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not just an inequation as (110.231) . The reason for the formulation (110.231) - (|10.241) is 
that it cannot be proved without those additional assumptions. We refer to [35] and 

m- 

Remark 10.5 In Theorem 1 10. 11 the viscosity coefficients p, A are independent of 
the temperature d. Instead several models use temperature dependent coefficients 
p(d), \(d). To handle that case, the proof given below would have to be modified; 
the compactness of the temperature would have to be established first. 

Sketch of the proof. The proof of Theorem 110.11 follows the same steps as the 
proof of Theorem 13.11 For this reason we only sketch briefly the points which are 
similar and insist more on the differences induced by the presence of the temperature. 
The general method is taken from [35] and is adapted to take advantage of our new 
compactness estimates. As usual one starts from an approximate model 

d t pk + di v(pkUk) = a k A p k , 

dt(pkUk) + di v(p k Uk 0 u k ) - pAu k - (A + /u)Vdiv« t + VP E ,s(pk,d k ) 

+ a k Vp k ■ Vu fc = 0, 

Kr(dk) 

dt(Pk s E (pk,dk )) + div (p k s E (pk,d k ) «fc) - div ( £ \7d k ) (10.27) 

+ a k ^f^(d k s e {pk,dk) - e E (p k ,dk) - P ^ Pk,dk ^ ^ 

v k Pk 

with again S k = p(S7u k + Vm^) + Adi vu k Id and where s ak G M 1 ([0,T] x II d ) 
vaguely converges to a], s G M 1 ([0,T] x II d ). The entropy actually only satisfies an 
inequality since a\ s ak is unknown and s e is hence not fully defined from (110.271) . 
However we supplement this by imposing the conservation of the total energy 

/ P (^T + e e,s{P^)j dx = J P° + e ^s(p°,d°) S j dx. 

Note that the last quantity in the entropy relation was already introduced to com¬ 
pensate some terms in the energy equality. Imposing the conservation of energy is a 
simple way of fixing the remaining additional terms (which are all positive) instead 
of writing them explicitly. 

Remark that contrary to the barotropic case, we always need here two steps of 
approximations in the pressure (and hence energy), indexed by e and <5. The existence 
of such solutions may be done adapting the procedure for instance that we can find 
in [36] using appropriate Galerkin scheme, penalization terms and precise estimates. 
This important part is out of the scope of our paper but will be fully described with 
more complicated pressure law in [18] (replacing the assumption (110.181) by (110.261 ) ). 

The first level of approximation is obtained by adding a barotropic correction 
P s (p,d) = P(p,d) + 5(p^ + p 2 ). (10.28) 
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From the approximated pressure, one defines as usual the approximate energy 


es(p,$) =m(i?)+ f 

J O 




n r 2 


dp' 




7-1 


P^ + Sp, 


(10.29) 


and the entropy also has the straightforward expression 


s,(e,<,)= f^ds- r = (1 o.30) 

J 0 s J p* P 


In the second stage of approximation, indexed by e, we first replace «;($) by 
K e ( , d). We assume that k e = k + k E with k E a positive function which vanishes for 
2 e<d<l/2 e and satisfies 


7t e (i9) —> +oo asil->eor£ ', J k E (d) dfl < oo. (10.31) 

This is slightly different from [36j but gives similar bounds on the temperature while 
being simpler to handle in our case. Then we define P e as an approximate pressure 
law of Ps truncated for large densities, namely for p E = e^ 1 , we define 

Pe{p,$) = P(p,fi) if P<Pe, Pe(p,ti) = P(Pe,d) if P > Pe (10.32) 

Note that since the temperature will be bounded by (110.3111 . it is not necessary to 
truncate it in (110.3211 . This leads to 


Pe,s(pA) = p s(pA) +< 5 ( p 7 +p 2 )- 

We also use p E to truncate the initial density p°: We define p E = p° if p° < p E and 
p E = p E otherwise. We impose p E1 u°,d° as initial data for (110.2711 . 

The energy density is still defined from the pressure P Ey s or P e , 


e e (io,d) = m(' i?) + f 

dp' 

e s ,s(p,&) = + f 

J p* 


(P e (p'A)-ddoP s (p'A) J , 
-- d P > 


(Pe,8{p',&) -ftd#P e Ap'A) ,, 
-35- d P 


= e s (p,$) H-7 P 1 1 +Sp 

7-1 

and so is the entropy which is actually independent of 6 (just as ss is independent 
of <5 in (110.301) 1 


SE,s{p,d) = s £ {p,&) = [ m ^ d s - f 
J o s J 0 


r p doP E (p',d) J , 

-72- d P ■ 

p 


(10.33) 
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Remark that for p < p e then e s (p, d) = e(p , d) and similarly s e (p, d) = s(p, d). This 
trivially implies in particular that e e (p®,d°) > 0 and that 

[ e E (p°, d°) dx —> f p° e(p° , d°) dx as e — > 0 , 

Jn d “ “ J n d 

which enables to propagate the energy. At later times, one may not have pk < p E 
and e E (pfc,'dfc) may be negative at some points. Nevertheless from (110.14 I) k. one has 
that for p > p e 

,9^4 + 1 

\e e (p,d) - e e {p s ,d)\ <Cd^ + C -. (10.34) 

P P 


The first step is to pass to the limit otk —> 0 in the system (110.2711 to obtain a 
global weak solution to 

dtp + div(pit) = 0 , 

d t (pu) + di v(pu ® u) — pAu — (A + /x)Vdivtt + = 0, 

'n e {d) 


d t (ps e (p,d)) + div (ps £ (p,d) u) - div 




(10.35) 


= vis 


1 


<S : Vit 


k e |W | 5 


where < 7 ^ G M+([0,T] x II d ) with a total mass bounded uniformly with respect to 
e and 8. The system is still supplemented by p®, u°, d° as initial data. 

Note that the entropy still only satisfies an inequality, given the unknown positive 
measure, however this is supplemented by the preservation of total energy 

/ P (^T + e e,s(p,^) S j dx = J p° +e Ej 5 ( / 9 °,'d 0 )^ dx. 

As we remarked above, this identity occurs also for ctk fixed as mentioned in [36]; 
this is the reason for the appropriate corrections in cufc in the momentum and entropy 
equations. 

Then one considers a sequence £> —> 0 and the corresponding sequences pk, Uk , 
of solution to (110.351) . Passing again to the limit, we obtain global weak solutions to 


dtp + di v{pu) = 0 , 

d t {pu) + di v(pu ®u) — p A u — (A + ^)Vdivw + \/Ps(p, d) = 0, 


d t {p s(p, d)) + div (p s(p, d) u) - div 


fn{d) Vd\ 

V # J 

a s + ^ ^5 : Vzi + 


k |Vd | 2 \ 

* )' 


(10.36) 


with <rf S M + ([0,T] x 12) with total mass bounded uniformly with respect to 6 and 
with now p°, u°, d° as initial data.. We still conserve the energy equality 

J p + dx = J p ° (yr - + e < 5 (P°^°)^) dx - 
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Finally, we take a sequence <5*, —>• 0 and obtain global weak solutions as announced 
by Theorem llO.il 


The three limits are somewhat similar with slightly different a priori estimates 
and in the case of (110.271) to (110.351) additional difficulties in obtaining those a priori 
estimates. This limit uses point i) of Theorem 15.11 and the other two point ii). For 
this reason, we sketch the limit (110.271) to (110.351) separately and both limits, (110.351) 
to (110.361) and (110.361) to (I10.22D - (I10.23I) at the same time. 

The limit from (110.271) to (110.351) . We use the a priori estimates described above, 
in particular in subsection 110.41 The entropy estimate though requires more care 
for System (110.271) . We emphasize here that most of the complications in this limit 
were already present in [36] and are solved here in a very similar manner. The novel 
features of our work are mostly present in the other two limits, from (110.351) to 
(110.361) and from (110.361) to the limit system. 

However because we handle pressure terms that are more general, we can unfor¬ 
tunately not simply use the result of to skip this first limit. Instead one has to 
check carefully that indeed at this level, estimates work the same. 

First of all, given assumption (110.311) . standard parabolic estimates imply that 
any solution to (110.271) is either constant in x or satisfies 

e < d k (t,x) < (10.37) 

£ 


Remark that since 


k E < oo, we do not have strict inequalities in (110.371) . In 


addition, the conservation of energy yields 


sup / pj.it, x) dx < oo. 

t,k Jn d 


(10.38) 


It also useful to write a partial energy estimate, based only on the barotropic part 
of the potential energy 


Pk 


¥+4 2+ ^t 


dx 




dx 


7 - 1 , 

J J {P e {p k ,d k )d\\u k - S k : Vu fc ) 

-ot k 5 [ ( (2 + 7Pfe“ 2 ) |Vp fc | 2 . 


(10.39) 
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Since P e < C e (8k is bounded and P e is truncated in pk), we deduce that 

\u k \ 2 


Pk 


+ 6 P 2 + 


1 

Pk 

7-1 


dx 


^ C e Oik & 


f j (2 + 7 pT 2 ) |vp fc | 2 . 


(10.40) 


Now observe that 
A p k 


a {pk s e e e ) 

Vk Pk 


= div 


(8 k s E - e e - Vpfc 


Pk '&k -I 


+ 


|Vp fc | 2 

Pk 8 k 


(10.41) 


{pe + pd p e e )^ 


V p k ■ Vd fc 


where the functions P e , e £ and s e are all taken on the points pk(t,x), 8 k (t,x). The 
identity (110.411) is obtained through the general relations between s, e and P and in 
particular 

<9p p p dp 


and 


d 


P 


dd S 8 p8 


1 / de 

) = #l e + V' 


Integrating in space and time the entropy equation in (110.271) . we therefore get 


< Oik 


n d 8k 

( 


0 Jn d 


+ 

~\~ Pk 


A8k) 


|Vd fc | ; 


8k 

v p k ■ vd fc 

#1 


Ot k 


0 7n d 


, p IVpfcl 2 
p e Pfctffc 


(10.42) 


+ / PfcS E (t) - / (pfcs £ )(0). 

in d Jn d 


Therefore for ak fixed we have extra terms in the right-hand side with respect to the 
case of System (110.351) or (110.361) as described in subsection 110.41 namely 


Oik I (&e 4“ Pk dpO-e') 


V pk ■ w fe 


/n d 


and an additional term on the left-hand side 


81 


Oik 


hi d 


dp \j£hl 

P 5 8 k ' 


On the other hand by the bounds (110.371) on 8k and the truncation in P £ imply that 
s E is bounded by some constant depending on s, hence trivially in that case 



Pk Se{Pki8k) dx C £ . 
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From the definition (110.321) of P E , we obtain that 


d p P e (p,i9) = 0, Vp>p E , Ve < d k (t,x) < 


(10.43) 


From the bounds (110.371) on d k , we deduce that P E satisfies 

|Vp fe | 2 . 


fa o |Vp fc | 2 ^ ^ /' 

&k I Q — C OLk I 

Jn d V k Jp k <Pe 


On the other hand, by the dissipation (110. 401) . we get a control on 

a k 6 f [ |V Pfc | 2 < C e 
Jo Jn d 


and thus 


/' n p |Vp *| 2 ^ r 

I Op± } £ ^ Cs,e- 

Jn d v k 


(10.44) 


Turning to the last term in the entropy estimate, we remark that e E and pd p e E are 
bounded still from the truncation of P £ , so that one obtains 


I (&e T Pk (Jp&e) 

Jn d 


^Pk ■ Vtffc 


#1 


<a k C E / Vp fc ■ Vd fc 
Jn d 


< 4 /2 v / i v ^i 2 + a * /2 [ W 

L Ju d J n d 


(10.45) 


, l/2 n 
— ^k 

Combining (I10.44D - (I10.45I) with (110.421) concludes the bound on the entropy. 

The rest of the a priori estimates described in subsection 1 10. 41 follow in a straight¬ 
forward manner. We hence summarize the uniform bounds here 

£ < Vk < 

£ 


sup f (pl + Pk w\ 2 + tfk) dx + a k 6 [ f \S7 p k f 
t,k Jn d Jo Jn d 

sup / f (p k +a + |V?ife| 2 ) dxdt < oo, a < 1/d, 

k Jo Jn d 

sup / f (|W“ /2 | 2 + |Vlogd fe | 2 ) dxdt< oo, 
k Jo Jn d ' ' 

sup \\dt(Pk s £ {p k ,'dk))\\M 1 ([o,T]-,w- 1 - 1 (n d )) < 
k 


< OO, 


(10.46) 
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Combining those bounds with the continuity and momentum equations shows that 
( 15 . 41 ) . ( 15 . 51 ) . ( 15 . 61 ) . and ( 15 . 71 ) are satisfied. Moreover taking the divergence of the 
momentum equation in ( 110 . 271 ) and inverting the Laplacian, one obtains 

(A + 2 n) divwfc =P e ,s{Pk^k) + A -1 div ( d t (p k u k ) + div ( p k u k ®u k )) 

+ a k A -1 div (Vp fe • Vu k ). 

This is exactly the identity ( 15 . 21 ) with P k (p : t,x) = P et $(p 1 'd k (t,x)) and p k = A + 2p 
which satisfies m • Observe that as a consequence if r is strictly larger than p E 
then Pfc(f, x, p) is only 5{p 1 + p 2 ) plus a given function of thus satisfying the first 
part of ( 15 . 81 ) . In addition if r and s are less than p E 


\Pk{t, x, r) — P k (t, y, s)| < max 

#fc(i,2/)], p<p, 


(\d p P e , s (p,d)\\r - s\ 
\doPe,s{p, I?) 11 •dk{t,x) - 'd k {t 1 y)\ s j 


< C E (|r - s| + - ’d k {t 1 y )|). 


Since i) k belongs uniformly to an appropriate Sobolev space, that means that the 
last part of (15.81) is also satisfied. We may hence apply the variant of Theorem 15.11 
yielding the compactness of p k in L [ . 

Since u k € L 2 H x , we get the compactness on ^ [p k u k in L 2 L 2 from the momentum 
equation. The passage to the limit in the continuity equation is therefore as usual. 
We also have compactness in space for u k and $k respectively from the viscosity and 
conductivity. 

The next step is to obtain the compactness of the temperature. This relies on 
the entropy inequation in (110.851) but follows the now classical approach and for this 
reason we only sketch the procedure: 

• Observe that p k s E (p k ,^k) in some L^ x with p > 1 from (110.461) . In fact 
PkSeipkjflk) is in L^°L], uniformly using the definition of s E (p k , r d k )• 

• Observe that d t (p k s E (pk,^k)) is uniformly in M 1 ([0, T]; W _1,1 (II d )). Rewrite 
the entropy equation as 


dt{p k s e (p k ,-d k )) +div(Rl k ) +R 2 Eik = <rl Aak , 

where R\ k , R 2 k contains all the other terms. By our previous estimates and 
in particular from (110.461) . we know that ||i?* is uniformly bounded in 

k for i = 1,2. Hence integrating the previous equation 



hi d 


Pk &£ (Pk i^k) 


R 2 


o Jn d 


e.k 


I [pk &£ (pk i ^/c)] \t—0 5: 


/n d 


Recalling that er l $ ak >0 then this implies that \\al s ak || M i < C E . Finally 
this bounds d t (p k s e {p k , tffc)) in M^ x + C 
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• Obtain the a.e. convergence of pk s E (pk, Ok)(t, x). From the first point and the 
compactness in space of Ok and pk, we deduce that pk s e {pk,Ok) is compact in 
space. From the second point, we deduce that it is compact in time. There¬ 
fore after extraction pk s E (pk,Ok) converges strongly in L\ x and after possibly 
further extraction, one has that there exists l(t,x) s.t. 


Pk{t,x) s s (pk(t,x),’&k{t,x)) — >l(t,x), for a.e. t, x. 


• Obtain the p a.e. convergence of Ok{t,x). Extracting again, we have from 
the compactness of pk(t,x), that pk(t,x) converges a.e. to p(t,x). By the 
assumption on the specific heat (110.161) . s E is invertible in Ok (P E and so s E 
is not truncated in 0 only in p). Therefore this implies that, p a.e., Ok(t,x) 
converges. 

• Obtain the a.e. (not just p a.e.) convergence of Ok- The integral defining s 
and s E is singular at p = 0 (because of the p'~ 2 factor). In fact ps(p, $)| P =o = 
d$P{p = 0,i?). Assumption (110.171) on the radiative part of P guarantees that 
ps is still invertible in 0 even if p = 0 and hence that Ok converges a.e. even 
on the set where p vanishes. 

• Dominated convergence then implies the compactness of 'Ok in L p ([0, T] x II d ) 
for any p < oo (recall that for fixed e, we have Ok uniformly bounded in time 
and space). 

This allows us to pass to the limit in every term of the momentum equation, including 
the barotropic term 5 p 7 using the extra-integrability pk G L' 1+as in (110.461) . We can 
also simply pass to the limit in the entropy s E (pk , Ok)- 

Define R e by k' s = k e . By (110.311) . we know that R e (iO) is bounded on [e, e -1 ]. 
This lets us pass to the limit in K e (0k)V0k = \7K e (Ok) by the compactness of Ok- 
Let us pass to the limit next in 



By relation (110.411) . the bounds on s E , e E , P El the a priori estimates in (110.451) and 
(110.461) and the dissipation term in (110.401) . this term will converge to 0 if we can 
prove that 



On the other hand from the continuity equation, one has that 


/ 



pk log pk (ti x) dx 


L 
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By the compactness of p k , we can pass to the limit in every term in the l.h.s. so 


/ p log p(t, x) dx — / p° log p°(x) dx + / / p div u dx dt 

Jn d J n d Jo Jn d 


= — lim a k 


r r ivpfcl 2 

/o Jn d Pk 


dx dt. 


On the other hand, p solves the continuity equation without any diffusion and since 
it belongs to L 2 by (110.461) . one also has that 

/ p log p(t, x) dx — / p° log p° (x) dx + / / p div udxdt = 0 , 

J n d Jn d Jo J n d 

which has for consequence the required property 


limafc [ ( -— ^—dxdt = 0. 

Jo Jn d Pk 

Note now that is not possible to pass to the limit in the r.h.s. of the entropy 
equation (I10.27IU without the use of appropriate defect measure. For instance the 
r.h.s. contains p |Vu/c | 2 which does not in general converge to p |Vit| 2 , as this would 
require the compactness of Vn*, in L 2 . Instead one uses convexity properties and 
compactness on the temperature dt to prove that 


— lim ■ 


Sk ^ Vllfc + 


|V0* 


= cr, 


£,<5 


+ 


S : V'U 


|V^»p 


where d 2 s £ M+([0, T] x fl). This leads to the equality in the entropy equation in 
(UnSSDs with a positive (and a priori unknown) measure cr 2 s which is the sum of 
a 2 s and the vague limit of a\ s Qfc when ak tends to zero. Note we can also write 
an inequality ignoring the positive measure. Finally we conclude the proof of this 
limit by noticing that we may pass to the strong limit in all the terms of the energy 
conservation. Therefore we indeed obtain the equality in (110.241) . 

Now use (110.14l) i (which implies (110.81) 1 to write the inequalities 

s(p,d) < C[e(p,d) + P~ l ], s e (p,d) < C[e e {p,tf) + f e {p)\, 


where f e {p) = 1 /p E + 1/p. 

We can hence integrate in t, x the entropy equations to find 

/ / < sup / p k s{p k ,d k ) < C sup / (pke{p k ,d k ) + 1) < C, 

Jo Jn d t J n d t Jn d 

f f cr 2 g < sup [ pk s e {pk, dk) (10.47) 

Jo Jn d t Ju d 

< C sup / (pk e E {pk, dk) + 1 + Pk) < C. 
t Ju d 
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This ensures that a'l 5 , erf G M 1 ([0,T] x Q) uniformly with respect to £ and <5. We 
again refer to [331 EH] for the details of similar procedures. 

It only remains to pass to the limit in the energy equality: For all t 

(pk e s ,6(pk, &k) + Pk ^y-^ dx = (p° ^APe > d°) + Pe yy) dx - 


It is straightforward to show the convergence of Pk^eAPkAk) to pe £ j{p,d) hi 
Lj oc t x , from the compactness of pk, 'Ok, the a priori estimates and the bounds 
on e Sj s that can be derived from (110.141b . (110.181) (see (110.341) for instance). Simi¬ 
larly the compactness of pk , Pk Uk and the a priori estimates imply the convergence 
of Pk \uk \ 2 in L\ t . Since this convergence is only L 1 in time and not L°°, we need 
one more step: Extract another subsequence s.t. for a.e. t 



^Pk {.Pki dfc) T Pk 





s{p, d ) + p- 


dx , 


which finally gives the conservation of energy a.e. t since the left-hand side is con¬ 
stant. 


The limits (110.351) to (110.361) and (110.361) to the final system. First we consider a 
sequence £k —> 0 and corresponding sequences of solutions pk, Uk, dk to (110.351) . 
Some of the a priori estimates for this sequence of solutions to (110.351) are obtained 
by simply keeping only the estimates in (110.461) uniform in e. The others, such as the 
entropy estimate, have to be derived again. But one proceeds exactly as described 
in subsection 110.41 leading to 


sup [ (pi + p k \u k \ 2 + dx < oo, 
t,k Jn d 

sup / / ( Pk +as + |37«fc| 2 ) dxdt < oo for all as < ~ v 

k Jo Jn d d 

„t „ 

7 log| 2 ^ dxdt< oo, 


(10.48) 


a/2,2 


k Jo J n d 
Similarly for (110.361) . we obtain 


8 ->■ 0 < 


sup / (pi + Pk \u k \ 2 + dff) dx < oo, 
t,k J n d 

sup / f (Pk +a + 137u/-) 2 ) dxdt < oo for all a < —, (10.49) 

k Jo Jn d d 


sup / / (|Vd ^ 2 | 2 + |V log^fc| 2 ') dxdt< oo. 

„ k Jo Jn d 3 > 

The bounds (110.491) will of course imply (110.211) after passage to the limit. As before 
combining those bounds with the continuity and momentum equations shows that 
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(l5~4l) . (1531) . (EH), and EH are satisfied with p = 7 + as or p = 7 + a for System 
(110.351) and for System (110.361) . 


We now have to check that the pressure defined by P k {t,x,pk) = Pi'dkit, x), pu) 
if p < e^T 1 or P k (t,x,pk) = P{dk{t,x),£ ^T 1 ) if p > e// 1 satisfies the stability assump¬ 
tions, namely EH - 

First property. First of all by the bounds in (110.141) 

\P{’d k (x), p k (x)) - P(’d k {y),Pk(y))\ < Qk(t,x,y) + I 


with 

Qk{t,x,y) = C{p k {x ) 03 + Pk{yf 3 + + {®k{t,y)) 0i ) | •&k{t,x) - i?k(i, 2 /)| 

and 

I = \P{’d k {t,x),p k (x)) - P(d k {t,x),p k {y))\. 

Note that 1 ? is compact in x in L\ x and uniformly bounded in LfLx ^ 1 2 ^ d \ Since 
Pa < a — 1 < a /2 as a > 2 , one has that for some 8 < 1 

[ [ K h (x- y) (i? fe (t,a;)) /34 + (^(t, y))^ 4 ) \$k(t,x) - &k{t,y)\ 

J 0 J u 2d 

<C\\K h \\ e L1 ( [ T [ I< h (x-y)\Mt,x)-Mt,y)\ 

\J 0 Jn 2d 

Similarly p k is uniformly bounded in L/~/. a so that for example 

as @3 a/(a — 1) < 7 + a. Indeed @3 < (7 + a + l)/2 and a > 4. This implies that 

sup/ f K »r\\ ^ Qk(t,x,y) —s> 0, as h —>■ 0. 
k Jo Jn 2d ll^hlU 1 

Note that if d p Pk(t , x, s) = 0 if s > ej“ 1 . Hence we can use (110.181) on to write 

| P(iJk{t,x),p{x)) - P(tf k (t,x),p k (y))\ ^Cipkixy - 1 + Pfc(y) 7-1 + iJk(t,x) b ) 

I Pk{x) - p k {y )|- 

Denote P k (t, x) = , dk{t,x) b . Let us now use again that d is compact in x in L] x 
and uniformly bounded in L/Lx ^ 1 and the assumption on b < a /2 in (110.181) 
to get the convergence 

sup/ / ^ r-Zi ^ I Pk(t,x) - P k {t,y) 1 -> 0 as h 0 , 

k Jo Jn 2d II-MIl 1 !! 
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together with the uniform bound 


sup||P fc || L 2 < sup||'d fc (t,a ;) b || I 2 < sup \\d k {t, x)\\ L? < oo. 
k ’ k t ' x k 


Second property. Following the same procedure as in the previous bounds, we obtain 
from (110.141) and (110.181) 


Pk{t,x,p k (t,x)) < C (p k (t,xy + (‘d k {t,x)) b p k (t,x) + ('d k {t,x)) 

+ p k {x) 03 

By Young’s inequalities with 1 /y* + l/y = 1 

P k {t,x,p k (t,x)) < C (p k {t, X y + {’d k (t,x)) b ' y * + 0 d k (t,x )) 1+ > 34 

Define 


1+04 


R k (t,x,y ) = 0? fc (t,:r)) b7 * + x)) 1+Pi + {■d k {t,x)) l/( - 1 

and observe that R k is compact in x in L 1 since d k is compact in x and b 7 * < a by 
(IIP.181) . 1 + /?4 < a since a > 2 , with finally 7/(7 — ^ 3 ) < a as fa < (7 + a + l)/ 2 , 
a < 2"f/d— 1 and a > 4 > d. 

Thus we conclude that the pressure law satisfies m for both limits with 7 = 7 . 
We can again apply Theorem 15.11 point ii in this case. We again need that p > 7 . 
Therefore we obtain the compactness of p k and from the momentum equation, the 
compactness of ypkUk- As in the limit of (110.271) to (110.351) . the next step is to 
obtain the compactness in time of d k - We follow the same procedure which however 
now requires more work on the first two points: 

• Bound pk s ek (pk,d k ) or Pk s(p k , , dk) uniformly in some Lf x with p > 1. This 
is not immediate anymore since s is not bounded and s Ek is not uniformly 
bounded. One uses the expression of the entropy in (110.331) and the bounds on 
\d$P\ and to in (110.141) to find 

Pk S Ek (p k ,Vk)<C {pi 3 + + Pk 1 ? “(7+«-l)/(7+a) )) (10 .50) 

and similarly for p k s. Then the a priori estimates on p k £ L 1+a or p k £ L 2+as 
and d k £ L^ x £\LyLiy are enough provided /?3 <7 + 0 and /I 4 < a or @4 < 7 ^; 
both are ensured by the stronger condition in (110.151) . 

• Bound d t {p k s ek {pk,’d k )), d t (p k s{p k , d k )) uniformly in M 1 ([0, T], W _1,1 (II d )). 
One first has to bound p k s ek u k in L\ x . For instance 

/ [ pi 3 u k < f f [p k \u k \ 2 + pi 03 " 1 ), 

J 0 Jn d Jo Jn d 
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which is bounded provided 2f3 3 — 1 <7 + 0 and 2 /3 3 — 1 < 2 + 0,5 leading to the 
assumption on j3 3 in (110.1511 . As for f 9^9 u k , we have by Sobolev embedding 
that Uk G 2 ) and therefore need that 9^ G L 2 L 2 d A d + 2 )_ Given 

that 9 k G PI L 1 ® , this is ensured by the condition on /? 4 (note that 

the condition could in fact be improved by using that 9 k G 
The last term is treated in a similar manner. From the entropy dissipation 
Ke(,$k) |V0fc| 2 /0* is uniformly in Lj^ as well. Therefore, we just have now to 
use that <r 2 5 is uniformly in M 1 ([0, T] x 0) proved previously through relation 
(110.471) . Thus we conclude the bound on d t (pk s Ek (pfc, $&))• Similar calculations 
may be done concerning dt(pk s(pk,9k)) with the defect measures appearing 
from the limit £*, —> 0. 

To conclude and pass to the limit in every term of the momentum equation, we 
use the integrability properties on p*,, 9 k (110.141) .1. (110.14k . (110.14k and (110.181) to 
provide an integrability with p > 1 of the pressure and more precisely 

[ \P(Pk,#k)\<C [ U k p^+pl+4p k )<CL- e , 

Jp k >LOXti k >L Jpk>L or 1 +>L v y 

for some 9 > 0 since p k is uniformly in L 1+a , 9k in L a and from the assumptions 
(110.151) on the /3j. Since P e ^ contains the additional barotropic part 6 ( pl+pl) which 
is also controlled, one has that 

[ (\Pe,s(Pk,9k)\ + \Ps(pk,9k)\) < C L~ e . 

JPk>L or 1+>L 

This lets us pass to the limit in Pg and P e ^g by using the pointwise convergence of 
Pk and 9 k - 

The type same control applies to pk s Et s(p k ,9k) or pk sg(pk,9k) so that we may 
pass to the limit in the entropy equation. The right-hand side requires the use of 
convexity just as for the limit of System (110.271) to (110.351) and the properties of the 
total mass defect measures control follow the same calculations than (110.471) . 

11 Models occurring in other contexts 

Macroscopic models in various biological settings involve a density p that is trans¬ 
ported by a velocity vector field u with source term, such as 

d t p + di v(pu) = p G(P(p), c) 

for some functions G and P. The function G may include birth and death terms and 
it could also depend on other quantities such as nutrients concentration denoted c, for 
example Oxygen in cancer modeling. This concentration is typically governed by a 
parabolic equation with right-hand side modeling the consumption of the resource(s) 

d t c— Ac = —pH(P(p),c) (11.1) 
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with 


d P H < 0, d c H > 0, H(P, 0) = 0. 

The velocity field is described through a constitutive law for instance 

- vA^ + ctf = P(p) - S, (11.2) 

where S' is a given source term. 

There have been several studies of such systems with applications to crowd mo¬ 
tion, traffic jams, cancerology using specific reformulations: gradient flow or kinetic 
descriptions and appropriate choices for / and G. 

For instance the special case c = 0, a = 0, G = 0 has been studied by B. Maury, 
A. Roudneff-Chupin and F. Santambrogio (see [54]) through the framework of 
optimal transportation reformulating the problem as a gradient flow in the Wasser- 
stein space of measures. Other examples concern the reformulation through a ki¬ 
netic formulation. For instance, very recently, B. Perthame and N. Vauchelet 
( see m) have studied the case c = 0 in the whole space with the pressure law 
P(p) = (7 + l)p 7 /7 with 7 > 1 , G satisfying 

G € C 1 (M), G'(-) < —r) < 0 , G(Pm) = 0 for some Pm > 0 

and a,v > 0 . 

The main result of this last paper is the ’’stiff pressure law” limit, namely the limit 
7 -A + 00 , leading to a free boundary model which generalizes the classical Hele-Shaw 
equation. Such kind of limit has also been performed for the compressible Navier- 
Stokes equations by P.-L. Lions and N. Masmoudi (see [52]) with P(p) = ap 1 
with G = 0. 

Recently C. Perrin and E. Zatorska (see [5B]) have studied the singular limit 
e —> 0 for a singular pressure law P(p) = ep 1 / (1 — p) 13 with 7 , /3 > 3. The advantage 
of such pressure law is that 0 < p < 1 for a fixed e which is important for some 
applications as mentioned by B. Maury in his review paper 153] . 

Pressure laws which blows-up for a critical density are of course the exact analo¬ 
gous of the Van der Waals equation of state for compressible fluid dynamics. They 
are also encountered in other setting such as crowd motion, granular flow, sedimen¬ 
tation problems. 

An other possibility to describe the velocity field is to consider the Brinkman 
equation instead of Eq. (111.21) . namely 

—vA u + au + \7P(p) = S , 

or with a Stokes viscosity term 

—vAu — (A + ^)Vdivu + au + VP(p) = S. 

This type of correction accounts for flow through medium where the grains of the 
media are porous themselves and has been justified in [2]. If the velocity is ir- 
rotational then this model is exactly reduced to Eq. (111.21) . Note however that 
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viscoelastic models for tumor growth may allow for instance to observe a lemon-like 
shape tumor, whereas with a Newtonian model an ovoid is obtained. Even limited, 
this difference can eventually lead to bigger ones as the outer rim is composed of 
proliferating cells with exponential growth. The kind of shape obtained for instance 
in m is observed in in-vitro experiments. 

What our method can bring: 

We do not try to state a theorem here given the large variety of possible models. 
Instead we give a few elements for which the method introduced here could prove 
crucial 

• More complex pressure laws, attractive and repulsive, could be considered. This 
would be the exact equivalent of Theorem 13.11 Note that biological systems 
frequently exhibit preferred ranged of densities for instance with attractive 
interactions for low densities and repulsive at higher ones. 

• More importantly the transition from attractive to repulsive interactions may 
depend on the concentration c of nutrients or other bio-cliemicals. This is 
similar to the dependence on the temperature in the state laws for the Navier- 
Stokes-Fourier system. For example if the pressure blows-up at some thresh¬ 
olds, enforcing a maximal density, then this threshold and the maximal density 
will depend in general on c. Because c is not necessarily uniformly bounded, 
the range of attractive interactions (where G(P(p),c) is decreasing in p) is not 
compactly supported and classical approaches may fail. 

• There can be several nutrients or bio-chemicals. That means that in general 
one has several ci, C 2 , ... with several equations (111.11) (or a vector-valued one 
if the diffusion speeds are the same). If chemotaxis is considered, some of those 
bio-chemicals may be attractive while other are repulsive. This may lead to a 
complicated pattern of interactions which again cannot be handled by classical 
approaches. 

• Many of this models are posed in porous media which are inherently anisotropic. 
In biology for instance the tissue or the porous matrix is heterogeneous. There¬ 
fore the equation for u should read 

—div(A(t, x)D(u )) + a(t, x)u + X7P(p) = S 

for some A and a and with in general a non-monotone pressure law. Our new 
approach could for instance help to enrich the model mathematically studied 
recently in [29] . 

12 Appendix: Notations 

For the reader convenience, we repeat and summarize here some of our main nota¬ 
tions. 

Physical quantities. 
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• p(t,x), or pk(t,x) denotes the density of the fluid. 

• u(t,x), or iik(t,x) denotes the velocity field of the fluid. 

• P(.), or Pfc(.) denotes the pressure law. 

• e(p) is the internal energy density. 

In the barotropic case, e(p) = j p p P(s)/s 2 ds. 

• E(p,u) = Jp( | u\ 2 /2 + e(p)) is the total energy of the fluid. 

• p, A and pk denote various viscosity coefficients or combination thereof. 

• S denotes the viscous stress tensor. 

In the simplest isotropic case: S = 2 pD{u) + Adivuld. 

• D is the diffusion term related to the viscous stress tensor by T>u = div<S. 

In the Navier-Stokes-Fourier case: 

we have the additional notations, appearing only in section [TOl 

• d(t, x) is the temperature field of the fluid. 

• s(p, d) is the entropy of the fluid. 

• n(d) is the heat conductivity coefficient. 

• C v is the specific heat of the fluid. 

Technical notations. 

• d is the dimension of space. 

• k as an index always denotes the index of a sequence. 

• h and ho are scaling parameters used to measure oscillations of certain quan¬ 
tities such as the density. 

• Kh is a convolution kernel on II d , 

Kh(x) = (h + |x|) _a for x small enough and with a > d. 

• K h is equal to K h /\\K h \\ L i. 

• ICh 0 = fh Kh(x) ^ is the weighted average of K h . Note that ||/C||x,i ~ | log /j-oI - 

• wq, w\ and w a are the weights and Wi^ = Kh * Wi their regularization with 
i = 0,1, a. 

• C is a constant whose exact value may change from one line to another but 
which is always independent of k, h or other scaling parameters. 

• e(h) is a smooth function with e( 0 ) = 0 . 

• 9 is an exponent whose exact value may change as for C but in (0,1). 

• The exponent p is most of the time such that p £ x . 

• q and r are other exponents for L p type spaces that are used when needed. 
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• I, II ,... and A, B, D, E... are notations for some intermediary quantities 
used in the proofs. Their definitions may change from one proof to another. 

• x, y, w, z are typically variables of integration over the space domain. 

• 6pk{x,y) = pk(x ) — pk{y) the difference of densities. 

• pk{x,y) = pk(x) + p k {y) the sum of densities. 

• D puk = A~ x div ( dt{pk Uk) + div (pk Uk <8> Uk)) denotes the effective flux. 

• The individual weights w(t,x) = wq, uq, w a satisfy Eq. (17.21) or 

dtW + Uk ■ Vw = — D w + afc A x w, 

where D = D 0 , D±, D a are the penalizations in (17.31) . (17.41) . (17.51) . 

• The weights ivq or w a may be convolved to give Wh = Kh*wo, w a ,h = Kh*w a - 

• The weights are then added or multiplied to obtain the composed W(t, x, y) = 
Wo, W u W 2 , W a with 

W 0 = w 0 (x) + w 0 (y), Wi=wi(x)+wi(y), 

W 2 = w 1 (x)w 1 (y), W a = w a (x) +w a {y). 

The main properties of the weights are given in Prop. 17.21 

13 Appendix: Besov spaces and Littlewood-Paley 
decomposition 

We only recall some basic definitions and properties of Besov spaces for use in Lemma 
16.31 We start with the classical Littlewood-Paley decomposition and refer to the 
readers for instance to and 0 for details and applications to fluid mechanic. 

Choose any family £ 5(13 d ) s.t. 

• Its Fourier transform is positive and compactly supported in the annulus 

{ 2 k ~ 1 < |£| < 2 k+1 }. 

• It leads to a decomposition of the identity in the sense that there exists $ with 
$ compactly supported in {|£| < 2} s.t. for any £ 

1 = 6 ( 0 + £**(£)■ 

fc>l 

• The family is localized in H d in the sense that for all s > 0 

sup ||^fc|| L i < 00 , sup 2 ks f \z\ s \ty k (z)\dz < oo. 

k k JH d 
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Note that in ]R d , one usually takes ^k(x) = 2 kd Af(2 k x) but in the torus, it can be 
advantageous to use a more general family. It is still necessary to take it smooth 
enough for the third assumption to be satisfied (it is for instance the difference 
between the Dirichlet and Fejer kernels). 

For simplicity, we then denote 'I'o = $ for k = 0 and for k > 1, = 

2 kd 5'( 2~ k x). For any / € <S'(K. d ), we also write fk = ^k * f and then obtain the 
decomposition 

OO 

/ = £/*■ ( 13 - 1 ) 

fc=o 

From this decomposition one may easily define the Besov spaces 


Definition 13.1 The Besov space B q)q is the space of all f G L\ oc fl S'{R d ) for 
which 


= 2 


s k 


Ukhih = 5Z 2 


1/9 


s k q 


WMh 


< OO. 


\k =0 


The main properties of the Littlewood-Paley decomposition that we use in this article 
can be summarized as 


Proposition 13.2 For any 1 < p < oo and any s, there exists C > 0 s.t. for any 
f€L} oc ns'(R d ) 


^\\fk\\L, <w' 2 fk\\ L , <cr k \\fk\\ L ,, 



( o° \ !/ 2 


/ °o \ !/ 2 

C - 1 

E 2 “'IAN 

<\\f\\w‘»<C 



\k~0 ) 

LP 

\k=0 ) 


And as a consequence for 1 < p <2 


C 


-l 


I Bp,2 ^ 


| w s ’ p — 


p<C I 


I B* 


In particular a consequence of Prop. 113.21 is the following bound on truncated Besov 
norm 

Lemma 13.3 For any 1 < p < 2, there exists C > 0 s.t. for any f £ Lj oc fl <S'(R d ) 
and any A'eN 

K 

Y / 2 sk \\f k \\ L P<CVK\\f\\wv 

k =0 
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Proof By a simple Cauchy-Schwartz estimate 


K 



1/2 


J22 sk \\h\\L* < Vk [J2 22sk WfxWh 



which concludes by Prop. 113.21 


□ 
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